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Preface 



Experimental results constituting the basis for this dis- 
sertation have been published in twelve articles [1-12] and 
eleven conference proceedings [13-23]. The research has 
been realized at the cooler synchrotron COSY by means 
of the COSY-11 facility [24], which was build to a large 
extent at the Jagellonian University. 

Ideas of experiments connected with this treatise, elab- 
orated and formulated as experimental proposals [25-31], 
have been positively judged and approved for realization 
by the Programme Advisory Committee of the COSY ac- 
celerator, and the measurements have been carried out by 
the COSY-11 collaboration. 

To render the reading of this dissertation more 
pleasant I have attempted to adduce at the beginnig of 
each section a few sentences from books of philosophers 
and scientists whom I immensely admire. 
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1. Introduction 



It is extremely advantageous to be able to bring 
a number of investigations under a formula of a 
single problem. For in this manner, we not only 
facilitate our own labour, inasmuch as we define it 
clearly to ourselves, but also render it more easy 
for others to decide whether we have done justice 
to our undertaking [32]. 

Immanuel Kant 

Pseudoscalar mesons (tt, K, rj, rj') constitute the lightest (basic) nonet of 
particles built out of quarks and antiquarks, and the investigation of their 
interaction with nucleons is one of the key issues in hadron physics. 

The interaction of vr and K mesons with nucleons has been deduced 
from experiments realized by means of charged pion or kaon beams. Such 
experiments are, however, generally not feasible in the case of flavour 
neutral mesons due to their too short lifetime. Therefore, although r] and 
T]' mesons were discovered over fourty years ago [33-35] their hadronic 
interaction with nucleons has not been established. The scattering length 
- the very basic quantity describing the low energy interaction potential - 
in the case of the rj meson is poorly estimated and in the case of the r]' 
meson it is entirely unknown. Estimated values of the real part of the 
proton-r^ scattering length varies from 0.20 fm to 1.05 fm depending on 
the approach employed for its determination [36-39]. 

It is the primordial purpose of this treatise to present methods devel- 
oped in order to study these interactions, to explain the applied experi- 
mental techniques and to discuss the progress achieved. We will address 
also the issue of the quark-gluon stucture of these mesons since in analogy 
to the connection between electromagnetic and Van der Vaals potentials, 
we may perceive the hadronic force as a residuum of the strong interaction 
that occurs between quarks and gluons - the constituents of hadrons. 

In the quark model the rj and rj' mesons are regarded as the mixture of 
the singlet (?7i = (mm + ss)) and octet (773 = -^{uu + dd — 2ss)) 
states of the SU(3)-flavour pseudoscalar nonet. A small mixing angle 
Q — — 15.5° [40] imphes that the percentage amount of strange and 
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non-strange quarkonium in both mesons is almost the same: 

77 = 0.77^{uu + dd) -OMss, rj' = 0. 63 ^ {uu + dd) + 0.77 ss. 
V 2 V 2 

Although, the r] and r]' mesons possess similar quark structure, their 
physical properties are unexpectedly difTerent. Let us adduce three ex- 
amples. The ?7'(958) meson mass is almost two times larger than the 
mass of the ri{5A7) meson. The branching ratios for the decays of B and 
Ds mesons into the rj' meson exceed significantly those into the 7] meson 
and the standard model predictions [41,42]. There exist excited states 
of nuclcons which decay via emission of the rj meson, yet none of the 
observed baryon resonances decays via the emission of the rj' meson [43] . 
Thus, it is natural to expect that the production mechanisms of these 
mesons in the collision of nucleons as well as their hadronic interactions 
with nucleons should also differ from each other. 

In spite of the fact that the strange and non-strange quarkonium con- 
tent of 1] and 1]' mesons is similar, the r] meson remains predominantly 
the SU(3)-flavor octet and the rj' meson the SU(3)-flavor singlet. This 
indicates that the 77' meson may be mixed with pure gluon states to 
a much larger extent than the rj and all other pseudoscalar and vector 
mesons. The anomalously large mass of the 1]' meson [44] and the anal- 
ysis of the decays of pseudoscalar, vector, and J/ip mesons indicate [45], 
that the rj' meson is indeed a mixture of quarkonium and gluonium. The 
gluonic admixture of 77' will influence the r^'-nucleon interaction via the 
U(l) anomaly [44]. The range of the glue induced r^'-nucleon interaction, 
if determined by the two-gluon effective potential, would be in the or- 
der of 0.3 fm [12]. This range is large enough to be important in the 
threshold production of the rj' meson e.g. via the pp — > pprj' reaction 
which occurs at distances of the colliding nucleons in the order of 0.2 fm. 
At such small distances the quark-gluon degrees of freedom may play a 
significant role in the production dynamics of the rj and rj' mesons. It 
is essential to realize that the investigation of the interaction between 
hadrons is inseparably connected with the study of their structure and 
the dynamics of the processes in which they are created. Therefore, in 
this treatise we present our study of the interaction between protons 
and T) or rj' mesons in strict connection with the investigation of their 
structure and production mechanism. 

The fact that fourty years after the discovery of the 77 and 7/' mesons, 
their interaction with nucleons remains so weakly established, in spite 
of its crucial importance for hadron physics, indicates that it is rather 
challenging to conduct such research. 
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From the methodological point of view, the first difficulty needed to 
be solved was to find a way of how the interaction with nucleons of such 
short-lived objects like rj and t]' mesons can manifest itself in a mea- 
surable manner. For, even when moving close to the velocity of light, 
both these mesons disintegrate on the average within a distance of tens 
of femtometers rendering their direct detection impossible. Not to men- 
tion, that it is completely unfeasible to accomplish out of them a beam or 
a target. However, the creation of these mesons in the vicinity of nucleons 
with low relative velocities appeared to us to be a promising tool which 
could facilitate the study of meson-nucleon interactions. For this purpose 
the best suited is the production close to the kinematical threshold or in 
kinematical regions where the outgoing particles indeed possess low rel- 
ative velocities and hence remain in distances of few Fermi long enough 
to experience the hadronic interaction. Investigations were started with 
measurements of the excitation function for the pp — > pprj and pp — * pprj' 
reactions, searching for its statistically significant distortion caused by 
the interacting ejectiles. We scanned the energy region close to the kine- 
matical threshold where the effect was supposed to be most significant 
and where the interpretation of the results is not burdened by the com- 
plexity resulting from contributions of many angular momentum states 
of initial and final systems. Indeed, we found out that the interaction 
between protons augments the total production cross section gradually 
with decreasing volume of the phase space. At the threshold the enhance- 
ment is most pronounced and corresponds to an enlargement of the cross 
section by more than an order of magnitude and it vanishes where the 
kinetic energy shared by ejectiles in the centre-of-mass system exceeds 
~ 100 MeV. It was most fascinating to see that the interaction of the 
rj meson with protons enhanced the cross section by a further factor of 
two. As it will be demonstrated in this dissertation, we were able to 
discern this from the effect caused by the overwhelming proton-proton 
interaction in a model independent manner. This was possible by making 
a comparison of the excitation functions of the pp —>■ pprj and pp —>■ ppn^ 
reactions. 

We would like to stress that the revealed phenomenon - viz the in- 
teraction of the created object with the outgoing nucleons affects the 
probability of its very production - is of a purely quantum mechanical 
nature. Our understanding of this peculiar effect will be presented in 
chapter 121 However, it is by far easier to imagine that the interaction be- 
tween ejectiles modifies the distributions of their relative momenta. This 
appears much more intelligible even on the analogy of the well known 
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effect that the momentum distributions of (3 particles are different for 
positrons and for electrons due to their Coulomb interaction with the 
residual nucleus. Therefore seeking for another measurable manifestation 
of the interaction between the t] meson and the proton we performed a 
measurement of the pp — > ppr] reaction with a statistics which permit- 
ted to derive distributions of a complete set of differential cross sections. 
The determined spectra of the invariant mass distributions for the two- 
particle subsystems of the ppr) final state revealed strong deviations from 
the predictions based on the assumption that all kinematically allowed 
momentum configurations of ejectiles are equally probable. Deviations 
at low relative proton-proton momenta can be explained satisfactorily 
well by the hadronic interaction between protons. Yet, a pronounced 
discrepancy observed in the region of large relative momentum between 
protons is by far larger than expected from an influence of the //-proton 
hadronic potential. This intricate observation will be examined in detail 
in section 14.51 Its explanation appears to be quite a challenge due to the 
need for the description of a three-body system subjected to the complex 
hadronic potentials. Later on in section 14.71 we will outline our ongoing 
investigations aiming at the determination of the hadronic interaction 
occuring between the rj' meson and the proton. In this case, contrary to 
the 7] meson, the excitation function for the pp — ^ ppr]' reaction did not 
reveal any statistically significant signal which could have been assigned 
to the ?7'-proton interaction. At present we are continuing these stud- 
ies by determining distributions of the invariant masses for two-particle 
subsystems of the pprj' final state which might deliver the first ever ex- 
perimental evidence for this still unknown potential. Such investigations 
constitute a much more difficult experimental task compared to the study 
of the rj meson. This is mostly due to the fact that the cross section for 
the pp —>■ pprj' reaction is by about a factor of thirty smaller than the 
one of the pp —>■ pprj reaction at the corresponding excess energy. 

Thorough theoretical investigations have shown that in order to un- 
derstand how the energy of motion is converted into a meson, it is of 
utmost importance to take into account also an interaction between nu- 
cleons before the very collision. We will demonstrate in section that 
the interaction between nucleons in the initial channel reduces the total 
cross section drastically. For the pp pprj reaction it is decreased by a 
factor of four and in the case of the pp pprj' reaction by a factor of 
three! 

As an essence of the above part of the introduction we would like to 
emphasize that the inferences about the hadronic interaction of r] and r]' 
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with nucleons will be based on comparisons of both the differential cross 
sections and the close-to-threshold excitation functions for the pp ppr] 
and pp ppr]' reactions with predictions based on the assumption that 
the kinematically available phase space is homogeneously populated. 

It cannot be, however, a priori excluded that possible deviations of 
the experimental distributions from the phase space governed calcula- 
tions may also be due to other physical effects caused by the specific 
production mechanism. Therefore, the discussion of the dynamics of the 
considered reactions, embodied in chapter El constitutes a part of this 
thesis. We will show that at threshold the entire production dynamics 
manifests itself in a value of a single constant which determines the mag- 
nitude of the total cross section. Consequently, the measurements of one 
reaction channel (eg. pp ppX) is by far not sufficient for a full under- 
standing of the reaction mechanism, since having only one observable at 
disposal it is impossible to establish contributions from many plausible 
production currents (as an example let us mention mesonic, nucleonic, 
or resonance currents). Therefore, an exploration of isospin and spin 
degrees of freedom is mandatory if we like to corroborate or refute the 
hypotheses proposed. In addition, the question of the reaction mecha- 
nism is connected with the very basic problem of whether hadronic or 
quark-gluon degrees of freedom are more appropriate for the description 
of phenomena occuring in the energy regime relevant to our studies. We 
address this issue in chapter El which comprises the description of cur- 
rently considered production models and gives account of present and 
planned future experimental studies of spin and isospin observables. For 
example, we will discuss to what extent one can infer a gluonic compo- 
nent of the 1]' wave function from a comparison of its creation in collisions 
of nucleons with different isospin combination. 

Consideration of the quasi-free meson creation on virtual nucleons 
bound inside a nucleus led us to the idea of studying the dependence of 
the meson production cross section on the virtuality of colliding baryons. 
The idea will be elucidated in section 15.71 and a method for the mea- 
surement of the close-to-threshold production of mesons in the quasi-free 
proton-neutron and even neutron-neutron collisions will be demonstrated 
in chapter ini 

From the experimental point of view, measuring the energy depen- 
dence of the total cross section for the production of mesons close to the 
kinematical threshold was quite a challenge. First of all, because these 
cross sections are by more than seven orders of magnitude smaller com- 
pared to the total yield of pp reactions, and secondly because they vary 
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by few orders of magnitude in a few MeV range of the excess energy. 
The afore mentioned relations are demonstrated in figure 11.11 It shows 
recently determined excitation functions for the close-to-threshold pro- 
duction of some of the pseudoscalar mesons, and in particular it shows 
the data for the rj and rj' mesons, which constitute an experimental basis 
of this treatise. 
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Figure 1.1: Close-to-threshold 
cross section for the proton-proton 
interaction leading to the produc- 
tion of mesons whose wave function 
comprises a significant amount of 
strangeness. For comparison, the 
total cross section of proton-proton 
collisions is also shown. The 
filled points depict data taken at 
COSY [3,10,11,18,46-49], and open 
symbols show results from other 
laboratories [50-60]. 



Conclusive inferences concerning the interaction between ejectiles could 
be drawn only due to the unprecedented experimental precision obtained 
using proton beams available at storage ring facilities and in particu- 
lar at the cooler synchrotron COSY. These beams are characterized by 
low emittance, small momentum spread, and well defined absolute mo- 
mentum. Experiments which will be discussed in this dissertation have 
been carried out by means of the COSY-11 facility [24,61] installed at 
COSY [62]. Details of the experimental method are elucidated in chap- 
ter IHl In short, the technique is based on the determination of the four 
momentum vectors of colliding and outgoing nucleons, and on the deriva- 
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tion of the mass of the produced meson or group of mesons employing 
the principle of energy and momentum conservation. The registered par- 
ticles are identified by independent measurements of their velocity and 
momentum, where the former is determined from the time-of-flight be- 
tween scintillation detectors and the latter is deduced from the curvature 
of the particle trajectories in a magnetic field. Figure IT!^ presents an ex- 
ample of the mass spectrum determined for a system X produced via 
the pp ppX reaction at an excess energy of 1.5 MeV with respect 
to the kinematical threshold for the production of the r]' meson. A sig- 
nal from the creation of the rj' meson can be well recognized and its 
width of 0.7 MeV (FWHM) indicates the precision achieved using the 
stochastically cooled proton beam of COSY and the COSY-11 detection 
setup. When examining the error bars, it is evident that if the accuracy 
was twice worse the signal would be disputable. The method invented 
to monitor the beam geometry at internal target facilities [8,23], which 
is described in section 16.21 enabled us to improve the accuracy of mass 
determination by correcting for variations of the avarage beam character- 
istics. On the cognitive side, its application led to the understanding of 
to what extent the stochastic cooling improves the quality of the beam. 



Figure 1.2: Missing mass spectrum for 
the pp ppX reaction measured by 
means of the COSY-11 detection facil- 
ity at the cooler synchrotron COSY. The 
measurement was performed at a beam 
momentum corresponding to an excess 
energy of 1.5 MeV for the pp — > pp-q' re- 
action. More details concerning the eval- 
uation of the data will be explained in 
section 16.11 It is worth noting that the 
attained experimental mass resolution is 
comparable with the natural width of the 

r?' meson {V„, = 0.202 MeV [431). 
0.952 0.955 0.958 0.961 0.964 ' ' 

missing mass [ GeV/c^ ] 

Necessary corrections for a non-ideal detection geometry were per- 
formed in an utterly model independent manner. To arrive at this aim 
we expressed the experimental acceptance as a function of a complete set 
of mutually orthogonal variables describing the reaction, and to facilitate 
calculations we digitized the phase space volume into a few thousand par- 
titions in which the acceptance could be safely regarded as constant. In 
the analysis, we assigned to each measured event a weight equal to the 
inverse of the acceptance according to the phase space bin to which it 
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belonged. This allowed us to elaborate the differential distributions in- 
dependently of the reaction model used in simulations. Employing the 
missing mass technique it is in principle impossible to discern between 
the multi-pion and t] or t]' meson production on the event-by-event basis. 
Therefore, to derive the differential distributions we divide the range of 
the studied variable into bins with a width corresponding to the achieved 
experimental resolution and separate the multi-pion background from 
the 7] or 1]' meson signal for each bin. The effort was undertaken in or- 
der to determine background- and model-free distributions of differential 
cross sections which in turn enable to test various hypotheses without an 
obscurity due to unnecessary assumptions. Hence we can benefit thor- 
oughly from simplifications gained due to the threshold kinematics. The 
most crucial facilitation and attractiveness when interpreting the meson 
production at the vicinity of the threshold is the fact that the relative 
angular momenta larger than I = Oh play no role due to the short range 
of the strong interaction and small relative momenta of the produced 
particles. It can be inferred from parity and angular momentum conser- 
vation laws that the production of a nucleon-nucleon-meson system (for 
pseudoscalar or vector mesons) in relative S-waves may only occur if the 
nucleons collide in P-wave. Thus at threshold the transition P — Ss is 
the only possible one, with capital letters denoting the angular momen- 
tum between nucleons and the small letter describing the meson angular 
momentum with respect to the pair of nucleons. Section will be de- 
voted to the circumstantial discussion of this issue. At this stage we 
would only like to note, that the dominance of the P — > Ss transition 
at threshold will not be unquestioningly assumed but rather the energy 
range of its applicability - which varies strongly with the mass of the 
created meson - will be established in section 14.11 

The present work is divided into seven chapters. After this intro- 
duction we will elucidate our understanding of the manifestation of the 
hadronic interaction via the excitation function and the phase space 
abundance of the studied reactions. Further, in the third chapter we 
will introduce the definitions of observables used throughout this work. 
The fourth chapter comprises a discussion of the experimental evidence 
concerning the interaction between the pseudoscalar isosinglet mesons [r], 
rj') and the proton. In the fifth chapter we explain the reaction dynam- 
ics and discuss the relevance for the production process of hadronic and 
quark-gluon degrees of freedom. The final conclusion is preceded by the 
sixth chapter which contains the detailed description of the experimental 
method. 



2. Manifestation of hadronic 
interactions 

... what it can mean to ask what is the probability of past 
event given a future event! There is no essential problem 

with this, however. Imagine the entire history of the uni- 
verse mapped out in space-time. To find the probability of 
p occurring, given that q happens, we imagine examining 
all occurrences of q and counting up the fraction of these 
which are accompanied by p. This is the required probabil- 
ity. It does not matter whether q is the kind of event that 
would normally occur later or earlier in time than p [63]. 

Roger Penrose 

The interaction of liadrons - being tiie reflection of tlie strong force 
acting between their constituents - provides information about their 
structure and the strong interaction itself. In the frame of the opti- 
cal potential model the hadronic interaction can be expressed in terms of 
phase-shifts, which in the zero energy limit are described by the scatter- 
ing length and effective range parameters. These are quite well known 
for the low-energy nucleon-nucleon interaction [64,65], yet they are still 
poorly established in the case of meson-nucleon or meson-meson inter- 
actions. This is partly due to the absorption of mesons when scattering 
on a baryon. To account for this effect the scattering length becomes a 
complex quantity where the imaginary part - e.g. in case of the nucleon-77 
interaction - describes the rjN nN and ^^A^ multi-T: N processes. 
Moreover, the short life-time of all neutral ground state mesons pro- 
hibits their direct utilization as secondary beams and therefore the study 
of their interaction with hadrons is accessible only via their influence 
on the cross section of the reactions in which they were produced (eg. 
A^A^ NN Meson). When created close to the kinematical threshold 
with the relative kinetic energy being in the order of a few MeV, the 
flnal state particles remain much longer in the range of the strong in- 
teraction than the typical life-time of N* or A baryon resonances with 
10~^^ s, and hence they can easily experience a mutual interaction before 
escaping the area of an influence of the hadronic force. This - as in- 
troduced by Watson [66] - flnal state interaction (FSI) may signiflcantly 
modify both the original distributions of relative momenta of the outgo- 
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ing reaction products and the magnitude of the production cross section. 
Although it is easy to understand that the FSI changes the distributions 
of the differential cross sections it is rather difficult to cope with the in- 
fluences on the magnitude of the total reaction rate, since one tends to 
separate the primary production from the final state interaction in space 
as well as in time [66]. The influence of the FSI on the absolute scale 
of the cross section is a pure quantum mechanical effect and our classi- 
cal imagination must inevitably fail in this respect. However even when 
considering the primary production as a separate process it is well worth 
trying to understand the phenomenon qualitatively. If there were no fi- 
nal state interactions the total cross section would be fully determined 
by the kinematically available phase space volume Vps, where each inter- 
val is populated with a probability governed by the primary production 
amplitude only: 



The approximation in the equation results from the assumption that 
|Mp = constant in a few MeV range above the production threshold [67- 
69]. F denotes the flux factor of the colliding particles. 

In the classical picture we might imagine that the reaction particles are 
created together with their appropriate force field and when escaping the 
interaction region they acquire a potential energy increasing or decreasing 
their kinetic energy depending whether the interaction is repulsive or 
attractive. For an attractive interaction they could be created also in 
those phase space partitions which are not available for non-interacting 
particles and subsequently be "pulled down" to the energetically allowed 
regions by final state interaction. The temporary growth of the primary 
production phase space would then increase the reaction rate. Contrary, 
in the case of a repulsive interaction the particles must be produced 
in the lower phase space volume since leaving the interaction area they 
will acquire additional kinetic energy. The reduction of the total cross 
section, for example in case of the repulsive Coulomb force, can easily 
be understood when considering the production in a coordinate space 
of point-like objects. Here - in contrast to non-interacting particles - 
a strongly repulsive object cannot be produced at appropriately small 
distances since their later acceleration would violate energy conservation 
and thus the space available to primary production is reduced. 

Another more commonly anticipated presupposition of changes of the 
reaction yield due to the final state interaction is based on the time 




const. ■ Vr, 



(2.1) 
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invariance principle. It asserts that the cross section for the production 
of repulsive particles is smaller in comparison to non-interacting particles, 
all other features being equal, because in the time inversed process the 
repulsive interaction hinders the slow particles from reaching distances 
at which a particular reaction may occur. 
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Figure 2.1: (left) The solid line indicates the kinematically available area for the ppr} sys- 
tem with a total centre-of-mass energy — 2433.8 MeV (excess energy Q ~ 10 MeV). 
The dashed line depicts the range assuming a reduction of the proton and ?7-meson 
masses by 2 MeV. The phase space volume results from an integral over the closed area: 
Vps = / / dnipj^ri d^p2V (Terms of (27r)" are skipped here and will be included in the 
flux factor F according to the convention introduced by Byckling and Kajantie [70]). 
(right) Distribution of the phase space for the ppr] system modified by the proton-proton 
interaction and calculated for an excess energy of Q = 16 MeV. The area of the squares is 
proportional to the number of entries and is shown in a linear scale. 

One can also argue, that relativistically the primary mechanism cre- 
ates the particles off the mass shell and subsequently they are lifted onto 
the mass shell by the final state interaction. The solid line in figure 12.11 
(left) depicts the boundary of the Dalitz plot in case of the pp ppr] 
reaction calculated at the total centre-of-mass energy ^s = 2433.8 MeV 
exceeding the threshold energy by 10 MeV. The area surrounded by that 
curve is a direct measure of the kinematically available phase space vol- 
ume. The dashed line shows the corresponding plot at the moment of 
the primary creation if the mass of each particle was reduced by 2 MeV, 
demonstrating that now the available phase space grows significantly. 
Indeed, as shall be inferred from the experimental results presented in 
subsequent sections, at excess energies of a few MeV above threshold, the 
mutual interaction among the outgoing particles enhances drastically - 
by more than an order of magnitude - the total cross section and modifies 
appreciably the occupation of the phase space. Figure 12.11 (right) indi- 
cates the phase space distribution expected for the pprj system at an 
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excess energy of Q = y/s — 2 nip — = 16 MeV, assuming a homoge- 
neous primary production and taking into account the S-wave interaction 
between the protons. The proton-proton FSI modifies the homogeneous 
Dahtz plot distribution of "non-interacting particles" , enhancing its pop- 
ulation at a region where the protons have small relative momenta. The 
interaction of the proton-?] system would manifest itself at low invariant 
masses m^^ corresponding to small relative momenta between the proton 
and the r] meson. Such effects observed in the experiments are presented 
in section 

The effect of the nucleon-nucleon FSI diminishes with increasing ex- 
cess energy since it significantly infiuences only that partition of the phase 
space at which the nucleons have small relative momenta. While this 
fraction stays constant, the full volume of the phase space grows rapidly: 
An increase in the excess energy from Q = 0.5 MeV to Q = 30 MeV cor- 
responds to a growth of Vps by more than three orders of magnitude. As 
a result the S-wave nucleon-nucleon FSI is of less importance for higher 
excess energies where it affects a small fraction of the available phase 
space volume only. A more quantitative discussion about the influence 
of the FSI and ISI (initial state interaction) effects at close-to-threshold 
production cross sections will be presented in section I^!T1 

In this section we would only like to state that an interaction of the 
nucleons in the entrance channel, similarly to that among ejectiles, in- 
fluences the production process noticeably [71]. For example, for rj pro- 
duction in the pp ppt] reaction the initial state interaction reduces the 
total cross section by a factor of about 3-5 [71, 72] due to the repulsive 
proton-proton ^Po-wave potential. This factor remains rather constant 
in the range of a few tens of MeV [72] and hence, does not influence the 
energy dependence of the total cross section of the meson production, 
which is mostly determined by the final state interaction. In particu- 
lar, the dominant S-wave nucleon-nucleon final state interaction is by far 
stronger than any of the low-energy meson-nucleon ones. The reduction 
of the cross section by the initial state interaction we may intuitively 
conceive as a loss of the initial fiux in favour of elastic, and a plethora of 
other - than the studied one - inelastic reactions, which may occur when 
the colliding nucleons approach each other. 

The above considerations point out that the hadronic interaction be- 
tween nucleons and the short-lived mesons like rj and t]' may be investi- 
gated by the observation of the excitation function for the production of 
these mesons off nucleons, as well as by the studies of the phase space 
abundance. 
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Hermogenes. Suppose that we make Socrates a 

party to the argument? 
Cratylus. // you please. 

Socrates. Then the argument would lead us to infer 
the names ought to he given according to a natural 
process, and with a proper instrument, and not at 
our pleasure: in this and no other way shall we 
name with success. 

Hermogenes. / agree. 

Socrates. But again, that which has to he cut has to 

he cut with something? 
Hermogenes. Yes. 

Socrates. And that which has to he woven or pierced 

has to he woven or pierced with something? 
Hermogenes. Certainly. 

Socrates. And that which has to he named has to he 

named with something? 
Hermogenes. True [73]. 

Plato 

3.1 Cross section 

Why is it that, in most cases, the definitions which 
satisfy scientists mean nothing at all to children? [74] 

Henri Poincare 

Investigations of the production of mesons and their interactions with 
nucleons are based on measurements determining the total and differen- 
tial production cross sections and their dependence on the energy of the 
interacting nucleons. Therefore, to enable a quantitative discussion on 
the mechanisms leading to the transformation of the cncrgy-of-motion of 
nucleons into matter in the form of mesons let us recall the formula of 
the reaction cross section. In the case of the NN — > NN Meson process 
- with the four-momenta of the colliding nucleons denoted by Pa and 
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and with n = 3 particles in the exit channel - it reads: 




where {Mab^usl"^ denotes the square of the Lorentz- invariant spin aver- 
aged matrix element describing the probability to create two nucleons and 
a meson with four-momenta of = (Ej,Pj) and i — l..n, respectively. 
The energy and momentum conservation as well as the on-shellness of the 
created particles is ensured by the Dirac-S and the Heaviside-G functions. 
The formula holds also for n 7^ 3. 

The total cross section is then defined as an integral of the probabilities 
to populate a given phase space interval over the whole kinematically 
available range - determined by energy and momentum conservation - 
normalized to the flux factor F of the colliding nucleons. In the case of 
non-interacting flnal state particles the matrix element close to threshold 
|-^a6^i23p is nearly constant [67-69] and hence the allowed phase space 
volume, Vps, is the decisive quantity which governs the growth of the 
cross section with increasing excess energy Q. The latter - deflned as 
the total kinetic energy - is shared among the outgoing particles in the 
reaction centre-of-mass frame ^: 

n 

Q^v^-J^m,, (3.2) 

i=l 

where s = |Pa + Pf,p = | X;r=i P denotes the square of the total centre- 
of-mass energy. Exactly at threshold, where the particles' kinetic energy 
in the centre-of-mass system equals zero (Q = OMeV), the total reac- 
tion energy ^/s amounts to: ^/sth = nij. Before writing explicitly 
the formula for Vpg let us introduce the kinematical triangle function A 
defined as [70]: 

X{x, y, z) — -\- -\- — 2xy — 2yz — 2zx. (3-3) 

It enables to formulate many useful kinematical variables in a very com- 
pact and Lorentz-invariant form. In particular, the momenta of particles 



^ Traditionally by the centre-of-mass system we understand a frame in which the momenta of all 
particles add to zero, called sometimes more explicitly centre-of-momentum frame. 
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i and j in their centre-of-mass frame may be expressed as follows: 



A(sij,m2,m|) 

where Sij = |Pj + P^p stands for the square of the invariant mass of 
the ij system considered as one quasi-particle. The above relation gives 
an expression for the flux factor F from equation p.lj) in terms of the 
colliding masses of the nucleons and the total energy s only, namely: 



F = 4 {211 f^-' pl = 2 (27r)3"-4 ^A(s, m^, m^), (3.5) 

where we have chosen the convention introduced by Byckling and Ka- 
jantie [70] and included the (27r) factors for the phase space (27r)^" and 
for the matrix element (27r)~^ into the definition of F. It is important to 
note that at threshold, for an excess energy range of a few tens of MeV, 
the small fractional changes of total energy = — ; — —) causes 

o '=>•' y mi + m2 + ni3 + Q ^ 

weak variations of the flux factor and influences the shape of the energy 
dependence of the total cross section only slightly. 

For a two-particle final state ab — >■ 12 (for instance for the reaction 
pn — > drj) the phase space integral defined in equation (j3.ip reduces to 
Vps := / dVps = -^Pi = ^\/A(s, mf , m|) and the total cross section 
for such reactions (when neglecting variations due to dynamical effects 
{{Mab^ul = const.)) should increase linearly with the centre-of-mass 
momentum of the produced meson in the vicinity of the reaction thresh- 
old. The total cross section can be expressed analytically as a function 
of the masses of the particles participating in the reaction and the square 
of the total reaction energy s: 

Vps const pI 1 a/A(s, ml,ml) 

o'ab-^u = const ■ = = . 3.6 

F IGttsp: 167rs^A(s,m2,m^) ^ ^ 

Near threshold, at a given excess energy Q = ^/s — mi — m2, the emis- 
sion of the reaction products in the centre-of-mass frame is isotropic and 
the whole dynamics of the process manifests itself in the absolute value 
of the transition matrix element {Mab^ul- The underlying production 
mechanisms can also be extracted from the deviations of the total cross 
section energy dependence following the prediction of relation ()3.fij) . A 
visualization of possible differences in the dynamics of the production of 
various mesons can be obtained by comparing the total cross sections of 
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the studied reactions at the same value of the phase space volume nor- 
malized to the flux factor (Vps/F). It is important to recognize that the 
comparison at the same centre-of-mass meson momentum (P2) gives di- 
rect information about the amplitude differences only if the production 
of mesons with the same masses is concerned (for example pp d7r~^ 
and nn dii~), yet the reactions pn dt] and pn — >■ dii^ have by far 
different Vps at the corresponding P2. 

Prior to the comparative analysis of the transition amplitudes for the 
pp —>■ ppr], pp —>■ pprj', and pp — > ppir^ reactions, a corresponding formula 
of the phase space volume Vps for the three particle final state will be 
introduced in section 



3.2 Partial waves — selection rules 

Not for nothing do we call the laws of nature 'laws ': 
the more they prohibit the more they say [75]. 

Karl Raimund Popper 



If one is interested in the decomposition of the production amplitude 
according to the angular momenta of the final state particles then indeed 
the appropriate variable for a qualitative comparison even for particles 
with different masses is the meson momentum in the reaction centre-of- 
mass frame. Correspondingly, for the more than two-body final state, 
the adequate variable is the maximum meson momentum, since it is 
directly connected with the maximum angular momentum by the inter- 
action range. In the case of a three-body exit channel {ab 123) the 
meson (ms) possesses maximum momentum (q^aa;) when the remaining 
two particles are at rest relative to each other. Hence, employing defini- 
tion ()3.4p one obtains: 

_ VA(s, (mi + m2)^m|) 

Contrary to the two-body final state, at a fixed excess energy the 
dynamics of the meson production associated with two or more particles 
reflects itself not only in the absolute value of the square of the matrix 
element but also in distributions of variables determining the final state 
kinematics. Usually, in non-relativistic calculations of the total cross 
section, one takes the Jacobi momenta, choosing as independent variables 
the 'q'-meson momentum in the reaction centre-of-mass frame and 'k'- 
momentum of either nucleon in the rest frame of the nucleon-nucleon 
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subsystem. By means of the A function they can be expressed as: 



q= V^S^l^, and k=V^Sl^M;M), (3.8) 
2 vs 2yM2 



with Si2 = |Pi + P2P denoting the square of the invariant mass of the 
nucleon-nucleon subsystem. In a non-relativistic approximation the ex- 
pression of the total cross section defined by formula (jH.lll for a meson 
production reaction in nucleon-nucleon interactions of the type A^A^ — >■ 
A^A^ Meson simplifies to: 

(Tcx / kq2|M,fe^i23|'f^q. (3.9) 
Jo 

Denoting by L and / the relative angular momenta of the nucleon-nucleon 
pair and of the meson relative to the A^A^ system, respectively, and 
approximating the final state particles by plane waves (case of non- 
interacting objects), whose radial parts ipi{q, r) are given by the spherical 
Bessel functions: 

V^/(q,r)(XjKqr)^-^°^^|^, (3.10) 

an expansion of the transition amplitude jM^,;]^ for the LI partial wave 
around qr = leads to 

|Mi,poc|M°,pk2^q2', (3.11) 

where, M^i denotes the matrix elements responsible for the primary pro- 
duction with the angular momenta L and /. Thus, the total cross section 
can be expressed as the following sum of the partial cross sections an- 

c7 = J2^u^Yl dq (3.12) 

L,l L,l 

Furthermore, if the strength of the the primary production partial 
amplitudes |M°;| were constant over the phase space then the energy 
dependence of the partial cross sections, obtained by solving the above 
equation, would be given by: 

^L.ocq^^+/+^ocr^^.r^+^ (3.13) 

where rjM = q^^x/^M with irm denoting the mass of the created me- 
son ^. Thus - at threshold - for the Ss partial wave the cross section for 



^ In former works [76, 77] dealing only with pions this parameter is denoted by 7], here in order 
to avoid ambiguities with the abbreviation for the eta-meson, we introduce an additional suffix M. 
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the A^A^ NN Meson reaction should increase with the fourth power of 
rjM- The dimensionless parameter tjm was introduced by Rosenfeld [76] 
as a variable allowing for qualitative estimations of the final state partial 
waves involved in pion production. He argued that if the phenomenon 
of pion production takes place at a characteristic distance R from the 
centre of the collision, with R in the order of h/Mc, then the angular 
momentum of the produced meson is equal to / = Rq = hq/Mc. Hence, 
r]M denotes the classically calculated maximum angular momentum of 
the meson relative to the centre-of-reaction. The same arguments one 
finds in the work of Gell-Mann and Watson [77], where the authors do 
not expect the range of interaction to be larger than the Compton wave- 
length of the produced meson [h/Mc). However, it is rather a momentum 
transfer Ap between the colliding nucleons which determines the distance 
to the centre of the collision R ^ h/Ap at which the production occurs. 
Based on indications from data, it is emphasized in the original article of 
Rosenfeld [76] that R is slightly less than h/{2Mc), which is numerically 
close to the value of h/Ap. Directly at threshold, where all ejectiles are 
at rest in the centre-of-mass frame, Ap is equal to the centre-of-mass mo- 
mentum of the interacting nucleons and hence, exploring equation ()3.4|) 
it can be expressed as: 



Aptf, = —— \ mpmx + — . 3.14 

2 V 4 



Though the present considerations are limited to the spin averaged pro- 
duction only, it is worth noting that very close-to-threshold - due to 
the conservation laws and the Pauli excluding principle - for many reac- 
tions there is only one possible angular momentum and spin orientation 
for the incoming and outgoing particles. The Pauli principle for the 
nucleon-nucleon system implies that 

where L, S, and T denote angular momentum, spin, and isospin of the 
nucleon pair, respectively. For example if the nucleon-nucleon wave func- 
tion is symmetric in the configuration-space {L = 0) as well as in the 
spin-space (5 = 1) then it must be antisymmetric in the isospin-space 
(T = 0) to be totally antisymmetric. 

In the conventional notation [78, 79] the transition between angular 
momentum parts of the initial and final states of the NN — >■ NNX 
reactions are described in the following way: 

'S'+^Vj.-^''+'Lj,l (3.16) 
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where, superscript i indicates the initial state quantities, / denotes the 
angular momentum of the produced meson with respect to the pair of 
nucleons, and J stands for the total angular momentum of the N-N sys- 
tem. The values of angular momenta are commonly being expressed 
using the spectroscopic notation (L = S,P,D,... and / = s,p,d,...). The 
parity conservation implies that + L + 1 must be an even number if the 
parity of the produced meson is positive or odd if this parity is negative. 
Employing additionally the conservation of the total angular momentum 
and Pauli principle one can deduce that at threshold the A^A^ NNX 
reactions will be dominated by the transitions listed in table IHTTl 



Table 3.1: Partial wave transitions for the pp pp Meson and nn —> nn Meson reactions 
at threshold 



type 


meson 


spin and parity 


transition 


pseudoscalar 
vector 
scalar 


n,r],r]' 
p,uj,4> 
O'Oi fo 


0- 

1- 
0+ 


'Po ^ 'Sos 

^Pl 'SqS 



For example, the production of neutral mesons with negative parity 
- as pseudoscalar or vector mesons - may proceed in the proton-proton 
collision near threshold only via the transition between ^Pq and ^SqS 
partial waves. This means that only the collision of protons with rela- 
tive angular momentum equal to 1 ^ may lead to the production of such 
mesons. Moreover, an inspection of the corresponding Clebsch-Gordan 
coefficients reveals that it is four times more probable that the orienta- 
tion of spins of the colliding nucleons is parallel than anti-parallel. 
In the case of the production of neutral scalar mesons, the protons or 
neutrons must collide with anti-parallel spin orientations which remains 
unchanged after the reaction. These simple considerations imply that 
in close-to-threshold measurements with polarized beam and target one 
should see a drastic effect in the reaction yield depending whether the 
polarization of the reacting protons is parallel or anti-parallel. Indeed, in 
reality, strong differences in the yield for various combinations of beam 
and target polarization have been determined in the pioneering measure- 
ments of the reaction pp — > ppn^ at the lUCF facility [80] . 

Although close-to-threshold higher partial waves are not expected, 
detailed studies of meson-nucleon interaction require a careful determi- 
nation of their contributions, in order to avoid false assignment of the 
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observed signals. If we allow for the production of protons with the rela- 
tive angular momentum L = or 1 and confine the t] meson production 
to the / = only, then the pp pprj reaction may proceed via three 
possible transitions listed in table I3.21 



Table 3.2: Transitions for the pp — > ppMeson{0~) reaction in the case where I — and 
L = or 1. More comprehensive list of allowed transitions for the low partial waves in the 
A'^A'' NNX reaction can be found in references [81,82] 

'So ^ =^Pos 



It is worth to note that the conservation rules mentioned above forbid 
the production of the t] meson with / = 1 if the protons are in the S-wave 
(^Sqp). In other words, the r] meson can be produced in the p-wave only 
if L is larger than zero. This suggests that with an increase of the excess 
energy the higher partial waves in the ppr/ system should first appear in 
the proton-proton subsystem. In fact such an effect was already observed 
in experiments at the SATURN laboratory in the case of the (j) meson 
production via the pp ppcp reaction [83]. Since the rj production is 
dominated by the formation and deexcitation of the S-wave baryonic 
resonance (S'ii(1535)) this supposition becomes even more plausible. 

Another important feature of the pp ppX and also nn — >■ nnX 
reactions, pointed out in reference [82], is the fact that the interference 
terms between the transitions with odd and even values of the angular 
momentum L of the final state baryons are bound to vanish. This char- 
acteristic is due to the invariance of all observables under the exchange 
of identical nucleons in the final state. Thus, for instance, out of three 
possible interference terms for transitions listed in table 13.21 only one - 
namely the interference between the ^Pq s and ^P2 s final states - will 
contribute to the production process significantly simplifying an inter- 
pretation of the experimental data. In general, conservation of the basic 
quantum numbers leads - as derived in reference [81] - to the following 
selection rule for the A^A^ —>■ NNX reaction: 

(_^)(A5+AT) ^ (-1)', (3.17) 

where irx describes the intrinsic parity of meson X, AS denotes the 
change in the spin, and AT in isospin, between the initial and final 
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A^A^ systems. A possible non-zero angular momentum between the out- 
going particles can manifest itself in an unisotropic population of the 
appropriate angles. In particular, the angular distribution of the rj me- 
son observed in the overall centre-of-mass frame should reflect a possible 
non-zero value of /. However, for the study of the angular momentum L 
we need to use another reference system, since for a three body final state 
the beam axis is not a good reference direction to look for the angular 
distributions relevant for the relative angular momenta of the two parti- 
cles [84] . This is because for the fixed relative settings of the two protons 
the angle between the beam and the vector of relative protons momen- 
tum may acquire any value. Searching for an appropriate variable let us 
consider a two body scattering. In this case, the beam line, which is at 
the same time the line along which the centre-of-mass system is moving, 
constitutes a reference frame for the angular distributions. Therefore, by 
analogy to the two body system, an instructive reference axis for angular 
distributions in the proton-proton subsystem is now the momentum of 
the recoil t] meson, since the direction of that meson is identical to the di- 
rection of the movement of the proton-proton centre-of-mass subsystem. 
Figure ITT] visualises both above introduced angles 9* and 9** the latter 



center of mass system proton-proton rest frame 




Figure 3.1: Definition of tlie angles used in tfie text. Tlirougliout tlie wfiole work tlie angles 
in the overall centre-of-mass frame will be denoted by one asterisk, those in the two-particle 
subsystems by two asterisks, and the angle in laboratory system will be left without any 
superscript. 



being defined as the angle between the relative proton-proton momentum 
and the recoil particle (r/) seen from the di-proton rest system [84]. 
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3.3 Phase space population 

It is essential to science that its matter should be 
in a space, but the space in which it is cannot be 
exactly the space we see or feel [85]. 

Bertrand Russell 

Due to both the strong nucleon-nucleon low-energy interaction in the ^Sq 
state and the meson-nucleon forces in the exit channel, the assumption of 
non-interacting plane waves leading to equation ()3.13|) failed when con- 
fronted with the experimental data. However, these deviations offer the 
possibility to determine still poorly known nucleon-meson interactions. 
The interaction between particles depends on their relative momenta. 
Consequently, for investigations of final state interactions, more instruc- 
tive coordinates than the q and k momenta are the squared invariant 
masses of the two-body subsystems [84]. These are the coordinates of the 
Dalitz plot. In the original paper [86] Dalitz has proposed a representa- 
tion for the energy partitions of three bodies in an equilateral triangle 
whose sides are the axes of the centre-of-mass energies. He took advan- 
tage of the fact that the sum of distances from a point within the triangle 
to its sides is a constant equal to the height. Therefore, the height of the 
triangle measures the total energy -^s = -|- E2 -|- Eg and interior points 
- fulfilling four-momentum conservation - represent energy partitions. 
For a constant ^/s, due to the energy conservation, without loosing any 
information, it is enough to consider the projection on any of the Ej Em- 
planes. The linear relation between E* and Sjk {sjk = s -|- mf — 2 V^E*) 
allows to use sjk Sik or Ej Ej coordinates equivalently, with the following 
relation between the phase space intervals: dE* dE* = ^ dsjkdski- Ac- 
cording to Kilian's geometrical representation [84] a Dalitz plot lies on a 
plane in the three dimensional space (si2, S13, S23) orthogonal to the space 
diagonal. The plane including a Dalitz plot corresponding to a fixed total 
energy ^/s is then given by the following scalar product [84]: 

(1, 1, l)(si2, Si3, S23) = S12 + Si3 S23 = s -h m^ -h m2 + mg. (3.18) 

The second equality of equation ()3.18|1 means that there are only two 
independent invariant masses of the three subsystems and therefore a 
projection onto any of the Sjj Sjk planes still comprises the whole princi- 
pally accessible information about the final state interaction of the three- 
particle system. In the case of no dynamics whatsoever and the absence 
of any final state interaction the occupation of the Dalitz plot would be 
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fully homogeneous since the creation in each phase space interval would 
be equally probable. The final state interaction would then appear as a 
structure in that area. 



-2 




Spp[GeVV] c) Sp^[GeVV] d) 



Figure 3.2: Monte-Carlo simulations: (a) Phase-space distribution for the pp pp-q reac- 
tion at Q = 16 MeV modified by the proton-proton final state interaction, (b) The dotted 
line shows the projection of the pure phase space density distribution onto the Spp axis and 
the solid curve presents its modification by the proton-proton FSI. (c) The square of the 
scattering amplitude for the pp — » pp elastic scattering as a function of the proton-proton 
invariant mass in the range 2 mp < yspp < 2 mp + Q, calculated according to formula 14.61 of 
section (d) Phase-space density distribution modified by the proton-r; interaction, with 
a scattering length equal to ap,, — 0.7fm-|-i0.3fm. The proton- scattering amplitude has 
been calculated according to equation 1)4. 9|l . A detailed discussion of the nucleon-nucleon 
and nucleon-meson interaction will be presented in section 

Figure shows - for the example of the pp — ^ ppr] reaction - how the 
uniformly populated phase space density is modified by the S-wave (^Sq) 
interaction between outgoing protons. An enhancement in the range cor- 
responding to low relative momenta between protons is clearly visible. A 
steep decrease of the occupation density with increasing invariant mass of 
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the proton-proton subsystem is even better seen in figure presenting 
the projection of the phase space distribution onto the Spp axis indicated 
by an arrow in figure 13.2b . This is a direct reflection of the shape of the 
proton-proton (^Sq) partial wave amphtude shown in figure 13.2b . Fig- 
ure I3.2t i shows the Dahtz plot distribution simulated when switching off 
the proton-proton interaction but accounting for the interaction between 
the rj-meson and the proton. Due to the lower strength of this interaction 
the expected deviations from the uniform distributions are by about two 
orders of magnitude smaller, but still one recognizes a slight enhance- 
ment of the density in the range of low invariant masses of proton-?] 
subsystems. However, due to weak variations of the proton-?] scatter- 
ing amplitude the enhancement originating from the ?7-meson interaction 
with one proton is not separated from the ?7-meson interaction with the 
second proton. Therefore an overlapping of broad structures occurs. It 
is observed that the occupation density grows slowly with increasing Spp 
opposite to the effects caused by the S-wave proton-proton interaction, 
yet similar to the modifications expected for the P-wave one [87]. From 
the above example it is obvious that only in high statistics experiments 
signals from the meson-nucleon interaction can appear over the over- 
whelming nucleon-nucleon final state interaction. It is worth noting, 
however, that the Dalitz plot does not reflect any possible correlations 
between the entrance and exit channels [84]. 

The Dalitz plot representation allows also for a simple interpretation 
of the kinematically available phase space volume as an area of that plot. 
Namely, equation (|3.1|) becomes: 



where the limits of integrations defining the boundaries of the Dalitz plot 
can be expressed as [70]: 




(3.19) 



max / \ 2,2 
S23 (Sl2)=m2 + m3 



(si2 - s + m3)(si2 + mi - m?) - ^ A(si2, s, m§)A(si2, m^, m^) 



min / \ 2,2 

S23 (S12) = m2 + mg 



(si2 - s + m§)(si2 + mj - mi) + yA(si2, s, m§)A(si2, m|, mf ) 

2S12 
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Thus, the phase space volume kinematically available for the three-body 
final state can be written by means of only one integral: 



Vps= dVps = — I dSi2 I dS23 



4s 

(mi+m2)2 



(v^-m3)2 
TT^ f dSi2 

Is 



[ ^^/A(sl2,s,mi)A(sl2,mi,m?), (3.20) 

J Si2 * 

(mi+m2)2 

whose solution leads, in general, to elliptic functions [70]. However, in 
the nonrelativistic approximation it has the following closed form: 



TT^ ^mi m2 mg ^2 
2 (mi + m2 + m3)2 



Vps = Q^ (3.21^ 



where the substitution of the non-relativistic relation between 77^3 and Q 

_ m| + 2m| (mi +m2) 2 



gives the S s partial cross section of equation ()3.13p . On the basis of 
formula 1)3.211) the kinematically available phase space volume {Vps) can 
be as easily calculated as the excess energy Q. Close-to-threshold - in the 
range of few tens of MeV - the non-relativistic approximation differs only 
by a few per cent from the full solution given in equation (|3.20p . which 
in fact with an up-to-date computer can be solved numerically with little 
effort. Therefore, in the following chapters we will describe the data as 
a function of Vps as well as of Q or r^^, if it is found to be appropriate. 



3.4 Orientation of the emission plane 

Angles are the invention of Man; 
God uses vectors! [88] 

Colin Wilkin 

For the full description of the three particle system five independent 
variables are required. In the center-of-mass frame, due to the momentum 
conservation, the momentum vectors of the particles are lying in one 
plane often referred to as the emission-, reaction-, or decay plane. In this 
plane (depicted in figure 13. 3|) a relative movement of the particles can be 
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described by two variables only. As anticipated in the previous section, 
the square of the invariant masses of the di-proton and proton-?] system 
denoted as Spp and Sp^, respectively, constitute a natural choice for the 
study of the interaction within the pp?7 system. This is because in the 
case of non- interacting objects the surface spanned by these variables is 
homogeneously populated. The interaction among the particles modifies 
that occupation density and in consequence facilitates an easy qualitative 
interpretation of the experimental results. 











COS 6n = ^ 





3e: = i 



Figure 3.3: Definition of the ccntro-of-mass kinematical variables used for the description 
of the pp?7 system. In the centre-of-mass frame, the momenta of ejectiles lie in an emission 
plane. Within this plane the relative movement of the particles is fixed by the square of the 
invariant masses Spp and Sp,,. As the remaining three variables needed to define the system 
uniquely we use either 0*, 9*, and i/j shown in the upper panel or 0^, 9lr, and tpN defined 
in the lower panel. iV is a vector normal to the emission plane, which can be calculated as 
the vector product of the centre-of-mass momentum vectors of the outgoing protons. As an 
example two extreme orientations of the emission plane are shown in the right-lower panel. 
For further descriptions see the text. 
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The remaining three variables must define an absolute orientation of 
the emission plane in the distinguished coordinate system. This may be 
attained - for example - by defining the orientation for the momentum of 
the arbitrarily chosen particle in the center-of-mass frame and the angle 
which describes the rotation around the direction fixed by that particle. 
In one of our choices (following references [83,89]) the corresponding vari- 
ables are the polar and the azimuthal angle of the rj momentum vector, 
depicted in figure IX^ as 0* and 9*, respectively, and the angle ip describ- 
ing the rotation around the direction established by the momentum of 
the r] meson. Figure ESI demonstrates that such rotation neither affects 
the rj meson momentum nor changes the configuration of the momenta 
in the emission plane. For experiments with unpolarized beams and tar- 
gets the only favoured direction is the one of the beam. Therefore, as 
a zero value of the ip angle we have chosen the projection of the beam 
direction on the plane perpendicular to the momentum vector of the r) 
meson. Note that we identify the z-axis with the beam direction. The 
0*, 9*, and ip variables can also be interpreted as Euler angles allow- 
ing for the rotation of the emission plane into a xz-plane. The angle ip 
may be calculated as an angle between the emission plane and the plane 
containing momentum vectors of the r) meson and the beam proton, or 
correspondingly as an angle between the vectors normal to these planes. 
The angle i/j is equal to zero when these two normals are parallel. 

As a second possibility we will describe an orientation of the emis- 
sion plane by the azimuthal and polar angle of the vector normal to that 
plane [261]. These angles are shown in figure IT^ as 0^ and 9*j^, respec- 
tively. Further the absolute orientation of the particles momenta in the 
emission plane will be described hj ipM, the angle between the t] meson 
and the vector product of the beam momentum and the vector A^. 

Obviously, the interaction between particles does not depend on the 
orientation of the emission plane, and therefore, it will fully manifest it- 
self in the occupation density of the Dalitz plot which in our case will be 
represented in terms of the square of the invariant masses of the two par- 
ticle subsystems. Yet, the distribution of the orientation of the emission 
plane will refiect the correlation between the initial and final channels 
and hence its determination should be helpful for the investigation of the 
production mechanism. 



4. Low energy interaction within 

the ppr] and pprj' systems 

It frequently happens that, when particles are pro- 
duced in a nuclear or elementary particle reaction, 
some of these interact among themselves so strongly 
that they influence appreciably the properties of the 
reaction cross section. These interactions we 
shall call "final state interaction" [66]. 

Kenneth Watson 

In general, for a three-body exit channel one expects an energy depen- 
dence of the total cross section which can be described by the linear 
combination of partial cross sections from equation (j3.13|) . Therefore, 
to extract information about the final state interaction of the outgoing 
particles the contributions originating from different partial waves have 
to be known precisely. Appropriately, close-to-threshold there is only one 
important combination of angular momenta of emitted particles (Ss) and 
in this region the energy dependence of the total cross section is uniquely 
determined and hence the interpretation of results is significantly simpli- 
fied. This is one of the most important advantages of the meson creation 
at its production threshold, however, in order to exploit this fact prop- 
erly it is essential to determine the range of its applicability, which - as 
will be demonstrated in the next section - changes significantly with the 
mass of the produced meson. 

4.1 Range of the dominance of the ^Pq -^^Sqs 
transition 

... the most exact of sciences are those which are 
most concerned with primary considerations: for sci- 
ences based on few assumptions are more exact than 
those which employ additional assumptions ... [90]. 

Aristotle 

Investigations with polarized beams and targets [78,91] of the pp — ^ ppji^ 
reaction allowed to deduce that the Ss partial-wave accounts for more 
than 95 % of the total cross section up to rjM ~ 0.4, as can be seen in 
figure HH^, where the Ss contribution is indicated by the dotted line. The 
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Ss contribution was inferred assuming the ri% and rj^^ dependence for Ps 
and Pp partial waves, respectively. These are power-laws taken from 
proportionality ()3.13|) . which was derived under the assumption of non- 
interacting particles. Relatively small values of ^Po-wave nucleon-nucleon 
phase-shifts at low energies (compared to ^Sq phase-shifts in figure l^?Tb ) 
and similarly weak low-energy interactions of P-wave protons in other 
spin combinations [92] justify this assumption. 




Figure 4.1: (a) Decomposition of the total cross section of the pp ppn'^ reaction into 
Ss, Ps, and Pp final state angular momenta. The dashed and solid lines represent the rjf j 
and rj^j dependence of Ps and Pp partial cross sections, respectively. The remainder is 
indicated as the dotted line. The Sp partial wave is forbidden by the conservation laws and 
the Pauli excluding principle. Note that at rjM ~ 1 the Pp and Ps partial waves seem to 
dominate. However, the analysis of the differential cross section measured at CELSIUS at 
Vm ~ 0.449 [93, 117] showed that also a d-wave pion production - due to the interference 
between Ss and Sd states - constitutes 7 % of the total cross section, when a meson-exchange 
model is assumed. The figure has been adapted from [78]. (b) The ^So and ^Pq phase-shifts 
of the nucleon-nucleon potential shown versus the centre-of-mass kinetic energy available in 
the proton-proton system. The values have been extracted from the SAID data base [94] 
(solution SM97). For higher energies the S- and P-wave phase-shifts are nearly the same. 
This is because the collision parameter required to yield the angular momentum of 1 ft 
diminishes significantly below 1 fm with increasing energy and consequently the interaction 
of nucleons - objects of about 1 fm size - becomes almost central. 

In accordance with the phenomenology of Gell-Mann and Watson [77] 
described in section lT^ one expects that also in the case of heavier mesons 
the Ss partial wave combination will constitute the overwhelming fraction 
of the total production cross section for tjm smaller than 0.4. This implies 
- as can be deduced from the relation between tjm and Q illustrated in 
figure 14.2k - that mesons heavier than the pion are produced exclusively 
via the Ss state in a much larger excess energy range and hence larger 
phase space volume (see figure 1^^ ). Thus, whereas for 7r° production 
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the onset of higher partial waves is observed at Q around 10 MeV it is 
expected only above 100 MeV and above ~ 40 MeV for rj' and rj mesons, 
respectively. 




10 

Q [ MeV : 



a 




b) 



Figure 4.2: (a) The variable tjm as a function of the excess energy for tt", rj, and rj' 
mesons produced via pp — > pp Meson reactions, (b) The phase space volume Vps (defined by 
equation 13.201 ') versus the excess energy Q. The picture indicates that for the production of 
"heavy mesons" in the nucleon-nucleon interaction at a given Q value there is only a slight 
difference of the available phase space volume on the produced meson mass, which is larger 
than that of tt" production by about 30 % only. Therefore, for the comparative studies of 
the production dynamics of different mesons, Q is as much a suitable variable as Vps- Note, 
that the dashed and dotted lines are almost undistinguishable. 
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Figure 4.3: Differential cross section of the pp — > pprj reaction as a function of the meson 
centre-of-mass polar angle. Dashed lines indicate the isotropic distribution. Shown are 
results of measurements for Q = 15.5 MeV [1,18] (a) and Q = 41 MeV [95] (b). Only 
statistical errors are plotted, which in figure (a) are smaller than the size of symbols. The 
distribution presented in picture (a) is consistent with a measurement performed at an excess 
energy of Q = 16 MeV by means of the PRO MICE/ WAS A detector [96], whereas the data 
at Q = 37 MeV also from reference [96] deviate significantly from isotropy. However, data 
shown in picture (b) have been taken with a detector of much higher angular acceptance. 
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Figures 14.3b and 14.3b present the angular distributions of the cre- 
ated 7] meson in proton colhsions. It is evident that at Q = 15.5 MeV 
and still at Q = 41 MeV the production of the rj meson is completely 
isotropic within the shown statistical errors. Although at Q = 41 MeV 
the accuracy of the data does not exclude a few per cent of contribu- 
tions originating from higher partial waves, the dominance of the s- 
wave creation is evident. Similarly, the measurements of the differen- 
tial cross section (figure 14 .41) for the pp pprj' reaction performed at 
SATURNE [97] at Q = 143.8 MeV and at COSY [5] at Q = 46.6 MeV 
are still consistent with pure Ss-wave production, though the relatively 
large error bars would allow for other contributions on a few per cent 
level (~ 10%). The observation, reported in section that distri- 
butions of the orientation of the emission plane (fig. 16. 23^ and angle 
ip (fig. I6.2(jp are anisotropic, also suggests a small contribution of higher 
order partial waves. Specifically, the modulus of the cosine of the po- 
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Figure 4.4: Differential cross section of the pp — > ppr]' reaction as a function of the meson 
centre-of-mass polar angle. Dashed lines indicate the isotropic distribution. Shown are 
results of measurements for the pp — > ppr]' reaction taken at Q = 143.8 MeV [97] (a) and at 
Q = 46.6 MeV [5] (b). 

lar angle of the vector normal to the emission plane, whose distribution 
is found to be anisotropic, is a function of the sixth power of the final 
momenta {\cos{e*,^)\ ~ ^J\{pprotor^^_ x Pprotona) • pieamP ), aud even more pecu- 
liar the sine of the angle which is a function of the eighth power of 
the final momenta (|sm(^/))| ~ VKft^-otoni x Pprotor,^) ■ {p-n x VbeamW is found 
to be anisotropic too. This could be a sign of an interference between 
partial waves higher than Ss-wave ^. 



^ We are grateful to Colin Wilkin for sharing with us this astonishment. 
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4.2 Influence of the pp?7 and ppr]' interaction on 
the excitation function of the pp pprj{rj') 
reaction 

The only proof capable of being given that an object 
is visible, is that people actually see it. The only 
proof that a sound is audible, is that people hear it: 
and so of the other sources of our experience [98]. 

John Stuart Mill 

Let us now consider to what extent the energy dependence of the 
total cross section in the estimated range of the dominance of the Ss 
partial waves can be understood in terms of the phase space variation 
and the interaction between the particles participating in the reaction. 
Watson [66] and Migdal [99] proposed the factorization of the amplitude 
when the production is of short- and the interaction among the outgoing 
particles of long range. This requirement is well fulfilled for the close-to- 
threshold meson production due to the large momentum transfer (Ap) 
between the interacting nucleons needed to create the considered mesons 

Table 4.1: Momentum transfer Ap calculated according to equation 13.141 and the corre- 
sponding distance R ~ ft/Ap probed by the A'^A'' N N Boson reaction at the kinematical 
threshold for different particles produced. The table has been adapted from [100]. 



particle 


mass [MeV] 


Ap [fm-i 


R [fm] 


7 





0.0 


oo 


TT 


140 


1.9 


0.53 




550 


3.9 


0.26 




780 


4.8 


0.21 


Tj' 


960 


5.4 


0.19 




1020 


5.6 


0.18 



According to the Heisenberg uncertainty relation the large momen- 
tum transfer brings about a small space in which the primary creation 
of the meson takes place. In table 14.11 the distance probed by the 
A^A^ — >■ NN Meson reaction at threshold is listed for particular mesons. 
It ranges from 0.53 fm for pion production to 0.18fm for the meson, 
whereas the typical range of the strong nucleon-nucleon interaction at 
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low energies determined by the pion exchange may exceed a distance of 
a few Fermi and hence is by one order of magnitude larger than the values 
listed in table lTTl Thus in analogy to the Watson-Migdal approximation 
for two-body processes [66] the complete transition matrix element of 
equation (j3.ip may be factorized approximately as ^ 

\Mpp^ppx\^ ^ \Mfsi\' ■ |Mop ■ Fisi, (4.1) 

where Mq represents the total short range production amplitude, Mpsi 
describes the elastic interaction among particles in the exit channel and 
Fisi denotes the reduction factor accounting for the interaction of the 
colliding protons. Further, in the first order approximation one assumes 
that the particles are produced on their mass shell and that the created 
meson does not interact with nucleons. This assumption implies that the 
\Mpsi\'^ term can be substituted by the square of the on-shell amplitude 
of the nucleon-nucleon elastic scattering: 

\Mfsi\^ = \Mnn^nn\^. (4.2) 

Effects of this rather bold assumption will be considered later, when 
comparing the estimation with the experimental data. 

In the frame of the optical potential model the scattering amplitude 
is determined by phase-shifts. Particularly, the ^Sq proton-proton partial 
wave - relevant for further considerations - can be expressed explicitly 
as follows [102]: 

M = - sm Opp^ 5oj 

where C denotes the square root of the Coulomb penetration factor. 

determines the ratio of the probability of finding two particles close 
together to the probability of finding two uncharged particles together, 
all other things being equal [103] and can be expressed as [104]: 

= (4.4) 

where rjc is the relativistic Coulomb parameter, which for the collision of 
particles i,j reads: 

q, q . a s^ - mf - mj 

Vc = = qiq,a- 



A(sij,m2,m2) 



^ For a comprehensive discussion of tlie FSI and ISI issue including a historical overview and a 
criticism of various approaches the reader is referred to [101]. 
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with the fine structure constant a, the relative velocity v of the colliding 
particles and with q^, q^- denoting their charges ^. 

The variable k in equation 14.31 stands for either proton momentum in 
the proton-proton rest frame and the phase-shift is indicated by Spp. 
The phase-shifts Spp(}So) can be extracted from the SAID data base 
(see fig. 14.1b ) or, alternatively, can be calculated according to the modi- 
fied Cini-Fubini-Stanghellini formula including the Wong-Noyes Coulomb 
correction [107-109] 

k ctg{6pp) + 2kvM = _± + ^ _ (4.5) 

dpp ^ i + y^pp K 

where /i(r/,) = -ln{r],) - 0.57721 + r^^ ^ [103]. 

The phenomenological quantities app = — 7.83fm and bpp = 2.8 fm 
denote the scattering length and effective range [107], respectively. The 
parameters Ppp = 0.73 fm^ and Qpp = 3.35 fm^ are related to the detailed 
shape of the nuclear potential and derived from a one-pion-exchange 
model [107]. Substituting equation (|4.5p into equation (|4.3p allows to 
calculate the low-energy amplitude for the proton-proton elastic scatter- 

4 

mg : 

\M P = - 

(4.6) 

The result is presented as a solid line in figure 14.5b . and is in good agree- 
ment with the values obtained from the phase-shifts of the VPI partial 
wave analysis [94], shown as solid circles and with the phase-shifts of 
the Nijmegen analysis [110], shown as open squares. The factor is 
always less than unity due to the Coulomb repulsion between protons. 
At higher energies, where is close to unity, the nuclear scattering will 
be predominant and for the very low energies the Coulomb and nuclear 
interactions are competing. The Coulomb scattering dominates approx- 
imately up to about 0.8 Me V of the proton energy in the rest frame of 
the other proton, where equals to one-half [103]. In the case of the 
pp pp Meson reaction the maximum possible energy of a proton seen 
from another proton is equal to 0.8 MeV already at an excess energy of 
about Q = 0.4 MeV. Therefore, a significant influence of the Coulomb 



^ For collisions at an angular momentum I larger than h the of equation 14.41 needs to be 
multiplied by a factor of Yln=i (l + iVc/n)^) [105,106]. 

^ In principle the formula is valid for k< 133 MeV/c [108]. \Mpp^pp\'^ is taken to be constant 
for larger values of k. 
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repulsion on the energy dependence of the total production cross sec- 
tion is expected only at very low excess energies, i.e. conservatively for 
Q < 2 MeV. 




Figure 4.5: (a) Square of the proton-proton scattering amphtude versus k, the proton 
momentum in the proton-proton subsystem from [102,107] (solid line), [111] (dashed line), 
and [112,113] (dotted line). The filled circles are extracted from [94] and the opened squares 
from [110]. The curves and symbols have been arbitrarily normalized to be equal at maximum 
to the result from reference [111], shown as the dashed line, (b) Total cross section for the 
pp ppir" reaction as a function of the centre-of-mass excess energy Q. Data are from 
refs. [114-117]. The dashed line indicates a phase space integral normalized arbitrarily. The 
phase space distribution with inclusion of proton-proton strong and Coulomb interactions 
fitted to the data at low excess energies is shown as the solid line. The dotted line indicates 
the parametrization of reference [118] written explicitly in equation 14.7|l . with e — 0.3. 

Assuming that the on-shell proton-proton amplitude determines ex- 
clusively the phase space population one can obtain the total cross sec- 
tion energy dependence substituting equation ()4.6p into formula ()3.19|) . 
The solid line in figure l^3b represents the determined dependence for the 
pp — i> ppiT^ reaction. The absolute scale of the calculations has been fixed 
by normalizing to the data. One recognizes the good agreement with the 
experimental points in the excess energy range up to Q ~ 10 IvIeV in 
agreement with the previous conclusions based on the polarisation ob- 
servables. 

The dotted line in figures 14.5b which is practically indistinguishable 
from the solid line, presents the excess energy dependence of the total 
cross section taking into account the proton-proton FSI effects according 
to the model developed by Faldt and Wilkin [118,119]. Representing 
the scattering wave function in terms of a bound state wave function 
the authors derived a closed formula which describes the effects of the 
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nucleon-nucleon FSI as a function of the excess energy Q only. This ap- 
proach is specifically useful for the description of the spin-triplet proton- 
neutron FSI, due to the existence of a bound state (deuteron) with the 
same quantum numbers. Though a bound state of the proton-proton 
system does not exist, the model allows to express the total cross section 
energy dependence for a pp ^ pp Meson reaction by a simple and easily 
utilizable formula: 




(4.7) 



where the parameter e has to be settled from the data. The flux factor F 
and the phase space volume Vps are given by equations ()3.5p and ()3.2H) . 
respectively. The normalization can be determined from the fit of the 
data which must be performed for each reaction separately. 

The determined energy dependences of the total cross section for 
Tj' [3,5,11,53,97,120] and r] [10,50-53,121] mesons production in proton- 
proton collisions are presented in figure 14.61 Comparing the data to the 
arbitrarily normalized phase space integrals (dashed lines) reveals that 
the proton-proton FSI enhances the total cross section by more than an 
order of magnitude for low excess energies. 

One recognizes also that in the case of the rj' the data are described 
very well (solid line) assuming that the on-shell proton-proton amplitude 
exclusively determines the phase space population. This indicates that 
the proton-r]' interaction is too small to manifest itself in the excitation 
function within the presently achievable accuracy. In the case of i] meson 
production the interaction between nucleons is evidently not sufficient to 
describe the increase of the total cross section for very low and very high 
excess energies, as can be concluded from the comparison of the data and 
the upper solid line in figure This line was normalized to the data at 
an excess energy range between 15 MeV and 40 MeV. The enhancement 
of the total cross section for higher energies can be assigned to the outset 
of higher partial waves. As expected from the previous considerations, 
this is indeed seen at Q ~ 40 MeV where the energy dependence of the 
total cross section starts to change its shape. On the contrary, the close- 
to-threshold enhancement - being by about a factor of two larger than 
in the case of the 7r° and rj' mesons - can be assigned neither to the 
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contribution from other than Ss partial waves nor to the variation of the 
primary production amphtude \Mq\. The latter is expected to change 
at the most by a few per cent for excess energies below 20MeV [124]. 
Instead, this discrepancy can be plausibly explained by the influence of 
the attractive interaction between the f] meson and the proton. 




Q [ MeV] 

Figure 4.6: Total cross section for the reactions pp —> ppr/ (circles) and pp —> pprj (squares) 
as a function of the centre-of-mass excess energy Q. Data are from refs. [3, 5, 10, 11, 50-53, 
97,120,121]. The dashed lines indicate a phase space integral normalized arbitrarily. 
The solid lines show the phase space distribution with inclusion of the ^So proton-proton 
strong and Coulomb interactions. In case of the pp — > ppr) reaction the solid line was 
fitted to the data in the excess energy range between 15 and 40 MeV. Additional inclu- 
sion of the proton-?; interaction is indicated by the dotted line. The scattering length of 
o-pri = 0.7fm + i0.4fm and the effective range parameter hpr, = — 1.50 fm — i 0.24 fm [36] have 
been chosen arbitrarily. The dashed-dotted line represents the energy dependence taking into 
account the contribution from the ^Po -^^Sos, ^So ^^Pqs and ^'^P2S transitions [79]. 
Preliminary results for the ^Po —^''Sos transition with full treatment of three-body effects are 
shown as a dashod-double-dotted line [122, 123]. The absolute scale of dashed-double-dotted 
line was arbitrary fitted to demonstrate the energy dependence only. 
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Note, that the real part of the scattering length of the r^-proton poten- 
tial - depending on the analysis method and the studied reaction - is 3 
to 10 times [38, 125] larger than the scattering length for vr^-proton scat- 
tering (cp^ ~ 0.13 fm) [126]. Hence, the modifications of the total cross 
section energy dependence due to vr^'-proton and r^'-proton interaction are 
too weak to be observed within the up-to-date accuracy of measurements 
and calculations. However, the influence of the ?7-proton interaction pre- 
sented in figure l^ni is evident and hereafter it will be considered whether 
it may serve for the estimation of the 77-proton scattering parameters. 

A strict quantitative calculation requires the evaluation of the three- 
body Faddeev equation which is out of the scope of the present work^. 
Here, we will rather present a simple phenomenological treatment which 
shall lead to the qualitative understanding how the mutual interaction 
among three outgoing particles affects the excitation function. One of 
the simplest possibilities based on the naive probabilistic interpretation 
of the incoherent pairwise interaction would be to factorize the overall 
enhancement factor into corresponding pair interactions [68,130]: 
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where 



P denotes the square of the elastic scattering amplitude of 



particles i and j. The \Mpp^pp\ term can be evaluated according to the 
formula ()4.6|) . which for the s-wave r^-proton scattering, after substitution 
of = 1 and rjc = Ppp = 0, reduces to 
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where 



A(s^ 



m^, m; 



denotes the 77 momentum in the proton-?] rest frame. The scattering 
length ttpri and effective range are complex variables with the imagi- 
nary part responsible e.g. for the pr] pir^ conversion. 



^ An exact derivation of the Faddeev equation can be found for example in [127]. Presently few 
theory groups are carrying on the corresponding calculations [82,122,123,128,129], 

^ Note that the sign of the term —1/a from equation 14. 6i was changed because the imaginary part 
of the proton-?? scattering length is positive [36], as we also adopted here, whereas in the majority of 
works concerning nucleon-nucleon interaction, the scattering length is negative [64], We are grateful 
for this remark to A. Gasparyan. 
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The factorization of both i) the overall production matrix element 
(eq. (|4.1|) ) and ii) the three particle final state interactions (eq. (|4.8p ) 
applied to formula (j3.19|) gives the following expression for the total cross 
section of the pp pprj reaction: 



F 4s ' '^S"" 



pp \-ir±pp^ppy 



(4.10) 

where the p roto ns are distinguished by subscripts. Exploring formu- 
lae ()4.fjj) and ()4.9|1 gives the results shown as the dotted line in figure WM 

Evidently, the inclusion of the proton-?] interaction enhances the total 
cross section close-to-threshold by about a factor of 1.5 and leads to 
a better description of the data. A similar effect close-to-threshold is 
also observed in the data of photoproduction of r] via the ■yd — ^ pnrj 
reaction [131] indicating to some extent that the phenomenon is inde- 
pendent of the production process but rather related to the interaction 
among the rj meson and nucleons in the Sii(1535) resonance region. How- 
ever, although a simple phenomenological treatment [4, 68, 130] - based 
on the factorization of the transition amplitude into the constant pri- 
mary production and the on-shell incoherent pairwise interaction among 
outgoing particles - describes the enhancement close-to-threshold (dot- 
ted line) very well, it fails to describe the invariant mass distribution 
of the proton-proton and proton-r/ subsystems determined recently at 
Q = 15 MeV by the COSY-TOF [95] and at Q = 15.5 MeV by the 
COSY-11 [1] collaborations. Though the groups utilized entirely differ- 
ent experimental methods the obtained results agree very well with each 
other. The structure of these invariant mass distributions may indicate 
a non-negligible contribution from the P-waves in the outgoing proton- 
proton subsystem [79], which can be produced for instance via ^5*0 -^^PqS 
or -^^P2S transitions. This hypothesis encounters, however, difficul- 
ties in describing the excess energy dependence of the total cross section. 
The amount of the P-wave admixture derived from the proton-proton 
invariant mass distribution leads to a good description of the excitation 
function at higher excess energies while at the same time it spoils signif- 
icantly the agreement with the data at low values of Q, as depicted by 
the dashed-dotted line in figure 14.61 However, these difficulties in repro- 
ducing the observed energy dependence might be due to the particular 
model used in reference [79], and thus higher partial wave contributions 
cannot be excluded a priori. 

In contrast to the P-wave contribution the three-body treatment [122, 
123] of the ppr] system (dashed-double-dotted line) leads to even larger 
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enhancement of the cross section near threshold than that based on the 
Ansatz of the factorization of the proton-proton and proton-77 interac- 
tions. It must be kept in mind, however, that too strong FSI effect 
predicted by the three-body model must be partially assigned to the ne- 
glect of the Coulomb repulsion in this preliminary calculations [122,123]. 
The above considerations illustrate that the simple phenomenological ap- 
proach shown by the dotted line could fortuitously lead to the proper re- 
sult, due to a mutual cancellation of the effects caused by the approxima- 
tions assumed in the calculations and the neglect of higher partial waves. 
They show also unambiguously that for the complete understanding of 
the low energy ppr] dynamics, in addition to the already discussed exci- 
tation function of the total cross section, a knowledge of the differential 
observables is necessary. These will help to disentangle effects caused by 
the proton-T] interaction and the contributions from higher partial waves. 
This issue will be discussed further in section 031 where we present dis- 
tributions determined experimentally for two sets of orthogonal variables 
fully describing the ppr] system, which was produced at an excess energy 
of Q = 15.5 MeV via the pp —* ppr] reaction using the COSY-11 [8,61] 
facility at COSY [62]. 

4.3 Phenomenology of the proton-proton initial 
and final state interaction 

Thus, it is suggested that among created beings 
there must be some basic agent which will move 
things and bring them together [90]. 

Aristotle 

When reducing the proton-proton FSI effect to a multiplicative factor, 
one finds that it depends on the assumed nucleon-nucleon potential and 
on the produced meson mass [132]. 

This issue was recently vigorously investigated e.g. by authors of ref- 
erences [71,100,101,132,133] and we shall briefly report this here as well. 
Up to now we factorized the transition matrix element into a primary 
production of particles and its on-shell rescattering in the exit channel 
(eqs. (j4.1|)(j4.2j)(j4.8|) ). Though it is a crude approximation, neglecting 
the off-shell effects of the production process completely, it astoundingly 
leads to a good description of the energy dependence of the total cross 
section, as already demonstrated in figures 14.5b and 14.61 The off-shell 
effects, as pointed out by Kleefeld [101], could have been safely neglected 
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in case of the electromagnetic transitions in atoms or (3 decays, where the 
excitation energy of the involved nucleons is by many orders of magni- 
tude smaller than their masses and the initial and final states go hardly 
off-shell [101]. However, in the case of the A^A^ NN Meson process 
the large excitation energy of the colliding nucleons is comparable with 
the nucleon masses and the primary interaction may create the particles 
significantly far from their physical masses, so that a priori the off-shell 
effects cannot be disregarded. 



A 



a) 




b) 



Figure 4.7: Diagrammatic representation of the DWBA expressed by equation 14.1111 . (a) 
The primary production term, (b) The loop diagram including the nucleon-nucleon FSI. 
r(k',k) stands for the half-off-shell T matrix with k and k' denoting the on- and off-shell 
centre-of-mass momentum in the nucleon-nucleon system, respectively, q indicates the mo- 
mentum of the created meson in the reaction centre-of-mass system and p the momentum 
of the colliding nucleons. 

Generally, the decomposition of the total production amplitude into 
the primary production and the subsequent nucleon-nucleon interaction 
visualized in figure 14.71 is expressed by the formula: 

M = M°" + M;;"GTnn, (4.11) 

where the second term of the equation represents the integration over the 
intermediate (k') momenta of the off-shell production amplitude and the 
half-off-shell nucleon-nucleon T matrix [132] . Assuming that the primary 
production occurs in such a way that one nucleon emits the meson which 
then re-scatters on the other nucleon and appears as a real particle the 
authors of reference [132] found that the enhancement of the cross section 
due to the nucleon-nucleon interaction depends strongly on the mass of 
the created meson. This is because with the increasing mass of the 
produced meson the distance probed by the nucleon-nucleon interaction 
decreases (see table 14. 1|) and hence the relevant range of the off-shell 
momenta becomes larger. The effect for the pp pp Meson reactions is 
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presented in figure where one can see that, when utihzing the Bonn 
potential model for the nucleon-nucleon T matrix, the enhancement in 
case of the 7r° production is by about a factor of four larger than for the 
1] or 7]' mesons. A similar conclusion, but with the absolute values larger 
by about 40 %, was drawn for the Paris A^A^ potential [132]. 




k[MeV/c] a) kfMeV/cl b) k [ MeV/c 1 



Figure 4.8: The FSI factor for Bonn [134] a) and Yamaguchi [135] b) potentials. The 
solid, dotted and dashed lines correspond to n", rj, and rj' meson production, respectively. 
Thick solid lines indicate the inverse of the squared lost function (|fpp(k)|^ = j J(— k)|~^). 
c) The FSI factors for tt" production of Paris (dotted curve) [136], Bonn (dashed curve), 
and Yamaguchi (solid curve) potentials normalized to be equal at maximum to the inverse 
of the squared Jost function of the Yamaguchi potential (thick solid line). The shapes 
stemming from different potentials are almost indistinguishable. Note that the thick solid 
line corresponds to the dashed line in figure 14.5b . The figures have been adapted from 
reference [132]. 

On the contrary, when applying the Yamaguchi potential into calcu- 
lations the enhancement grows with the increasing mass of the meson, 
as shown in figure 14.8b . The thick solid curves in figures 14. 8h . and 14.8b 
show the results of the frequently applied approximation of the nucleon- 
nucleon FSI effects: 

M = M°" + M°f^GTNN ~ M°"-(l + GTnn) = 

= M°"J-i(-k) = M^F^^ik), (4.12) 

where the overall transition matrix element M is factorized to the pri- 
mary on-shell production and the A^A^ FSI expressed as the inverse of the 
Jost function J~^(— k) [113]. As can be seen in figures I^TSk and 14.8b . the 
variation of the absolute values - of such obtained enhancement factors - 
with the applied potential is significant. 

Since the physical value of the total cross section cannot depend on the 
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off-shell features of the potential used in calculations, which are in prin- 
ciple not measurable [137], the differences in the magnitude of the \Fpp\ 
factor must reflect itself in a corresponding dependence of the primary 
production amplitude on the potential used. Therefore, it is of vital im- 
portance to realize that the values of the threshold amplitudes \Mq\ are 
significant only in the context of the potential they were extracted from. 
Figure I^Ht demonstrates, however, that the shapes of the enhancement 
factors, with the meson exchange mechanism assumed for the primary 
production, are pretty much the same, independently of the applied NN 
potential and correspond to the form of the Jost function inferred from 
the Yamaguchi potential. This indicates that the energy dependence of 
the NN FSI factors is predominantly determined by the on-shell A^A^ T 
matrix. 

In references [71,100,101] the formula for the transition matrix element 
with explicit dependence on the considered off-shell features for the initial 
and final state interaction is derived: 

Uk)e2'*('^) - l] ■ [1 + P/(p,k)] ] ( kp)e^'*(P) - ll ■ [1 + P.ipM ] 

i 1 1^+^ i 1 ' 

(4.13) 

where subscripts i and / indicate the initial and final state, respectively. 
The functions P(p, k) exhibit all the off-shell effects of the NN inter- 
action and the primary production current [100] and 6 and t] denote 
the phase-shift and inelasticity, correspondingly. At threshold, the in- 
elasticity in the exit channel is equal to unity (rjik) = 1) due to the 
small relative momentum of the outgoing nucleons. The last term of 
the formula expresses the influence of the initial state interaction on the 
production process. Due to the large relative momenta of the colliding 
protons needed to create a meson it is characterized by a weak energy 
dependence in the excess energy range of a few tens of MeV. For example 
in figure ITTb one can see that the phase-shift variation of the ^Pq partial 
wave (having predominant contribution to the threshold production of 
pseudoscalar mesons) in the vicinity of the threshold for mesons heavier 
than 7r° is indeed very weak. Taking additionally into account that the 
initial state off-shell function Pi(p, k) is small (as it is the case at least for 
meson exchange models [71]) one reduces the influence of the A^A^ initial 
state interaction to the reduction factor Ffsi which can be estimated 
from the phase-shifts and inelasticities only: 

Fjsi = ] I r/(p) e^^^^p) + if . (4.14) 
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At the threshold for vr meson production this is close to unity since at 
this energy the inelasticity is still nearly 1 and the ^Pq phase-shift is close 
to zero (see figure ITTb ). However, at the t] threshold, where the phase- 
shift approaches its minimum, the proton-proton ISI diminishes the total 
cross section already by a factor of 0.2 [71]. A similar result was obtained 
using a meson exchange model for t] production in the pp ppt] reaction 
and calculating the proton-proton distortion from the coupled-channel 
ttNN model [72]. The authors of reference [72] concluded that the initial 
proton-proton distortion reduces the total cross section by about a factor 
of ~ 0.26, which keeps constant at least in the studied range of 100 MeV 
in kinetic beam energy. Hence, the closed formula ()4.14j) disregarding the 
off-shell effects (Pi(p, k)) permits to estimate the cross section reduction 
due to the initial state distortion with an accuracy of about 25 %. The 
shape of the ^Pq phase-shift shown in figure l^!Tb indicates that the effect 
is at most pronounced close to the 1] production threshold, yet for the rj' 
meson the formula (j4.14p leads to a factor Fjsi = 0.33 [138]. The primary 
production amplitude as well as the off-shell effects of the nucleon-nucleon 
FSI (-P/(p, k)) are also weakly energy dependent [100], since they account 
for the short range creation mechanism, which shall be considered in 
chapter El 

4.4 Qualitative comparison between ppr], ppr]', and 
ppiT^ interactions 

But what we must aim at is not so much to ascertain 
resemblances and differences, as to discover similar- 
ities hidden under apparent discrepancies [74]- 

Henri Poincare 

The accordance of the experimental data with the simple factorization 
represented by solid lines in figures 14.5b and 14.61 fully confirms the sup- 
positions considered in the previous sections which imply that the energy 
dependence of the total cross section is in the first order determined by 
the on-shell scattering of the outgoing particles. However, since the dis- 
tortion caused by the nucleons is by some orders of magnitude larger than 
that resulting from the meson-nucleon interaction, even small fractional 
inaccuracies in the description of nucleon-nucleon effects may obscure 
the inference of the meson-nucleon interaction. The differences between 
the square of the on-shell proton-proton scattering amplitude and the 
Jost function prescription are presented in figure 14.5b . To minimize the 
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ambiguities that may result from these discrepancies at least for the 
quantitative estimation of the effects of the unknown meson-nucleon in- 
teraction one can compare the spectra from the production of the meson 
under investigation to the spectra determined for the production of a 
meson whose interaction with nucleons is well established. For instance, 
to estimate the strength of the ripp and rj'pp FSI one can compare the 
appropriate observables to those of the vr^pp system. 




-1 2 3 

10 1 10 10 10 




Figure 4.9: The ratios of 

\M^j/\Mo°\ (upper panel) and 
\/\Mo I (lower panel) extracted 
from the data by means of equa- 
tion 14. ion , assuming a pp-FSI 
enhancement factor as depicted 
by the dotted line in figure I4.5h , 
and neglecting the proton-meson 
interaction [2]. 



Upper and lower panels of figure 14.91 show the dependence of |Mo| 
on the phase space volume for rj and rj' production normalized to \Mq \. 
The values of \Mq\ were extracted from the experimental data by means of 
equation ()4.10|) disregarding the proton-meson interaction (|Mp^(^/)^p^(^/) | 
was set to 1). If the influence of the neglected interactions were the same 
in the cases of the T]{ri') and 7r° production the points would be con- 
sistent with the solid line. This is the case for the pp pprj' reaction 
visualizing the weakness of the proton-r^' interaction independently of 
the prescription used for the proton-proton FSI [2]. In case of the rj' 
meson its low-energy interaction with the nucleons was expected to be 
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very weak since there exists no baryonic resonance which would decay 
into Nrj' [43] . Figure 14.9b shows - independently of the model used for 
the correction of the proton-proton FSI - the strong effects of the rjpp 
FSI at low Vpg. 

4.5 Dalitz plot occupation for the ppr] system 

... a knowledge of the effects is what leads to an 
investigation and discovery of causes [139]. 

Galileo Galilei 

The strength of the interaction between particles depends on their rel- 
ative momenta or equivalently on the invariant masses of the two-particle 
subsystems. Therefore it should show up as a modification of the phase 
space abundance in the kinematical regions where the outgoing particles 
possess small relative velocities. Only two invariant masses of the three 
subsystems are independent (see eq. 13.18^ and therefore the entire ac- 
cessible information about the final state interaction of the three-particle 
system can be presented in the form of the Dalitz plot. Figure lT.lUb shows 
distribution of the events as determined experimentally for the pprj sys- 
tem at an excess energy of Q = 15.5 MeV. This distribution originating 
from kinematically complete measurements comprises the whole experi- 
mentally available information about the interactions of the pprj-sjstem. 
In this figure one easily recognizes the growth of the population den- 
sity at the region where the protons have small relative momenta which 
can be assigned to the strong attractive S-wave interaction between the 
two protons. This is qualitatively in agreement with the expectation 
presented in figure l4.1Ub . which shows the result of Monte-Carlo calcu- 
lations where the homogeneously populated phase space was weighted 
by the square of the on-shell ^5*0 proton-proton scattering amplitude. 
However, already in this two dimensional representation it is visible that 
the experimentally determined distribution remains rather homogeneous 
outside the region of the small proton-proton invariant masses, whereas 
the simulated abundance decreases gradually with growing Spp (as in- 
dicated by the arrow). Figure l^.lOb shows the simulated phase space 
density distribution disregarding the proton-proton interaction but ac- 
counting for the interaction between the r^-meson and the proton. At 
this excess energy, corresponding to the small relative momentum range 
(k^"^ ^ 105MeV/c) the variations of the proton-r] scattering amplitude 
are quite moderate (see figure l^?T3b in the next section). 
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Figure 4.10: Dalitz plots. Lines surrounding the plots show the kinematical limits. 

(a) Experimental result determined for the pp — > pprj reaction at Q = 15.5 MeV. Data were 
corrected for the detection acceptance and efficiency. 

(b) Monte-Carlo simulations for the pp — > pprj reaction at Q = 15.5 MeV: Phase-space den- 
sity distribution modified by the proton - proton final state interaction. 

(c) Simulated phase space density distribution modified by the proton-r; interaction with a 
scattering length of flpr, = 0.7 fm -|- i 0.3 fm. Details of the calculations together with the 
discussion of the nucleon-nucleon and nucleon-meson final state interaction can be found in 
reference [4]. The scale in the figure is linear in contrast to the above panels. 

(d) Dalitz plot distribution of the pp — > ppr] reaction at an excess energy of Q = 37.6 MeV, 
corrected for the detection acceptance. Out of the two invariant masses corresponding to 
two p — rj pairs the one being larger is plotted along the x-axis. The figure is taken from 
reference [50]. A similar spectrum for the ppyr" system [114] reveals the infiuence of the pro- 
ton-proton interaction, yet again, as for the total cross section energy dependence (fig. 14. 5b ). 
the proton- tt" interaction is too weak to affect the density distribution of the Dalitz plot ob- 
servably. 

Due to the low strength of this interaction the expected deviations 
from a uniform distribution is smaller - from that caused by the proton- 
proton force - by about two orders of magnitude, yet an enhancement of 
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the density in the range of low invariant masses of proton-?] subsystems 
is clearly visible. Note that the scale in the lower panel of figure 14.101 is 
linear whereas in the upper panel it is logarithmic. Due to weak varia- 
tions of the proton-r; scattering amplitude the enhancement originating 
from the ?7-meson interaction with one proton is not separated from the 
?7-meson interaction with the second proton. Therefore an overlapping of 
broad structures occurs. It is observed that the occupation density grows 
slowly with increasing Spp, opposite to the effects caused by the S-wave 
proton-proton interaction, but similar to the modifications expected for 
the NN P-wave [87]. From the above example it is obvious that only 
from experiments with high statistics, signals of the meson-nucleon in- 
teraction can be observed on top of the overwhelming nucleon-nucleon 
final state interaction. The increase of the distribution density at regions 
of small invariant masses of the proton-proton and proton-77 subsystems 
is also visible in figure l^!TUl i showing experimental results determined at 
Q = 37.6 MeV by the WASA/PROMICE collaboration. At this excess 
energy these regions are quite well separated. However, since this is close 
to the energy where the advent of higher partial waves is awaited, the 
possible contribution from the P-wave proton-proton interaction cannot 
be a priori excluded. Specifically, the latter leads to the enhancement 
at large invariant masses of the proton-proton pair [87] and hence af- 
fects the phase space region where the modification from the proton-r/ 
interaction is expected. A deviation of the experimentally observed pop- 
ulation of the phase space from the expectation based on the mentioned 
assumptions is even better visible in figure 14.111 This figure presents the 
projection of the phase space distribution onto the Spp axis corresponding 
to the axis indicated by the arrows in figures 14.10b and 14.10b . The su- 
perimposed lines in figure I^TTl correspond to the calculations performed 
under the assumption that the production amplitude can be factorized 
into a primary production and final state interaction. The dotted lines 
result from calculations where only the proton-proton FSI was taken into 
account, whereas the thick-solid lines represent results where the over- 
all enhancement was factorized into the corresponding pair interactions 
of the pprj system. This factorisation Ansatz is only valid if the dif- 
ferent amplitudes are completely decoupled which is not the case here. 
Therefore, these calculations should be considered as a rough estimate 
of the effect introduced by the FSI in the different two body systems. 
As introduced in section 03 the enhancement factor accounting for the 
proton-proton FSI has been calculated [2, 4] as the square of the on- 
shell proton-proton scattering amplitude derived according to the modi- 
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fied Cini-Fubini-Stanghellini formula including the Wong-Noyes Coulomb 
corrections [108]. The homogeneous phase space distributions (thin solid 
lines) deviate strongly from the experimentally determined spectra. The 
curves including the proton-proton and proton-?] FSI reflect the shape 
of the data for small invariant masses of the proton-proton system, yet 
they deviate significantly for large Spp and small Sp^ values. An explana- 
tion of this discrepancy could be a contribution of P-wave proton-proton 
interaction [79], or a possibly inadequate assumption that proton-r] and 
proton-proton interaction modify the phase space occupations only as 
incoherent weights [101]. 




Figure 4.11: Distributions of the square of the proton-proton (spp) and proton-r; (sp^) 
invariant masses determined experimentally for the pp — > pprj reaction at the excess energy 
of Q = 15.5 MeV by the COSY-11 collaboration (closed circles), at Q = 15 MeV by the 
TOF collaboration (open circles) [95], and at Q = 16 MeV by PRO MICE/ WAS A (open 
triangles) [96]. The TOF and PRO MICE/ WAS A data have been normalized to those of 
COSY-11, since these measurements did not evaluate the luminosities but rather normalized 
the results to reference [50] (see also comment [140]). The integrals of the phase space 
weighted by the square of the proton-proton on-shell scattering amplitude FSIpp (dotted 
lines) , and by the product of FSIpp and the square of the proton-?7 scattering amplitude (thick 
solid lines), have been normalized arbitrarily at small values of Spp. The thick solid line was 
obtained assuming a scattering length of ap^ = 0.7 fm + i 0.4 fm. The expectations under 
the assumption of the homogeneously populated phase space are shown as thin solid curves. 

A slightly better description is achieved when the proton-proton in- 
teraction is accounted for by the realistic nucleon-nucleon potential. Fig- 
ure 14.12b depicts the results obtained using two different models for the 
production process as well as for the NN interaction [79,141,142]. The 
calculations for the ^Pq -^^Sqs transition differ slightly, but the differ- 
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ences between the models are, by far, smaller than the observed signal. 
Therefore we can safely claim that the discussed effect is rather too large 
to be caused by the particular assumptions used for the production op- 
erator and NN potential. 




Figure 4.12: (a) Distribution of the square of the proton-proton (spp) invariant mass for 
the pp ppri reaction at an excess energy of Q = 15.5 MeV. Solid and dashed lines corre- 
spond to the calculations under the assumption of the '^Pq Sqs transition according to 
the models described in references [141,142] and [79], respectively. The dotted curve shows 
the result with the inclusion of the ^So -^'^Pos contribution as suggested in reference [79]. 
(b) The same data as above but with curves denoting preliminary three-body calcula- 
tions [122, 123] of the final pprj system as described in [143]. At present only the dominant 
transition ^Pq -^^Sqs is taken into account and the production mechanism is reduced to the 
excitation of the 5'ii(1535) resonance via the exchange of the tt and rj mesons. The solid line 
was determined with the rigorous three-body approach [122, 123] where the proton-proton 
sector is described in terms of the separable Paris potential (PEST3) [144], and for the 
77-nucleon scattering amplitude an isobar model analogous to the one of reference [145] is 
used with a^jv — 0.5 f m -I- i • 0.32 fm. The dashed line is obtained if only pairwise interac- 
tions {pp+prf) are allowed. The effect of proton- proton FSl at small Spp is overestimated due 
to neglect of Coulomb repulsion between the protons. The lines are normalized arbitrarily 
but their relative amplitude is fixed from the model. 

As we will see in subsection 16 . 3 . 31 in figure lU^^Hb . the experimental dis- 
tribution of the 7] polar angle in the center-of-mass frame is fully isotropic. 
This is the next evidence - besides the shape of the excitation function 
and the kinematical arguments discussed in reference [4] - that at this 
excess energy (Q = 15.5 MeV) the r] meson is produced in the center- 
of-mass frame predominantly with the angular momentum equal to zero. 
Similarly, the distribution determined for the polar angle of the relative 
proton-proton momentum with respect to the momentum of the r] me- 
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son as seen in the di-proton rest frame is also consistent with isotropy. 
Anyhow, even the isotropic distribution in this angle does not imply di- 
rectly that the relative angular momentum between protons is equal to 
zero, because of their internal spin equal to |. Therefore, the contribu- 
tion from the ^Po-wave produced via the ^5*0 -^^Pos transition cannot 
be excluded. Moreover, as pointed out in reference [79], the isotropic 
angular distribution, can principally also be achieved by the destructive 
interference between the transitions ^Sq -^^PqS and ^^P2S. 

Figure 4.13: Distribution of square of 
the proton-proton invariant mass from 
the pp — > ppri reaction measured at 
COSY- 11 for the excess energy range 
4 MeV < Q < 5 MeV [1, 10, 146]. The 
superimposed line shows the result of simula- 
tions performed under the assumption that 
the phase space population is determined 
exclusively by the on-shell interaction between 
outgoing protons. The "tail" at large Spp 
values is due to the smearing of the excess 
energy, since this former COSY-11 data have 
not been kinematically fitted. Additionally to 
S [ GeV^/c^ ] ^ MeV range of Q a smearing of about 

0.3 MeV (a) should be taken into account. 

Since flat angular distributions do not preclude the occurrence of 
higher partial-waves an effect of their contribution to the production 
process was proposed to explain the structure observed in the invariant 
mass spectra [79]. In fact, as depicted by the dotted line in figure ll?T^ . 
an admixture of the ^Sq -^^PqS transition to the main ^Pq -^^SqS one 
leads to the excellent agreement with the experimental points. How- 
ever, at the same time, this conjecture leads to strong discrepancies in 
the shape of the excitation function as can be deduced from the com- 
parison of the dashed-dotted line and the data in figure 14.61 While it 
describes the data points in the excess energy range between 40 MeV 
and 100 MeV it underestimates the total cross section below 20 MeV 
by a factor of 2. This observation entails that either the contribution 
of higher partial waves, though conforming with invariant mass spectra 
and angular distributions, must be excluded or that the simplifications 
performed e.g. for the description of the final state three-body system 
were too strong. Interestingly, the enhancement at large Spp is visible 
also at much lower excess energy. This can be concluded from figure 14.131 
in which the COSY-11 data at Q ~ 4.5 MeV [1,146] are compared to the 
simulations based on the assumption that the phase space abundance is 
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due to the proton-proton FSI only. This observation could imply that the 
effect is caused by the proton-r] interaction rather than by higher partial 
waves, since their contribution at such small energies is quite improba- 
ble [4] . However, as shown in the lower part of figure 14.121 the rigorous 
three-body treatment of the ppi] system leads at large values of Spp to 
the reduction of the cross section in comparison to the calculation taking 
into account only first order rescattering {pp+pr]) [122,123]. Here, both 
calculations include only the ^Pq ^^SqS transition. Although the pre- 
sented curves are still preliminary, we can qualitatively assess that the 
rigorous three-body approach, in comparison to the present estimations, 
will on one hand enhance the total cross section near threshold as shown 
in figure 1321 while on the other hand it will decrease the differential cross 
section at large values of Spp. This is just opposite to the influence of 
P-waves in the proton-proton system. 

From the above presented considerations it is rather obvious that the 
rigorous three-body treatment of the produced ppr] system and the exact 
determination of the contributions from the higher partial waves may 
result in a simultaneous explanation of both observations: the near- 
threshold enhancement of the excitation function of the total cross sec- 
tion, and the strong increase of the invariant mass distribution at large 
values of Spp. For the unambiguous determination of the contributions 
from different partial waves spin dependent observables are required [79] . 
The first attempt has been already reported in [7] , and more comprehen- 
sive investigations are in progress [27] (see also section Ei3|) • 

4.6 Quasi-bound state 

... which exists either because of some impossible 
blip on the curve of probability or because the gods 
enjoy a joke as much as anyone [147]. 

Terry Pratchett 

... / am inclined to think that scientific discovery 
is impossible without faith in ideas which are of a 
purely speculative kind, and sometimes even quite 
hazy ... [75]. 

Karl Raimund Popper 

With the up-to-date experimental accuracy, from all meson- iV-iV systems 
the T]NN one reveals by far the most interesting features. Recently, the 
dynamics of the rjNN system has become a subject of theoretical in- 
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vestigations in view of the possible existence of quasi-bound or resonant 
states [148]. A direct measure of the formation - or non-formation - 
of an ?7-nuclear quasi-bound state is the real part of the r^-nucleon scat- 
tering length [150]. The determined values of Re (a^Tv) range between 
0.20 fm [37] and 1.05 fm [39] depending on the analysis method and the 
reaction studied [38], and at present an univocal answer whether the at- 
tractive interaction between the t] meson and nucleons is strong enough 
to form a quasi-bound state is not possible. 

Figure 4.14: Averaged squared pro- 
duction amplitudes of the reactions 
pd — > ^Herj and dd — » ^Her] as function 
of the centre-of-mass 77 momentum. The 
data are taken from references [151-153]. 
The amplitude /„ is related to the un- 
polarized centre-of-mass cross section by 
the equation ^ = ^l/nP, where p* 
and denote the centre-of-mass momen- 
tum in rj'^He and pd systems, respectively. 
Superimposed lines correspond to the fit 
of the s-wave scattering-length formula: 
f(p \^ coustans ^j^g^.^ deuotes the 

scattering length of the rj — ^He or ri^^He 
potential. 

The shape of the energy dependence of the pd ^Het] production 
amplitude (see figure implies that either the real or imaginary part 
of the 7] ^He scattering length has to be very large [176], which may be 
associated with a bound rj ^He system. Similarly encouraging are results 
of reference [154], where it is argued that a three-body rjNN resonant 
state, which may be formed close to the rjd threshold, may evolve into 
a quasi-bound state for Re (o^at) > 0.733 fm. Also the close-to-threshold 
enhancement of the total cross section of the pp pprj reaction was 
interpreted as being either a Borromean (quasi-bound) or a resonant 
ripp state [155] provided that Re (a,,Ar) > 0.7 fm. Contrary, recent cal- 
culations performed within a three-body formalism [148] indicate that 
a formation of a three-body riNN resonance state is not possible, inde- 
pendently of the rjN scattering parameters. Moreover, the authors of 
reference [156] exclude the possibility of the existence of an riNN quasi- 
bound state. Results of both calculations [148,156], although performed 
within a three-body formalism, were based on the assumption of a sep- 
arability of the two-body rjN and NN interactions. However, in the 
three-body system characterized by the pairwise attractive interactions, 
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the particles can be pulled together so that their two-body potentials 
overlap, which may cause the appearance of qualitatively new features 
in the rjNN system [148]. Particularly interesting is the rjd final state 
where the pair of nucleons alone is bound by the strong interaction. 




Figure 4.15: (a) Total cross section of the quasi-free pn — > d-q reaction as a function 
of the excess energy [121, 157]. The curve - fitted in amplitude - indicates the energy 
dependence proportional to the function -^Q • (1 + Q/83.5) which accounts for the s- and 
p-wave contribution. The ratio of the s- and p-waves magnitudes was taken the same as 
determined for the np — > dyr" reaction [158-160]. It was assumed that the p-wave to s-wave 
ratio is the same for the pn d-q and pn div^ reactions at a corresponding value of 
riM- (b) Arbitrarily normalized enhancement factors for pd-, d-q-, and pr;-FSI. The r;-proton 
factor is calculated according to equation l|4.9|l . with Opr, — 0.717 fm -I- i 0.263 fm [161] and 
bpr] ~ — 1.50fm — z0.24fm [36]. The enhancement factor for 77-deuteron has been extracted 
from the data of panel (a) parametrizing the ratio of the cross section to the phase space 
volume by the expression [162]: F^^ = 1 + 0.5/(0.5 + (Q/5)^). The proton-deuteron FSI 
factor is calculated according to ref. [163]. 

Figure 14.15b shows the total cross section for the pn — drj reaction 
measured close to the production threshold. For excess energies below 
10 MeV the data are enhanced over the energy dependence determined for 
the pn — *• dTT^ reaction indicated by the solid curve. This is in qualitative 
agreement with the calculations of Ueda [164] for the three-body riNN- 
ttNN coupled system, which predict the existence of an rjNN quasi- 
bound state with a width of 20 MeV. Ueda pointed out that the binding 
of the r]NN system is due to the Sn rjN and ^Si NN{d) interaction 
which is characterized by no centrifugal repulsion. Such repulsion makes 
the irNN system, in spite of the strong P33 vrA^ attraction, hard to be 
bound [165]. Whether the observed cusp at the pn dt] threshold is 
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large enough to confirm the existence of the rjNN bound state has been 
recently vigorously discussed [166,167]. 

The enhancement factor for the deuteron-i] interaction inferred from 
the data in figure I^T5k varies much stronger in comparison to the proton- 
Tj one, as demonstrated in figure ITT3b . This suggests that the effects of 
this interaction should be even more pronounced in the differential distri- 
butions of the cross section for the pd —>■ pdr] reaction as those observed 
in case of the pp — >■ ppr] process, especially because the "screening" from 
the proton-deuteron interaction is by more than an order of magnitude 
smaller compared to the proton-proton interaction as can be deduced 
from the comparison of the solid lines in figures 14.15b and 14.5b . Experi- 
mental investigations on that issue [162, 168, 169] as well as searches for 
77-mesic nuclei [170] by measuring proton-deuteron induced reactions in 
the vicinity of the rj production threshold are under way. Theoretically, 
the existence of mesic nuclei or quasi-bound meson-nucleus systems is 
not excluded, however, up to now a compelling experimental proof for 
the formation of such a state is still missing. The probability to create 
such states depends crucially on the sign and the strength of the meson- 
nucleus interaction. Since the vr-nucleon and i^^-nucleon interaction 
have been found to be repulsive, it is unlikely that they will form quasi- 
bound states. Contrary, there is evidence for an attractive i^'~-nucleon 
interaction [171], however, corresponding experiments would suffer from 
the low cross sections for the creation of K~ mesons [49]. Furthermore, 
the Coulomb interaction is expected to screen corresponding physical ob- 
servables and might lead to the formation of mesonic atoms bound by the 
electromagnetic interactions. In contradistinction to K~ , the 1] meson is 
uncharged and the observation of an attractive //-nucleon interaction led 
to speculations concerning the existence of ?7-nuclear quasi-bound states. 
In the absence of ?7-meson beams such states bound by the strong inter- 
action would offer a new possibility to study the ?7-nucleon interaction 
since the meson would be trapped for a relatively large time in the nu- 
clear medium. In a very first prediction of such states termed ?7-mesic 
nuclei Haider and Liu [172,173] in the framework of an optical model have 
estimated that it can be formed already for nuclei with atomic numbers 
of A > 12. Further on, the topic was vigorously discussed [149,174-178], 
the hmit was lowered considerably, and presently one considers seriously 
the existence of 77-nucleus bound states even for d, t, ^He, or '^He nu- 
clei [149,176]. 

Experimental evidences for bound states have been found in the near- 
threshold production of r] mesons in the reaction channel pd ^Herj. 
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The unexpected large production amplitude as well as its rapid decrease 
with increasing energy (figure I4.14|) is attributed to a strong s-wave fi- 
nal state interaction associated with a large rj- ^He scattering length 
{a^-me ~ (-2.31 + i 2.57) fm) [176]. 

Assuming that ^He and r] can form a quasi-bound state, it is interesting 
to compare pd — >■ ^He t] and dd — "^He t] production data, where for the 
latter, according to the higher mass number, the formation of a bound 
state should be even more probable. A comparison of corresponding 
near-threshold data is presented in figure 14.141 Interestingly, the slope 
of the fr){p-q) function of the dd -^^Herj reaction is smaller than the one 
of the reaction pd — >^ Her]. This observation suggests, though contra- 
intuition, that the existence of the ?7-mesic nucleus is more probable in 
the case of '^He than ^He. 

4.7 Search for a signal from ry'-proton interaction 
in the invariant mass distributions 

/ am sure that everyone, even among those who 
follow the profession, will admit that everything 
we know is almost nothing compared with what 
remains to he discovered... [281]. 

Rene Descartes 

A phenomenological analysis - presented in section 14.21 - of the de- 
termined excitation function for the pp pprj' reaction (lower part of 
figure 14. 6|) revealed no signal which could have been assigned to the 
proton-?]' interaction. Later on, in section 03 a comparison of the energy 
dependence of the production amplitudes for the pp pprj', pp pprj 
and pp — > ppTT^ enabled to conclude, in a model-free way, that the proton- 
rj' interaction is indeed much weaker than the proton-r] one, and a trial to 
determine quantitatively the scattering length of the proton-?]' potential 
resulted only in a very modest estimation of an upper limit of its real 
part {\Re apri'\ < 0.8 fm) [3]. Seeking for a visible manifestation of the 
proton-?]' interaction we have performed a high statistics measurement 
of pp —>■ pprj' reaction in order to determine a distribution of events over 
the phase space [26]. We expect that the invariant mass spectra for two 
particle subsystems of the pprj' final state can give the first ever experi- 
mental evidence for this still completely unknown interaction. The data 
are presently analyzed [179]. The missing mass spectrum determined 
on-line from 10% of the data written on tapes is shown in figure 14.161 
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From a clear signal visible at the mass of the 77' meson, one can estimate 
that the overall number of registered pp pprj' events amounts to about 
13000. The measurement has been carried out at a beam momentum 
corresponding to the excess energy of Q = 15.5 MeV for the pp pprj' 
reaction. Such value of Q was chosen to enable a direct comparison 
with the invariant mass spectra obtained for the pp pprj reaction (see 
figure I^TTj) ■ without a need for a correction of kinematical factors. 
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Figure 4.16: On-line missing-mass 
distribution of the pp — » ppX re- 
action measured by means of the 
COSY- 11 detection system in Octo- 
ber 2003 [179] at the beam momen- 
tum of 3.257 GeV/c, which in the 
case of the pp — > pp-q' reaction is 
equivalent to the excess energy of 
Q = 15.5 MeV. 



At present due to the need of the subtraction of the unavoidable multi- 
pion background, the available statistics allows only for the background- 
free determination of one-dimensional invariant mass distributions. Yet 
it would be also desired to determine an occupation density over the 
Dalitz-plot since this two-dimensional distribution comprises full empir- 
ically accessible information of the mutual interaction within a three 
particle system. However, the present experimental conditions render 
it impossible, unless the binning was made so large that the possible 
effect would be smeared out completely. In spite of all, we still endeav- 
our to measure this two-dimensional spectrum for both pprj and pprj' 
systems free from the background. A detection of protons and gamma 
quanta from the decay of the meson would allow us to reach this aim. 
Therefore, when the WASA [180] detector will be installed at the cooler 
synchrotron COSY [62] we plan to seize the opportunity of performing 
pertinent investigations [181]. 



5. Dynamics of the near threshold 
production of 77 and vl mesons in 
coUisions of nucleons 

Everything outside of the dynamics is just a verbal 
description of the table of data, and even then the 
data table probably yields more information than 
the verbal description can [182]. 

Werner Karl Heisenberg 



5.1 Comparison of the production yields 

All kinds of reasoning consist in nothing but a 
comparison, and a discovery of those relations, 
either constant or inconstant, which two or more 
objects bear to each other [183]. 

David Hume 

Considerations presented in chapter HI led to the conclusion, that close 
to the kinematical threshold, the energy dependence of the total cross 
section is in first approximation determined via the interaction among 
the outgoing particles and that the entire production dynamics manifests 
itself in a single constant, which determines the absolute scale of the total 
cross section. We demonstrated also that the interaction among created 
particles increases the production probability drastically, though in the 
naive time ordered approach the final state interaction takes place after 
the production. Therefore, the observed increase of the total cross sec- 
tions close to the kinematical threshold, for the pp pprj and pp — *• pprj' 
reactions (see figure . can be considered as a pure quantum mechan- 
ical effect, which clearly demonstrates that the production of particles 
and their mutual interaction must be treated coherently for the proper 
appraisal of the absolute production cross section. Thus, for the sake 
of completeness, we will devote one chapter to a brief description of the 
mechanisms responsible for the production of t] and r]' mesons in the 
collisions of nucleons. 
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As a first step towards the understanding of tlie creation meclianism 
underlying tlie production let us compare the total cross sections for these 
mesons with their flavour-neutral pseudoscalar partner - the meson 7r°. 

Since the masses of these mesons are significantly different^ the influ- 
ence of the kinematical flux factor F on the total cross section and the 
suppression due to the initial state interaction Fjsi depend substantially 
on the created meson. Therefore, for the comparison of the primary dy- 
namics we will correct for these factors and instead of comparing the total 
cross section we will introduce - according to reference [2] - a dimension- 
less quantity a ■ F/Fjsi, which depends only on the primary production 
amplitude Mq and on the final state interaction among the produced 
particles. 

Close-to-threshold the initial state interaction, which reduces the total 
cross section, is dominated by proton-proton scattering in the '^Pq state 
which may be estimated in terms of phase-shifts and inelasticities by 
employing equation ()4.14|) . The Fjsi factor is close to unity for pion 
production and amounts to ~ 0.2 [71] and ~ 0.33 [138] for the rj, and rj' 
meson, respectively, at threshold. 

A comparative study of the production of mesons with significantly 
different masses encounters the difficulty of finding a proper variable at 
which the observed yield can be compared. Seeking such quantity let 
us recall the meaning of the total cross section which is defined as the 
integral over the available phase space volume of transition probabilities 
- reflecting the dynamics of the process - from the initial to the final 
state as written explicitly in equation 13.11 Thus, if the dynamics of the 
production process of two different mesons were exactly the same then the 
above introduced yield would also be strictly the same for both mesons, 
provided it was extracted at the same value of the volume of the phase 
space Vps (eq. I3.20|) . This inference would, however, not be valid if the 
production yields were compared at the 7]m or Q variables. Therefore, 
the volume of the available phase space for the produced particles is the 
most suited quantity for the regarded comparison [2] . This could also be 
the best choice for the investigation of isospin breaking where the cross 
sections for the production of particles with different masses need to be 
compared (e.g. vr+d pprj and 7i~d —>■ nnrj [185]). Figure l^!T] shows 
the yield of vr*^, t], and t]' mesons in the proton-proton interaction as a 
function of the available phase space volume. 



1 The TT^, r), and r]' masses amount to 134.98 MeV/c^, 547.30 MeV/c^, and 957.78 MeV/c^, re- 
spectively [43]. 
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Figure 5.1: Total cross section (cr) multiplied by the flux factor F and divided by the 

initial state interaction reduction factor Fisi versus the available phase space volume for the 
reactions pp — > pprj (squares [10, 18,50-53, 121]), pp — > ppn° (diamonds [114-117, 184]), and 
pp^pprj' (circles [3,5,11,53,97]). 

The onset of higher partial waves is seen for tt" and 7] mesons in the 
Vps range between 10^ and lO^MeV^, whereas the whole range covered 
by the r)' data seems to be consistent with the pure Ss production. One 
can also recognize that the data for the Ss final state are grouped on 
parallel lines indicating a dependence according to the power law a ■ 
^/Fisi ~ a ■ Vpg^^ [186] and that over the relevant range of Vpg the 
dynamics for rj' meson production is about six times weaker than for the 
7r° meson, which again is a further factor of six weaker than that of the rj 
meson. This is an interesting observation, since the quark wave functions 
of T} and Tj' comprise a similar amount of strangeness (~ 70% [237]) 
and hence, in the nucleon-nucleon collision one would expect both these 
mesons to be produced much less copiously than the meson 7r° being 
predominantly built out of up and down quarks. On the hadronic level, 
however, one can qualitatively argue that the r] meson - contrary to 7r° 
and 77' - owes its rich creation to the existence of the baryonic resonance 
A'^*(1535) whose branching ratio into the Nr) system amounts to BO- 
SS % [43]. There is no such established resonance, which may decay 
into an s-wave 77' system [43], and the vr*^ meson production with the 
formation of the intermediate A(1232) state is strongly suppressed close- 
to-threshold, because of conservation laws. 



70 



Hadronic interaction of rj and r/' mesons with protons 



5.2 Mesonic degrees of freedom 

For what is determinate cannot have 
innumerable explanations [187]. 

Nicholas Copernicus 

One of the most crucial issue in the investigations of the dynamics of the 
close-to-threshold meson production is the determination of the relevant 
degrees of freedom for the description of the nucleon-nucleon interaction, 
especially in case when the nucleons are very close together^. The tran- 
sition region from the hadronic to constituent quark degrees of freedom 
does not have a well defined boundary and at present both approaches 
are evaluated in order to test their relevance in the description of close- 
to-threshold meson production in the collision of nucleons. Simple geo- 
metrical considerations presented in figure indicate that at distances 
smaller than 2fm the internucleon potential should begin to be free of 
meson exchange effects and may be dominated by the residual colour 
forces [189]. 



Figure 5.2: A cartoon illustrating 
in naive geometrical terms that for 
r < 2rjv + 2rM meson exchange is un- 
likely to be appropriate for the description 
of the internucleon potential. Figure and 
caption are taken from reference [189]. 




In section it was shown, that the close-to-threshold production of 
mesons occurs when the colliding nucleons approach distances of about 
0.5 fm in case of 7r° and of about 0.18 fm in case of production (see ta- 
ble |^^. This distance is about one order of magnitude smaller than 2 fm 
and it is rather difficult to imagine - in coordinate space - an exchange 



^ These investigations are listed as one of the key issues in hadronic physics [188]. 
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of mesons between nucleons as mechanism of the creation process. Thus, 
the meson exchange for the threshold meson production can be under- 
stood as an effective description of the process occuring on the deeper 
level. Such small collision parameters imply that the interacting nucleons 
- objects of about 1 f m - overlap and their internal degrees of freedom 
may be of importance. On the other hand, as demonstrated by Han- 
hart [81], if there is a meson exchange current possible at leading order 
it should dominate the creation process. This can be derived from the 
strong dependence of the production operator on the distance between 
the colliding nucleons. 
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a) b) 
Figure 5.3: Illustration of the momentum mismatch. The picture shows why one should 
a priori expect meson exchange currents to be very important close to the threshold. The 
figure and title are adapted from reference [81]. By courtesy of C. Hanhart. 
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Figure 5.4: Diagrams for the pp — » pp Meson reaction near threshold: (a) 

Meson-brcmsstrahlung (nucloonic current) (b) "rcscattcring" term (nucloonic current) (c) 
production via heavy-meson-exchange (d) emission from virtual meson (mesonic current) 
(e) excitation of an intermediate resonance (nucleon resonance current). 
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This particular feature of the meson exchange current will prevent the 
cancellation of the transition matrix element being a convolution of the 
rapidly and mildly oscillating initial and final state wave functions, as 
it is illustrated in figure 15.31 Consequently, during the last decade, the 
effective theory based on meson exchanges, which accounts for the size 
of the participating particles by introduction of the momentum transfer 
dependent form factors, has been extensively employed for the description 
of the creation process. Figure l^^ represents the mechanisms in question. 
In the next sections we will report on the stage of understanding of the 
rj and 7]' mesons creation on the hadronic level, and thereafter we shall 
briefly present first attempts to explain the meson threshold production 
in the domain of quarks and gluons. 

5.3 Baryonic resonance — a doorway state for the 
production of the meson r] 

...although these excitations of the soul are often 
joined with the passions that are like them, they 
may also frequently be found with others, and may 
even originate from those that are in opposition to 
them [190]. 

Rene Descartes 

It is at present rather well established [67, 72, 124, 191-196] that the r] 
meson is produced predominantly via the excitation of the Sn baryonic 
resonance iV*(1535) which subsequently decays into t] and nucleon, and 
whose creation is induced through the exchange of vr, rj, p, a, and u 
mesons, as shown in figure 15.4b . Although all the quoted groups re- 
produce the magnitude of the total cross section, their models differ 
significantly as far as the relative contributions from the tt, r], and p ex- 
change mechanisms are concerned. The discrepancies are due to the not 
well known strength of Meson-N-Su couplings and the r]N scattering 
potential. For example, while the dominance of the p meson exchange 
is anticipated by authors of references [67,124,191-193], it is rather the 
exchange of the rj meson which dominates the production if one takes 
into account the effects of the off-shell rjN scattering [198], or uses the 
multi-channel multi-resonance model [72]. In any case, the hitherto per- 
formed studies aiming to describe the total cross section show that the 
close-to-threshold production of rj mesons in nucleon-nucleon collisions 
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is dominated by the intermediate virtual Sn nucleon isobar whose width 
overlaps with the threshold. 

Figure 5.5: Total cross sections for the 
pp ppr] (upper panel) and pn pnrj 
(lower panel) reactions as a function of 
excess energy. The dashed curves corre- 
spond to the nucleonic current contribu- 
tion and the dash-dotted curves to the 
mesonic current consisting of rjpp, 'quooj, 
and TjaaTT contributions. The resonance 
current presented by the dotted line con- 
sists of the predominant Sii(1535) and of 
the Pii(1440) and Di3(1520) resonances 
excited via exchange of tt, rj, p, and uj 
mesons. The solid curves are the to- 
tal contribution. The deviation from the 
data at low Q in the upper panel re- 
flects the r]p FSI which was not included 
00 in the calculations. The data are from 
refs. [10,50-53,197]. The figure is adapted 
from [100]. 

In order to disentangle the various scenarios of the Sn excitation a 
confrontation of the predictions with other observables is needed. The 
interference between considered amplitudes causes a different behaviour - 
depending on the assumed scenario - e.g. of the r] meson angular distribu- 
tions. These differences, however, are too weak in the close-to-threshold 
region to judge between different models. Also the ratio of the r] meson 
production via the reactions pp —>■ pprj and pn pnrj can be equally 
well described by either assuming the p meson exchange dominance [193] 
or by taking pseudoscalar and vector mesons for exciting the Sn reso- 
nance [195]. In the latter case, shown in figure 1^31 the excitation of the 
resonance via the p meson exchange was found to be negligible. Yet, 
promisingly, the predictions of the analyzing power depend crucially on 
the assumed mechanism [193,195]. This fact has already triggered ex- 
perimental investigations which aim to determine the spin observables. 
First results of tentative measurements and perspectives of these studies 
will be presented in section 15.41 
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5.4 Spin degrees of freedom — a tool to study 
details of the pp pprj reaction dynamics 

/ thought that scientific theories were not the 
digest of observation, hut that they were inventions- 
conjectures boldly put forward for trial, to he 
eliminated if they clashed with observation; with 
observations which were rarely accidental but as a 
rule undertaken with definite intention of 
testing a theory by obtaining, if possible, a 
decisive refutation [183]. 

David Hume 



In spite of the precise measurements of the total cross section for the 
creation of rj meson in proton-proton [10,18,50-53,121] as well as proton- 
neutron [197] collisions there are still many ambiguities in the description 
of the mechanism underlying the production process. Calculations per- 
formed under different - often mutually exclusive - assumptions led to an 
equally satisfying description of the data. As we had already described in 
section EH it is generally anticipated [67, 72, 191, 192, 194, 196, 198-200] 
that the rj meson is produced predominantly via the excitation of the 
Sii baryonic resonance iV*(1535), whose creation is induced through the 
exchange of the virtual vr, rj, p, a, and u mesons, however, at present it 
has still not been established what the relative contributions originating 
from a particular meson are. Measurements of the total cross section in 
different isospin channels add some more limitations to the models, yet 
still the rj meson production in the pp pprj and pn pnrj reactions 
can be equally well described by e.g. assuming the p meson exchange 
dominance [193] or by taking contributions from the pseudoscalar meson 
exchanges [195]. Therefore, for a full understanding of the production 
dynamics the determination of more selective observables is mandatory. 
A natural choice seems to be an extension of the research to encompass 
the spin degrees of freedom. Figure 15.61 shows that predictions of the 
analyzing power Ay are strongly sensitive to the type of the exchanged 
particle. A discrepancy between considered models is so pronounced 
that even a measurement of the beam analysing power with a preci- 
sion of ±0.05 could exclude at least one of the above mentioned 
possibilities with a statistical significance better than three standard 
deviations. A tentative measurement of that quantity was conducted [7] 
and points in figure ESb indicate the obtained result. Unfortunately, the 
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present statistics is insufficient to allow discrimination between the alter- 
natives considered, and for a conclusive inference a better accuracy of the 
data is required. After successful measurement of Ay for the pp — > ppr) 
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Figure 5.6: Predictions for the angular dependence of the analysing power of the reaction 
pp pprj at Q = 15 (a) and at Q = 37MeV (b) compared to the data taken at 
Q = AQ MeV [7] . Solid lines present results of the model [195] characterized by the dominance 
of the exchange of the pseudoscalar mesons in the production process. Dashed lines show 
expectation obtained under the assumption that the production mechanism is predominated 
by the exchange of the p meson [193]. 

reaction at Q = 40 IVIeV, we continued the study at Q = 37 IVIeV and at 
Q = 10 MeV [27,201,202]. The measurement at Q = 37 MeV was per- 
formed with proton beam polarisation amounting to approximately 70 %, 
a value significantly larger than in the first measurement. Moreover, the 
luminosity integrated over the measurement period was larger by about 
a factor of 1.5. These two factors together, improve the accuracy of 
the measurement more than twice, in the sense, that the errors corre- 
sponding to those in figure 15.6b are expected to be two times smaller. 
The second measurement was performed at Q = 10 MeV since at this 
excess energy the maximum of the discrepancy between the predictions 
of the regarded models [193, 195] is expected. The data on the angular 
distributions of the analysing power would also enable to determine the 
relative magnitudes - or at least to set upper limits - for the contribution 
from the higher than s-wave partial waves to the production dynamics. 
This can be done with an accuracy by far better than the one resulting 
from the measurements of the distributions of the spin averaged cross 
sections. This is because the polarisation observables are sensitive to 
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the interference terms between various partial amplitudes, which may 
become measurable even if one of the interfering terms alone appears to 
be insignificant in case of the spin averaged cross sections. The analy- 
sis would be based on the formalism developed for the measurement of 
the tt" meson production with a polarized proton beam and target [78], 
which we have already used for the analysis of our first measurement 
with a polarized proton beam [7, 203, 204] . This procedure permits the 
derivation of all available information contained in the distributions of 
the cross sections and the analysing power. Information about the par- 
tial wave decomposition is of crucial importance for the interpretation of 
the production dynamics but also for the study of the mutual interaction 
of the produced T^-proton-proton system. The strength of the proton-r^ 
interaction cannot be derived univocally from the differential cross sec- 
tion distributions without the knowledge of the contribution of various 
partial waves. 

Measurements with the polarized beam and target would allow de- 
termination of these shares. As demonstrated in references [78, 91] the 
close-to-threshold contributions of the Ps and Pp partial waves can be 
determined in the model-free way from the measurements of the spin 
dependent total cross sections only. For example the strength of the Ps 
final state can be expressed as [78]: 



where atot denotes the total unpolarized cross section and Acxr and Actl 
stand for differences between the total cross sections measured with anti- 
parallel and parallel beam and target polarizations. Subscripts T and L 
associate the measurements with the transverse and longitudinal polar- 
izations, respectively. 

At present none of the experimental facilities enable to investigate the 
pp — > pprj reaction. However, an installation of the WASA detector at 
COSY will offer possibility to accomplish such studies. Corresponding 
measurements have been already proposed [205], and the predictions of 
the spin correlation functions are given e.g. in references [79,206]. 




(5.1) 
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5.5 Possible mechanisms responsible for the 
creation of the rl meson 

But in fact, even if all writers were honest and 
plain; even if they never passed off matters of 
doubt upon us as if they were truths, but set 
forth everything in good faith; nevertheless, since 
there is hardly anything that one of them says 
but someone else asserts the contrary, we should 
be continually uncertain which side to believe. It 
would be no good to count heads, and then follow 
the opinion that has most authorities for it; for 
if the question that arises is a difficult one, it 
is more credible that the truth of the matter 
may have been discovered by few men than by 
many [207]. 

Rene Descartes 



In case of the r]' meson the investigations of the mechanisms underly- 
ing the production process are even less advanced than these of rj and tt. 
This is partly due to the fact that there is no well established baryonic 
resonance decaying into rj'N channel, which could constitute a door- 
way state for the rj' production, and hence there is no mechanism whose 
strength could be regarded as dominant allowing for the neglect of other 
possibilities at least in the first order. There is also not much known 
about the rj'NN and other relevant coupling constants, making impos- 
sible inferences about the relative strength of different mechanisms on 
the level achieved in the understanding of the tt meson creation [81,208]. 
It is also partly because the experimental database is much poorer than 
this for rj and tt mesons. This holds not only for the data on the r)' 
production in the collisions of nucleons, but also for the data at relevant 
energies on nucleon-nucleon elastic scattering. The latter are needed for 
the estimation of the reduction of the production cross section due to 
the interaction between nucleons in the initial state. It must be stated, 
however, that in the recent years the database - relevant for the pro- 
duction of the rj' meson - for the proton-proton reactions [3,53,97,209] 
improved significantly, yet measurements with a corresponding accuracy 
for the proton-neutron channel are still missing. 
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5.5.1 Meson exchange models 

Salviati. But if, of many computations, not even two 
came out in agreement, what would you think of that? 

Simplicio. // that is how matters stand, it is truly a 
serious defect [139]. 

Galileo Galilei 

For the r] meson case a contribution from the mesonic current where 
the meson is created in the fusion of virtual e.g. p or u mesons emitted 
from both colliding nucleons is by a factor of thirty weaker in comparison 
to the overwhelming strength of the resonance current f figure . In 
contrast, the mesonic current suffices to explain the magnitude of the 
close-to-threshold rj' meson production in the proton-proton interaction, 
which is just by about a factor of thirty smaller compared to the rj meson. 
Among the mesonic currents regarded by authors of reference [138] the 
ppr]' gives the dominant contribution, by a factor of five larger than the 
arjri'- and uurj'-exchange. However, the understanding of the r]' produc- 
tion on the hadronic level is far from being satisfactory. The magnitude 
of the total cross section was also well reproduced in the frame of a 
one-boson-exchange model (nucleonic current) where the virtual boson 
{B = 7i,r],a, p,u,ao) created on one of the coUiding protons converts to 
the r]' on the other one [69]. Taking into account the off-shell effects 
of the Bp Tj'p amplitude it was found that the short range a and p 
meson exchanges dominate the creation process, whereas the tt exchange 
plays a minor role. On the contrary, other authors [210] reproduced the 
magnitude of the total cross section with vr exchange only and concluded 
that the 77, p, and uj exchange currents either play no role or cancel each 
other. However, in both of the above quoted calculations [69,210] the ini- 
tial state interaction between protons, which reduces the rate by a factor 
of about 3, was not taken into account and hence the obtained results 
would in any case reproduce only 30 % of the entire magnitude of the 
cross section and could be at least qualitatively reconciled with the men- 
tioned result of reference [138], where the nucleonic current was found to 
be small. Moreover, the choice of another prescription for the form fac- 
tors could reduce the one-pion-exchange contribution substantially [138]. 
However, the picture that the rj' meson is predominantly created through 
the mesonic current, remains at present unclear as well. This is because 
the magnitude of the total cross section could have also been described 
assuming that the production of rj' is resonant [138]. As possible interme- 
diary resonances the recently reported [211] Sii(1897) and Pii(1986) have 



5. Dynamics of the near threshold production of rj and r/' mesons 



79 



been considered. These resonances were deduced from i]' photoproduc- 
tion data, under the assumption that the close-to-threshold enhancement 
observed for the 7p rj'p reaction can be utterly assigned to resonance 
production. Further, as well strong, assumptions have been made in the 
derivation of the gNN*rj' and gNN*n coupling constants [138]. 

Hence, it is rather fair to state that in the case of the close-to-threshold 
TT^ and 1] meson production in nucleon-nucleon collisions the dynamics is 
roughly understood on the hadronic level, but the mechanisms leading 
to the f]' creation are still relatively unknown. 

Until now it has not been possible to satisfactorily estimate the rel- 
ative contributions of the nucleonic, mesonic, and resonance current to 
the production process. In fact, model uncertainties allow that each 
one separately could describe the absolute values of the pp pprj' total 
cross section. This rather pessimistic conclusion calls for further exper- 
imental and theoretical research. The understanding of the production 
dynamics of the rj' meson on both hadronic and the quark-gluon level 
is particularly important since its wave function comprises a significant 
gluonic component [212], distinguishing it from other mesons and hence 
the comprehension of the mechanism leading to its creation in collisions 
of hadrons may help to determine its quark-gluon structure. That, in 
turn will be helpful when investigating possible glueball candidates [212]. 
Hereafter we will briefly report of the gluonic mechanisms which may - 
along with the meson exchange processes discussed above - contribute 
to the T]' production. 

5.5.2 Approaches on the quark-gluon level 

We must clearly acquire knowledge of factors that are 
primary. For we claim to know a thing only when 
we believe that we have discovered what primarily 
accounts for its being [90]. 

Aristotle 

As afore mentioned, the close-to-threshold production of rj and rj' 
mesons in the nucleon-nucleon interaction requires a large momentum 
transfer between the nucleons and hence can occur only at distances of 
about 0.3 fm (see table ITT|1 . This suggests that the quark-gluon de- 
grees of freedom may indeed play a significant role in the production 
dynamics of these mesons. A possibly large glue content of the rj' and 
the dominant fiavour-singlet combination of its quark wave function may 
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cause that the dynamics of its production process in nucleon-nucleon col- 
hsions is significantly different from that responsible for the production 
of other mesons. In particular, the i]' meson can be efficiently created via 
a "contact interaction" from the glue which is excited in the interaction 
region of the colliding nucleons [44]. A gluon-induced contact interac- 
tion contributing to the close-to-threshold pp — ^ ppr]' reaction derived in 
the frame of the U(l)-anomaly extended chiral Lagrangian is discussed 
in references [44,213,214]. The strength of this contact term is related 
to the amount of spin carried by polarized gluons in a polarized pro- 
ton [214,215], thus making the study of the close-to-threshold rj' meson 
production even more interesting. 
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Figure 5.7: Diagrams depict- 
ing possible quark-gluon dynam- 
ics of the reaction pp pprj' ■ 
(a) production via a fusion of 
gluons [216] with rearrangement 
of quarks, (b) production via 
a rescattering of a "low energy 
pomeron" . 
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Figure 5.8: Coupling of rj and ri' to two gluons through (a) quark and antiquark triangle 
loop and (b) gluonic admixture, rji denotes the flavour-singlet quark-antiquark state. This 
indicates that gluons may convert into the 77 or rj' meson via a triangle quark loop only by 
coupling through their flavour singlet part. The flgure is taken from reference [217]. 



Figure 15.71 depicts possible short-range mechanisms which may lead 
to the creation of the t]' meson via a fusion of gluons emitted from the 
exchanged quarks of the colliding protons [216] or via an exchange of 
a colour-singlet object made up from glue, which then re-scatters and 
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converts into i]' [218]. The hadronization of gluons to the rj' meson 
may proceed directly via its gluonic component or through its over- 
whelming flavour-singlet admixture rji (see fig. 15. Sj) . Contrary to the 
significant meson exchange mechanisms and the fusion of gluons of fig- 
ure EIZI graph a), the creation through the colour-singlet object proposed 
by S.D. Bass (graph EHb) is isospin independent, and hence should lead 
to the same production yield of the rj' meson in both reactions {pp — *• pprj' 
and pn — > pnr]') because gluons do not distinguish between flavours. 
This property should allow to test the relevance of a short range glu- 
onic term [219] by the experimental determination of the cross section 
ratio Rrj' = a{pn pnri')/a{pp ppf]')y which in that case should be 
close to unity after correcting for the final and initial state interaction 
between participating baryons. The other extreme scenario - assum- 
ing the dominance of the isovector meson exchange mechanism - should 
result in the value of i?^/ close to 6.5 as had already been established 
in the case of the t] meson [197]. Perspectives of the experimental in- 
vestigations aiming to determine i?^/ and the discussion to what extent 
these studies may help to establish contributions from quarks and glu- 
ons in the 1]' meson will be presented in section 15.61 In addition to the 
investigation presented above [44,213], the interesting features of the 
close-to-threshold meson production, in particular the large cross section 
for the T) meson in proton-proton interactions exceeding the one of the 
pion, or even more surprisingly large cross section of the rj production 
in proton-neutron collisions, encouraged also other authors to seek the 
underlying - OZI rule violating - creation mechanisms in the frame of 
microscopic models of QCD [220-222]. The hitherto regarded processes 
are presented in figure 15.91 As indicated in the pictures the structure of 
participating baryons has been modeled as quark-diquark objects with 
harmonic confinement [223]. In the upper graphs the large momentum 
transfer is shared by the exchanged gluons with a subsequent interchange 
of quarks to provide a colourless object in the final state. The lower 
graphs depict the two examples of instanton induced interactions with a 
6-quark-antiquark ffig. I5.9l i) and two-gluon vertex ffig. 15.9b ). Adjusting 
the normalization to the cross section of the pp ppir^ reaction at a sin- 
gle energy point the model [220] accounts roughly for close-to-threshold 
cross sections of other pseudoscalar and vector mesons in proton-proton 
collisions. 
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Figure 5.9: Diagrams for the qq q<l{q<j) production operator: (a), (b) Two-gluon 
exchange and rescattering mechanism (c) Correlated, colourless two-gluon exchange. The 
dashed line indicates the excitation of an intermediate A'^A'^* system. (d) Instanton induced 
6-quark interaction. Figures (a-d) according to [220]. (e) The instanton contribution to 
the 77 meson production in the proton- neutron interaction. Figure copied from [221]. 



Though this approach is characterized by a significantly smaller num- 
ber of parameters than in the meson exchange models, their uncertainties 
allow for the description of the data equally well, either by the gluon ex- 
change or by the instanton induced interactions, at least in case of the tt", 
T]', u, and mesons. Yet for the resonance dominated 1] and meson 
production in proton-proton collisions it was found that the instanton 
induced interaction presented by graph is not sufficient. Similarly, 
the authors of reference [221] argue that the instanton induced interac- 
tion with a quark-gluon vertex (graph IHTHfe ) should be of no importance 
for the T) production in proton-proton collisions. This arises from the 
properties of the vertex which lead to the 1] production only in case of an 
interaction between quarks of different flavours [ud — *■ udgg) and in the 
quark-diquark model of baryons the proton consists of wd-diquark and 
u-quark and correspondingly the neutron is modeled as wrf-diquark and 
(i-quark. Although negligible in case of proton-proton interactions, this 
mechanism may contribute signiflcantly to the total cross section of the 
pn — *■ pnr] reaction and indeed as shown in reference [221] it reproduces 
the data. However, the magnitude of the cross section is very sensitive 
to the size of the instanton in the QCD vacuum and - similar to the 
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uncertainty of coupling constants in case of the hadronic approach - at 
present it makes precise predictions impossible. Thus, despite the partial 
successes, at the present stage of developments both approaches - on the 
quark-gluon and on the hadronic level - do not provide an unambiguous 
answer to the dynamics of the close-to-threshold meson production in 
the nucleon-nucleon interaction. 

5.6 Exploration of isospin degrees of freedom 

When a person is placed between two choices, the 
person bends himself or herself toward the choice he 
or she desires [224] ■ 

Hildegard von Bingen 

Treating proton and neutron as different states of nucleon distinguished 
only by the isospin projection, +| for the proton and — | for the neu- 
tron, we may classify the A^A^ NNX reactions according to the total 
isospin of the nucleons pair in the initial and final state. A total isospin 
of two nucleons equals 1 for proton-proton and neutron-neutron pairs, 
and may acquire the value of 1 or for the neutron-proton system. 
Since 77 and 77' mesons are isoscalars, there are only two pertinent tran- 
sitions for the NN —>■ NNX reaction, provided that it occurs via the 
isospin conserving interaction. These are CToo ^-nd an, where the first and 
second subscript - labeling the total cross section - describe the total 
isospin of the nucleon pair in the initial and final state, correspondingly. 

It is thus enough to measure two reaction channels for an unambiguous 
determination of isospin and 1 cross sections. In particular for the total 
cross section the following relations are satisfied: 

and 

Cpn-^pnr]' — ^('^00 + Cn). (5-3) 

In the case of the proton-neutron reaction the total cross sections add in- 
coherently since the different isospin of the NN pairs imphes that either 
total spin or angular momentum of these pairs must differ. This is be- 
cause the A^A^ system is bound to satisfy the Pauli Principle. Therefore, 
different isospin states do not interfere as far as the total cross section is 
concerned [81]. 
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It is needless to mention, that experimental determination of cross 
sections for different isospin configurations allows for deeper insight into 
the production mechanism. In sections 16.41 and 16.61 we will present how 
one can conduct the study of the meson production in the proton-neutron 
and neutron-neutron reactions and in the next section we will give an 
example of inferences which may be drawn from the comparison of the 
total cross sections for the pp — > pprj' and pn pnt]' reactions. 

5.6.1 Glue content of the t]' meson 

... even God does not know his own nature [225]. 

Bryan Magee 

Experience no doubt teaches us that this or that 
object is constituted in such and such a manner, but 
not that it could not possibly exist otherwise [32]. 

Immanuel Kant 

The most remarkable feature - in the frame of the quark model - dis- 
tinguishing the r]' meson from all other pseudscalar and vector ground 
state mesons, is the fact, that the r]' is predominantly a flavour- singlet 
combination of quark- ant iquark pairs and therefore can mix with purely 
gluonic states. A comprehensive discussions on the issue of how one 
can experimentally determine the differences in the quark-gluon struc- 
ture between t]' and other mesons can be found in the proceedings of 
the International Workshop on the Structure of the r]' Meson held in Las 
Cruces [226]. Unfortunately, except for the gedanken experiments, no ex- 
perimentally feasible test allowing for the model-independent conclusions 
was suggested. 

Due to the short life-time of this meson it is impossible to use it as 
a beam or target. This fact entails that the study must be performed 
by measurements of the reactions where the r]' meson is created in the 
collisions or decays of more stable particles, and preferentially the studied 
observables should depend - in a predictive manner - on properties of 
the 7]' meson. 

In this section we will argue that a comparison of the close-to-threshold 
total cross sections for the rj' meson production in both the pp — ^ pprj' 
and pn — > pnr]' reactions should provide insight into the flavour- singlet 
(perhaps also into gluonium) content of the rj' meson and the relevance 
of quark-gluon or hadronic degrees of freedom in the creation process. It 
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is of great credit to S. D. Bass for this idea to be shared with us [218]. 
However, prior to the formulation of the main conjectures we will recall 
the most important facts of the production mechanism necessary for the 
understanding of the further conclusions. 

Close-to-threshold production of rj and rj' mesons in the nucleon- 
nucleon interaction requires a large momentum transfer between the nu- 
cleons and occurs at distances in the order of ~0.3 fm. This implies 
that the quark-gluon degrees of freedom may play a significant role in 
the production dynamics of these mesons. Therefore, additionally to the 
mechanisms associated with meson exchanges it is possible that the rj' 
meson is created from excited glue in the interaction region of the collid- 
ing nucleons [44,214], which couple to the r]' meson directly via its gluonic 
component or through its SU(3)-flavour-singlet admixture. The produc- 
tion through the colour- singlet object as suggested in reference [44] is 
isospin independent and should lead to the same production yield of the 
1]' meson in the pn — >■ pn gluons pnrj' and pp — > pp gluons — pprj' re- 
actions after correcting for the final and initial state interaction between 
the nucleons. 

Investigations of the T^-meson production in collisions of nucleons al- 
lowed to conclude that, close to the kinematical threshold, the creation of 
rj meson from isospin 1 = exceeds the production with I = 1 by about a 
factor of 12. This was derived from the measured ratio of the total cross 
sections for the reactions pn pnrj and pp pprj [R^ = f^^^^^), 
which was determined to be i?^ ~ 6.5 in the excess energy range between 
16 MeV and 109 MeV [197]. The large difference of the total cross sec- 
tion between the isospin channels suggests the dominance of isovector 
meson (tt and p) exchange in the creation of rj in nucleon-nucleon colli- 
sions [197,227]. Since the quark structure of r] and rj' mesons is very 
similar we can - by analogy to the rj meson production - expect that in 
the case of dominant isovector meson exchange the ratio Rr^i should also 
be about 6.5. If, however, the rj' meson was produced via its flavour-blind 
gluonic component from the colour-singlet glue excited in the interaction 
region, the ratio i?^/ should approach unity after corrections for the in- 
teractions between the participating baryons. 

Figure EUni demonstrates qualitatively the fact that the production of 
mesons in the proton-neutron collisions may be more probable than in the 
proton-proton interaction if it is driven by the isovector meson exchanges 
only. This is because in the case of the proton-neutron collisions there are 
always more possibilities to realize an exchange or fusion of the isovector 
mesons than in the case of the reaction of protons. More precisely, for the 
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Figure 5.10: (left) Example of diagrams with the isovector (upper row) and isoscalar 
(lower row) meson exchange leading to the creation of the meson rj' in the proton-proton 
and proton-neutron collisions, (right) Fusion of the virtual uj (isoscalar) and p (isovector) 
mesons emitted from the colliding nucleons. 

one meson exchange diagrams, when considering vr, 77, p, and u mesons 
and neglecting initial and final state interaction, the matrix element has 
the following structure [227]: 



\M{pn pnr])\^ = - [\T^ + T^-T^~ T^p + |3T^ - T, + 3T, - T^p] 



\M{pp ^ ppr])\^ = \T^ + T^-T,- (5.4) 

where the amplitudes are labeled by the names of the exchanged mesons. 
In general, when regarding the isospin structure of the exchange cur- 
rent, as shown in reference [81], if the production operator accounts for 
the isoscalar meson exchange then aoo = and in case of isovector ex- 
change (Too = 9(Jii. This is also clear from the particular case represented 
by equations 15.41 when combined with formulae 15.21 and 15.31 Exploiting 
these equations we can estimate that, if the production of rj' meson via 
the pp ppr)' was governed only by the exchange of the isovector mesons 
the ratio i?^/ would be equal to five, and in the case of isoscalar mesons 
it would be equal to unity. 

A very important result of the theoretical investigation, relevant for 
further consideration, is that regardless of whether it is a mesonic, nucle- 
onic, or resonance current the contribution from the exchange of isovec- 
tor mesons (p or tt) is much larger from that of isoscalar ones (u; or 
rj) [138,141,227]. This unequivocally entails that if the ratio i?^/ - cor- 
rected for FSI and ISI distortions - will be found to be close to unity we 
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will have a clear indication that the rj' is produced directly by gluons. 
Gluons, as we discussed in chapter may hadronize to rj' either via its 
SU{3)f flavour-singlet component or via its gluonic content. In order to 
disentangle these two effects a substantial theoretical input is required. 
The beneflt of the invested effort will be the determination of the quark- 
gluon structure of the t]' meson. If, however, the measured ratio i?^/ will 
not be equal to one, a quantitative determination of the contribution 
from gluonic mechanism to the production process will require a better 
understanding of the meson exchange currents. Then if the contribution 
of meson exchange currents is once understood we will also be able to 
infer the one from gluons. The dynamics of the meson production in 
both hadro- and photo-production is at present vigorously studied. The 
recent results are reported e.g. in references [4,81,228,229]. 

The close-to-threshold excitation function for the pp — pprj' reaction 
has been already determined [3,11,53,97] whereas the total cross section 
for the rj' meson production in the proton-neutron interaction remains 
unknown. As a flrst step towards the determination of the value of Rrj' 
the feasibility of the measurement of the pn — »• pnrj' reaction by means of 
the COSY-11 facility was studied using a Monte-Carlo method [20,230]. 
As a second step, a test experiment of the pn —> pnr] reaction - suspected 
to have by at least a factor of thirty larger cross section than the one for 
the pn —y pnr]' reaction - was performed. In this test measurement, using 
a beam of stochastically cooled protons and a deuteron cluster target, we 
have proven the ability of the COSY-11 facility to study the quasi-free 
creation of mesons via the pn pnX reaction. Appraisals of simula- 
tions [20, 230] and preliminary results of the tentative measurements of 
the quasi-free pn — ^ pnr] reaction performed using the newly extended 
COSY-11 facility [21,25,28] will be presented in section IFTHl 

5.7 Transition probability as a function of 
particles' "virtuality" 

Es gibt keine freie Masse, genauso wie es kein leeres 
Vakuum gibt [231]. 

Kurt Kilian 

First of all the title of this section requires apology and the explanation 
what is meant under the notion of virtuality. Usually we consider and 
measure the reactions of free particles for which the difference between 
the squared energy and momentum is equal to the square of its mass. 
This is not the case for the nucleons bound inside the nuclei, whose 
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masses may differ from the mass of the free particles. In particular, in 
case of the deuteron, the entire difference (summed for both nucleons) 
amounts to at least a value equivalent to the binding energy. The inner 
motion of nucleons inside a nucleus makes this difference even larger. 
Thus principally the virtual nucleons inside a deuteron may not be iden- 
tical with their free equivalent. By virtuality we would like to name the 
difference between the mass of the real particle and the mass calculated 
as a difference between squared energy and squared momentum of its 
virtual counterpart. When using a deuteron as a source of neutrons for 
the measurement of e.g. proton-neutron reactions this difference is usu- 
ally regarded as an obstacle making difficult a direct comparison of free 
and quasi-free scattering. Here we would like to suggest how to take 
advantage of this fact. Namely we intend to make a systematic study of 
the dependence of the cross sections for a meson production in nucleon- 
nucleon collisions as a function of their virtuality, which in the first order, 
if a deuteron is used as a target, can be deduced from the momentum of 
the nucleon which does not take part in the reaction. 

There is a difference between the free and bound nucleon since they 
have different masses. Hence it is natural to expect that they will also 
differ in other aspects. Recently a chance appeared to realise such studies 
at the WASA facility [180], which is planned to be installed at COSY in 
the near future, and the corresponding letter of intent has been already 
prepared [181]. 

5.8 Natural width of mesons — does it depend on 
the momentum transfer? 

Throughout the history of spectroscopy, 
improved resolution has led to discoveries 
of finer and finer structures [232]. 

Bogdan Castle Maghch 

Although only loosly connected to the subject of this monograph, we 
would like to comment on a very recent and highly interesting hypothe- 
sis announced by Maglich at the HADRON'03 conference in Aschaffen- 
burg [232]. Based on the experience gained from measurements of the 
Oi and 02 mesons via the 7r~p — > pX~ reactions, he postulated that the 
width of mesons may be a function of the momentum transfer |t|, and 
in particular that the mesons width grows with decreasing value of \t\. 
According to his idea the intrinsic widths of all meson states are narrow 
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and can be observed only at high momentum transfers \t\ > 0.2. Since 
the production cross sections are inversely proportional to \t\, the ex- 
perimental data samples are dominated by events corresponding to low 
values of the momentum transfer, and therefore the observed widths of 
mesons appear to be broad [232] ^. 

Nothing can be a priori excluded, however we would like to indicate 
that the broadening of the width of mesons may only be an apparent 
effect caused by the error propagation through the formula exploited for 
the missing mass calculation. Such effects have been observed at the 
missing mass spectra from measurements of the rj meson at the COSY- 
11 [1] and TOF [95] facilities. 
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Figure 5.11: Distribution of the polar scattering angle of the X system created via the 
pp — > ppX reaction measured at Q = 15.5 MeV above the threshold for the production of 
the T] meson [18]. 



Figure 15.111 presents the dependence of the missing mass distribution 
as a function of the polar angle of the emission of the produced system X. 
At a mass value corresponding to the mass of the rj meson one observes 
a strip of enhanced population-density over a smooth multi-pion back- 
ground (figure inHH from section 1^7^ constitutes the projection of this plot 
onto a mass axis). Here it is incontestably evident that the width of the 
signal originating from the reaction pp pprj is not constant but rather 
varies with the polar angle G*. Specifically, the signal from rj meson is 



•^An alternate explanation proposed after completion of this treatise can be found in refer- 
ence [233]. 
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much sharper for cosG^ = 0.3 than for cosQ* = -0.95. The difference 
cannot be assigned to the variation of the width of the t] meson, since 
this is three orders of magnitude smaller (F^ = 1.18 ± 0.11 KeV). The 
effect is also seen by comparison of left and right panel in figure 16.221 

In general the ground-state pseudoscalar mesons are too narrow to 
study a possible variation of their width via the missing mass method 
at any of the up-to-date missing mass spectrometers. For this purpuse 
vector mesons would be much better suited, and in particular mesons 
uj and (f) possessing width of 8.4 MeV and 4.3 MeV, respectively. In 
figure IT^ we demonstrated that the COSY-11 detection setup combined 
with the stochastically cooled proton beam [234] of COSY results in 
a mass resolution of about 0.7 MeV (FWHM). In the case of or u; 
mesons, such a precision would enable to observe effects in the order 
of 10% without necessity of any sophisticated analysis. Unfortunately 
the limited acceptance of COSY-11 setup allows for the reliable (model 
independent) evaluation of the data only in the vicinity of the kinematical 
threshold, which in case of the broad^ mesons makes the derivation of 
their spectral function rather complicated due the strong variation of the 
phase space over a resonance range [6] . A comparable resolution but with 
by far higher acceptance shall be obtainable at the WASA detector. Its 
installation at cooler synchrotron COSY will allow for investigations of 
all ground state pseudoscalar and vector flavour-neutral mesons and we 
consider to perform such a study for the u meson. 

In case the bold hypothesis rised by Maglich [232] were wrong it can be 
easily falsified by the high statistics measurements of the u mesons at the 
upcoming facility WASA@COSY. If confirmed it would be a sensation. 



* For a comprehensive discussion of a notion of the threshold and the definition of the total cross 
section in case of broad resonances the reader is referred to reference [6] . 
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Science is the attempt to discover, by means of obser- 
vation, and reasoning based upon it, first, particular 
facts about the world, and then laws connecting facts 
with one another and (in fortunate cases) making it 
possible to predict future occurrences [235]. 

Bertrand Russell 

Major part of the experimental results, constituting the basis for the 
inferences considered in this work, has been obtained using the COSY-11 
facility. This chapter will be devoted to the description of this detection 
system placing emphasis on its use for the determination of the total 
cross section in the case of the pp ppr]' reaction and derivation of the 
differential cross sections for the reaction pp — > pprj. The description 
will be focused on the data analysis including the methods for the off- 
line monitoring of the beam and target parameters, model independent 
multi-dimensional acceptance corrections, the procedure of kinematical 
fitting and the subtraction of the background, especially in case of the 
derivation of the differential cross sections. 

For details of the functioning of the employed detectors, the descrip- 
tion methods of their calibration, as well as for the elucidation of the 
hardware event selection and the data acquisition the reader is referred 
to references [236,237]. 

6.1 Measurement of the total cross section on the 
example of the pp ppr]' reaction 

/ hold that most observations are more or less indi- 
rect, and that it is doubtful whether the distinction 
between directly observable incidents and whatever is 
only indirectly observable leads us anywhere [238]. 

Karl Raimund Popper 

The COSY-11 facility enables exclusive measurements of the meson pro- 
duction occurring via interaction of nucleons in the range of the kine- 
matical threshold. The coUisions of protons are realized using the cooler 
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synchrotron COSY-Jiilich [62, 239] and the H2 cluster target [240, 241] 
installed in front of one of the regular COSY dipole magnets, as shown 
schematically in figure 16.11 The target being a beam of H2 molecules 
grouped inside clusters of up to 10^ atoms, crosses perpendicularly the 
beam of ~ 2 ■ 10^° protons circulating in the ring. 



Si 




T = 25 K 

Figure 6.1: Schematic view of the part of the COSY-11 detection setup [61]. The cluster 
target [240] is located in front of the accelerator dipole magnet. Positively charged particles 
which leave the scattering chamber through the thin exit foil are detected in two drift cham- 
ber stacks Dl, D2 [242] and in the scintillator hodoscopes SI, S2 [243], and S3 [244,245]. 
Scintillation detector S4 and the position sensitive silicon pad detector Si [246] are used 
in coincidence with the SI counter for the registration of the elastically scattered protons. 
Elastic scattering is used for an absolute normalisation of the cross sections of the investi- 
gated reactions and for monitoring both the geometrical spread of the proton beam and the 
position at which the beam crosses the target [8]. In the left-lower corner a schematic view 
of the interaction region is depicted. The figure has been adapted from reference [237]. 
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The beam of accelerated protons is cooled stochastically during the 
measurement cycle [234]. Longitudinal and vertical cooling enables to 
keep the circulating beam practically without energy losses and without 
a spread of its dimensions even during a 60 minutes cycle, when passing 
more than 10^ times per second through the 10^^ atoms/cm^ thick target. 

If at the intersection point of the cluster beam with the COSY proton 
beam the collision of protons results in the production of a meson, then 
the ejected protons - having smaller momenta than the beam protons - 
are separated from the circulating beam by the magnetic field. Further 
they leave the vacuum chamber through a thin exit foil and are registered 
by the detection system consisting of drift chambers and scintillation 
counters [10,61]. The hardware trigger, based on signals from scintillation 
detectors, was adjusted to register all events with at least two positively 
charged particles. Tracing back trajectories from drift chambers through 
the dipol magnetic field to the target point allowed for the determination 
of the particles momenta. Having momentum and velocity, the latter 
measured using scintillation detectors, it is possible to identify the mass 
of the particle. Figure shows the squared mass of two simultaneously 
detected particles. 




Figure 6.2: (a) Squared masses of two positively charged particles measured in coincidence. 
Pronounced peaks are to be recognized when two protons, proton and pion, two pions, or pion 
and deuteron were registered. Note that the number of events is shown in logarithmic scale. 
(b)(c) Mass spectrum of the unobserved particle or system of particles in the pp ppX 
reaction determined at Q = 5.83 MeV above the rj' production threshold. 



A clear separation is seen into groups of events with two protons, 
two pions, proton and pion, and also deuteron and pion. This spectrum 
enables to select events with two registered protons. The knowledge 
of the momenta of both protons before and after the reaction allows 
to calculate the mass of an unobserved particle or system of particles 
created in the reaction. The four-momentum conservation applied to the 
pp ppX reaction leads to the following expression: 

™2 _ 77'2_ p2 _ / p I 771 _ TPP\'^ _ IP I p — pP^ pP\'^ 

'"'x — X X — {-'-^beami -l-^target -l-^l ^21 l-l^ beam 1 -l^ target 1 -'2 I ) 
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where the used notation is self-explanatory^. Figure 16.2b depicts the 
missing mass spectrum obtained for the pp —>■ ppX reaction at an excess- 
energy of Q = 5.8 MeV above the rj' meson production threshold. Most 
of the entries in this spectrum originate from the multi-pion produc- 
tion [11,237], forming a continuous background to the well distinguish- 
able peaks accounting for the creation of u and rj' mesons, which can be 
seen at mass values of 782 MeV/c^ and 958 MeV/c^, respectively. The 
signal of the pp —>■ ppr]' reaction is better to be seen in the figure 16.2b . 
where the missing mass distribution only in the vicinity of the kinemat- 
ical limit is presented. Figure 16.3b shows the missing mass spectrum for 
the measurement at Q = 7.57 MeV (above kinematical threshold of the 
pp — > ppr]' reaction) together with the multi-pion background (dotted 
line) as combined from the measurements at different excess-energies. 
Subtraction of the background leads to the spectrum with a clear signal 
at the mass of the t]' meson as shown by the solid line in figure IH73b . The 
dashed histogram in this figure corresponds to the Monte-Carlo simula- 
tions where the beam and target conditions were deduced from the mea- 
surements of elastically scattered protons in the manner described in the 
next section. The magnitude of the simulated distribution was fitted to 
the data, but the consistency of the widths is a measure of understanding 
of the detection system and the target-beam conditions. Histograms from 
a measurement at Q = 1.53 MeV shown in figures I6.3b .d demonstrate 
the achieved missing-mass resolution at the COSY-11 detection system, 
when using a stochastically cooled proton beam. The width of the miss- 
ing mass distribution (figure I6.3t i). which is now close to the natural 
width of the rj' meson (F^/ = 0.202 MeV [43]), is again well reproduced 
by the Monte-Carlo simulations. The broadening of the width of the t]' 
signal with increasing excess-energy (compare figures and I6.3l i) is a 
kinematical effect discussed in more detail in reference [10]. The decreas- 
ing of the signal-to-background ratio with growing excess-energy is due 
to the broadening of the rj' peak and the increasing background when 
moving away from the kinematical limit (see e.g. figure l^72b or 16.3b ). 
At the same time, the shape of the background, determined by the con- 
volution of the detector acceptance and the distribution of the two- and 
three-pion production [237], remains unchanged within the studied range 
of the beam momentum. The signal-to-background ratio changes from 
1.8 at Q = 1.53 MeV to 0.17 at Q = 23.64 MeV. The geometrical ac- 
ceptance, being defined by the gap of a dipole yoke and the scintillation 
detector most distant from the target [10,61], decreases from 50 % to 



^ The missing mass technique was first applied by authors of reference [247] . 
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4 % within this excess-energy range. However, in the horizontal plane it 
is still 100 % and thus covers the whole phase space. This issue will be 
considered in details in subsection 16.3.21 
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Figure 6.3: Missing mass distribution with respect to the proton-proton system: (a),(b) 
measurements at Q = 7.57 MeV above the threshold of the pp — » ppri' reaction and (c),(d) at 
Q = 1.53 MeV. Background shown as dotted lines is combined from the measurements at dif- 
ferent energies shifted to the appropriate kinematical limits and normalized to the solid-line 
histogram. Dashed histograms are obtained by means of the Monte-Carlo simulations. 



The absolute value of the excess-energy was determined from the posi- 
tion of the rj' peak in the missing mass spectrum, which should correspond 
to the actual mass of the meson rj' . The systematic error of the excess- 
energy established by this method equals to 0.44 MeV and constitutes 
of 0.14 MeV due to the uncertainty of the rj' meson mass [43] and of 
0.3 MeV due to the inaccuracy of the detection system geometry [248], 
whereby the largest effect originates in the inexactness of relative settings 
of target, dipole, and drift chambers. 
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6.2 Monitoring of beam and target parameters 

Since we can never know anything for sure, it is 
simply not worth searching for certainty; but it is 
well worth searching for truth; and we do this chiefly 
by searching for mistakes, so that we can correct 
them [249]. 

Karl Raimund Popper 

An exact extraction of absolute cross sections from the measured data 
demands a reliable estimation of the acceptance of the detection system. 
This in turn crucially depends on the accuracy of the determination of 
the position and dimensions of both the beam and the target. In this sec- 
tion we will describe a method for estimating the dimensions of the pro- 
ton beam, based on the momentum distribution of elastically scattered 
protons, which can be measured simultaneously with the investigated 
reaction. 

As already pointed out, the production of short-lived uncharged mesons 
{rj, uj, rj', or 0) is investigated at the COSY-11 facility (see figures lUTD 
and 16.4b ) by means of the missing mass technique via the pp — ppX 
reaction, with the four-momenta of protons being fully determined ex- 
perimentally. The accuracy of the extracted missing mass value, which 
decides whether the signal from the given meson is visible over a back- 
ground, depends on the precision of the momentum reconstruction of 
the registered protons, which in turn depends on the detector resolution 
and both the momentum- and the geometrical spread of the accelerator 
proton beam interacting with the internal cluster target beam. The mo- 
mentum reconstruction is performed by tracing back trajectories from 
drift chambers through the dipole magnetic field to the target, which is 
ideally assumed to be an infinitely thin vertical line. In reality, however, 
the reactions take place in that region of finite dimensions where beam 
and target overlap, as depicted in figure lOk . Consequently, assuming in 
the analysis an infinitesimal target implies a smearing out of the momen- 
tum vectors and hence of the resolution of the missing mass signal. The 
part of the COSY-11 detection setup used for the registration of elasti- 
cally scattered protons is shown in figure 16.4b . Trajectories of protons 
scattered in the forward direction are measured by means of two drift 
chambers (Dl and D2) and a scintillator hodoscope (SI), whereas the 
recoil protons are registered in coincidence with the forward ones using 
a silicon pad detector arrangement (Si) and a scintillation detector (S4). 
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Figure 6.4: (a) Schematic description of the relative beam and target setting. Seen from 
above (upper part), and from aside (lower part), ax and oy denote the horizontal and ver- 
tical standard deviation of the assumed Gaussian distribution of the proton beam density, 
respectively. The distance between the centres of the proton and the target beam is de- 
scribed as Ajf. 

(b) Schematic view of the COSY-11 detection setup. Only detectors used for the measure- 
ment of elastically scattered protons are shown. Numbers, at the silicon pad detector (Si), 
and below the scintillator hodoscope (SI), indicate the order of segments. Dl and D2 denote 
drift chambers. The Xsi axis is defined such that the first segment of the SI ends at 80 cm 
and the sixteenth ends at -80 cm. The proton beam, depicted by a shaded line, circulates in 
the ring and crosses each time the H2 cluster target installed in front of one of the bending 
dipole magnets of the COSY accelerator. 



The two-body kinematics gives an unambiguous relation between the 
scattering angles Gi and G2 of the recoiled and forward flying protons. 
Therefore, as seen in figure 16. 5b . events of elastically scattered protons 
can be identified from the correlation line formed between the position in 
the silicon pad detector Si, and the scintillator hodoscope SI, the latter 
measured by the two drift chamber stacks. For those protons which are 
elastically scattered in forward direction and are deflected in the mag- 
netic field of the dipole the momentum vector at the target point can 
be determined. According to two-body kinematics, momentum compo- 
nents parallel and perpendicular to the beam axis form an ellipse, from 
which a section is shown as a solid line in figure lfj.5b . superimposed on 
data selected according to the correlation criterion from figure 16.5b for 
elastically scattered events. 
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Figure 6.5: (a) Identification of elastically scattered protons from the correlation of hits 
in the silicon detector Si and the SI scintillator hodoscope. Note that the number of entries 
per bin is given in a logarithmic scale, ranging from 1 (smallest box) to 19000 (largest box). 

(b) Perpendicular versus parallel (with respect to the beam direction) momentum compo- 
nents of particles registered at a beam momentum of 3.227 GeV/c. The number of entries 
per bin is shown logarithmically. The solid line corresponds to the momentum ellipse ex- 
pected for elastically scattered protons at a beam momentum of 3.227 GeV/c, the dashed 
line refers to a beam momentum of 3.350 GeV/c, and the dotted line shows the momentum 
ellipse obtained for a proton beam inclined by 40 mrad. 

(c) The same data as shown in b) but analyzed with the target point shifted by -0.2 cm 
perpendicularly to the beam direction (along the X-axis in fiaure l^^lb ). The solid line shows 
the momentum ellipse at a beam momentum of 3.227 GeV/c. 

In figure lfj.5b . similarly as in figure If). 5k . it is obvious that data of 
elastically scattered events arise clearly over a certain low level back- 
ground. However, what is important here is that the mean of the elasti- 
cally scattered data is significantly shifted from the expected line, indicat- 
ing that the reconstructed momenta are on average larger than expected. 
This discrepancy cannot be explained by an alternative assumption of ei- 
ther the proton beam momentum or the proton beam angle because of 
the following reason: Trying to reach an acceptable agreement between 
data and expectation the beam momentum must be changed by more 
than 120 MeV/c (dashed line in figure IH3b ). which is 40 times more 
than the conservatively estimated error of the absolute beam momen- 
tum (± 3 MeV/c) [62,248,250]. Similarly, the effect could have been 
corrected by changing the beam angle by 40 mrad (see dotted line in 
figure 16.5b ). which also exceeds the admissible deviation of the beam 
angle (± 1 mrad [62]) by at least a factor of 40. 

However, the observed discrepancy can be explained by shifting the 
assumed reaction point relative to the nominal value by —0.2 cm per- 
pendicular to the beam axis towards the center of the COSY-ring along 
the X-axis defined in figure IH^t- The experimentally extracted momen- 
tum components obtained under this assumption, shown in figure 16.5b . 
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agree with the expectation depicted by the sohd hne, and the data are 
now spread symmetrically around the ellipse. This spread is essentially 
due to the finite extensions of the cluster target and the proton beam 
overlap, which corresponds to about ±0.2 cm as can easily be inferred 
from the value of the target shift required. Other contributions as the 
spread of the beam momentum and still some multiple scattering events 
appear to be negligible and therefore we consider further the infiuence of 
the effective target dimensions and the spread of the COSY beam. 

Assuming both the target beam to be described by a cylindrical pipe, 
with a diameter of 9 mm, homogeneously filled with protons [240,251], 
and the COSY proton beam density distribution to be described by Gaus- 
sian functions with standard deviations ox and oy for the horizontal and 
vertical directions, respectively [62,252], we performed Monte-Carlo cal- 
culations varying the distance between the target and the beam centres 
(Ajf), and the horizontal proton beam extension [px] (see figure IHilt)- 

In order to account for the angular distribution of the elastically scat- 
tered protons, to each generated event an appropriate weight w was as- 
signed according to the differential distributions of the cross sections 
measured by the EDDA collaboration [253]. The generated events were 
analyzed in the same way as the experimental data. A comparison of the 
data from figure 16.5b (data analyzed by using the nominal interaction 
point) with the corresponding simulated histograms allows to determine 
both Ox and Ax. For finding an estimate of the parameters ax and Ax 
we construct the statistic according to the method of least squares: 



where and = w denote the content of the i*'^ bin of the P_ 



^bin 



versus-P|| spectrum determined from experiment (figure IHISb) and sim- 
ulations, respectively. The background events hi in the i*^ bin amount 
to less than one per cent of the data [237] and were estimated by linear 
interpolations between the inner and outer part of a broad distribution 
surrounding the expected ellipse originating from elastically scattered 
protons. The free parameter a allows to adjust the overall scale of the 
fitted Monte-Carlo histograms. Thus, varying the a parameter the Xmm 
for each pair of ax and Ax was established as a minimum of the x^(c^) 
distribution. 



100 



Hadronic interaction of rj and r/' mesons with protons 




Figure 6.6: (a) Xmin a-s a function of ax and Ax- The number of entries is shown in 
a logarithmic scale. (b) Xmin a-s a function of Ax- (c) Xmin as a function of ax. 
The minimum value of Xmm is larger than two- However, it decreases when only a part of 
the 60 minutes long COSY cycle is taken into account (see subsection 16.2. H . This is due 
to the fact that the beam changes during the cycle and its shape, when integrated over the 
whole cycle time, does not suit perfectly to the shape assumed in Monte-Carlo simulations, 
(d) Distribution of the vertical component of the reaction points determined by tracing 
back trajectories from the drift chambers through the dipole magnetic field to the centre of 
the target (in the horizontal plane). "Tails" are due to secondary scattering on the vacuum 
chamber and were parametrized by a polynomial of second order. The solid line shows the 
simultaneous fit of the Gaussian distribution and the polynomial of second order. 

Figure 16.6b shows the logarithm of the obtained xLin ^ function 
of the values of ax and Ax used in the Monte-Carlo calculations: A 
valley identifying a minimum range is clearly recognizable, which gives 
the unique possibility to determine the varied parameters ax and Ax- 
The overall minimum of the Xmm('^x, Ax )-distribution, at ax = 0.2 cm 
and Ax = -0.2 cm, is obviously and better seen in figures and 16.6b . 
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which show the projections of the valley line onto the respective axis. 
The same results were obtained when employing the Poisson likelihood 
derived from the maximum likelihood method [254, 255] 

X' = 2-J2 ["^^ + h-N: + N: (6.2) 

The vertical beam extension of cry = 0.51 cm was established di- 
rectly from the distribution of the vertical component of the particle 
trajectories at the centre of the target. This is possible, since for the 
momentum reconstruction only the origin of the track in the horizontal 
plane is used, while the vertical component remains a free parameter. 
As shown in figure 16. (jd . the reconstructed distribution of the vertical 
component of the reaction points indeed can well be described by a 
Gaussian distribution with the a = 0.53 cm. The width of this dis- 
tribution is primarily due to the vertical spread of the proton beam. The 
spread caused by multiple scattering and drift chambers resolution was 
determined to be about 0.13 cm (standard deviation) [237]. Therefore, 

ay = - 0.13^ ^ 0.51 cm. 

The effect of a possible drift chamber misalignment, i. e. an inexact- 
ness of the angle (3 in figure 16.4b . was estimated to cause a shift in the 
momentum plane corresponding to a value of 0.15 cm for Ax- This gives 
a rather large systematic bias to the estimation of the absolute value of 
Ax, but it does not influence the implications concerning the parame- 
ter ax, and still allows for the determination of relative beam shifts Ax 
during the measurement cycle. 

The absolute value of Ax can be found if one additionally takes into 
account the width of the mass spectrum of the studied meson. Clearly, 
using in the analysis a wrong value of Ax broadens a missing mass dis- 
tribution. Therefore, the real value of Ax can be determined as the one 
at which the width of the mass spectrum of the investigated meson is the 
least. 

6.2.1 Influence of the stochastic cooling on the proton beam 
quality 

The cooling of a single particle circulating in a 
ring is particularly simple [256]. 

Simon van der Meer 

Meson production cross sections in the proton-proton collisions increase 
rapidly with the beam momentum, and change by more than one or- 
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der of magnitude when the excess energy grows from 1 MeV to 10 
MeV [2,3,10]. Due to the rapid change of the total cross section, stud- 
ies of the close-to-threshold meson production require a precise beam 
with both small momentum and geometrical spread, and an accurate 
determination of the momenta of the reaction products. The use of a 
stochastically cooled beam [62,234,256,257] and of internal cluster tar- 
get facilities allows to realize these conditions. Very low target densities 
10 atoms cm~^) [240] minimize changes in the ejectiles' momentum 
vectors due to secondary scattering in the target, and simultaneously al- 
low the beam to circulate through the target without significant intensity 
losses and changes of the momentum. An improvement of the experimen- 
tal conditions caused by the stochastic cooling can be deduced already 
from a time dependence of the trigger frequency presented in figure 16.71 




Figure 6.7: Rate of coincidences between the detectors SI and S4 shown in figure IS^lb . 
In the middle of the second spill the cooling system was switched on. 



This figure presents the rate of the coincidence between the scintil- 
lator detectors S4 and SI. The observed event rate is proportional to 
the number of reactions occuring in the region of the beam and target 
overlap. A steep decrease in the beginning of the first and second shown 
cycles is caused by the outward movement of the uncooled beam from 
the target center. Switching on the longitudinal and horizontal stochas- 
tic cooling, during the second cycle, clearly increases the luminosity and 
in the subsequent spills prevents the beam from spreading and shifting 
out of the target location. Consequently, a constant counting rate is ob- 
tained during the cycles with a cooled beam. In the following the beam 
monitoring technique introduced in the previous section will be applied 
to inspect in more details an influence of the stochastic cooling on the 
proton beam quality. 
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The horizontal stochastic coohng [62] is intended to squeeze the proton 
beam in the horizontal direction until the beam reaches the equilibrium 
between the cooling and the heating due to the target. 
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Figure 6.8: The vertical a) and horizontal b) beam width (one standard deviation) deter- 
mined for each five minutes partition of the COSY cycle. 

Please note that the beam width ax determined for each five minutes period is smaller than 
the (Tx averaged over the whole cycle (compare figure 15^ 1. This is due to the beam shift 
relative to the target during the experimental cycle as shown in the panel c). 
c) Relative settings of the COSY proton beam and the target centre versus the time of the 
measurement cycle. 

The vertical error bars in pictures b) and c) denote the size of the step used in the 
Monte-Carlo simulations (± 0.025 cm ; the bin width of figures l(T!Hb and 16.6b ). 
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Figure 6.9: Perpendicular versus parallel momentum components with respect to the beam 
direction of particles registered at a beam momentum of 3.227 GeV/c, as measured during 
the first (a), and the last minute (b) of the 60 minutes long COSY cycle. The data at 
the first minute were analyzed with Ajc = -0.15 cm, and the data of the 60's minute with 
Ax ~ -0.25 cm, see figure l6^ . The number of entries per bin is shown logarithmically. The 
solid line corresponds to the momentum ellipse expected for protons scattered elastically at 
a beam momentum of 3.227 GeV/c. 
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In case of an uncooled beam its size would increase during the cycle as 
can be seen in figure 16.8b . which shows the spreading of the beam in the 
vertical plane during the 60 minutes cycle. This was expected since dur- 
ing the discussed experiment the stochastic cooling in the vertical plane 
was not used. The influence of the applied cooling in the horizontal plane 
is clearly visible in figure IHiHb- During the first five minutes of the cycle 
the horizontal size of the beam of about 2 ■ 10^° protons was reduced by a 
factor of 2 reaching the equilibrium conditions, and remaining constant 
for the rest of the COSY cycle. Figure 16.91 depicts the accuracy of the 
momentum reconstruction, reflected in the spreading of the data, which 
is mainly due to the finite horizontal size of the beam and target overlap. 
The left and right panel correspond to the first and the last minute of 
the measurement cycle, respectively. The data were analyzed correcting 
for the relative target and beam shifts A^. The movement of the beam 
relative to the target, during the cycle, is quantified in figure lUTSh . The 
shift of the beam denotes also changes of the average beam momentum, 
due to the nonzero dispersion at the target place. After some minutes 
the beam remains unchanged, this can be understood as reaching the 
equlibrium between the energy losses when crossing 1.6 ■ 10^ times per 
second through the H2 cluster target, and the power of the longitudinal 
stochastic cooling which cannot only diminish the spread of the beam 
momentum but also shifts it whole. 

6.3 Detailed data analysis on the example of the 
pp ppj] reaction 

Precision and exactness are not intellectual values 
in themselves, and we should never try to he more 
precise or exact than is demanded by the problem in 
hand [249]. 

Karl Raimund Popper 

Significant part of the differential distributions of the cross sections used 
in this work was derived from a high statistics measurement of the 
pp —>■ pprj reaction at a nominal beam momentum of 2.027 GeV/c. The 
experiment performed at the COSY-11 facility was based on the four- 
momentum registration of both outgoing protons, whereas the rj meson 
was identified via the missing mass technique. The method allows for 
the extraction of the kinematically complete information of the produced 
pprj system, provided that the reaction was identified. This is not feasible 
on the event-by-event basis due to the unavoidable physical background 
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originating from multi-pion production. However, with a large number 
of identified events the background subtraction can be performed sepa- 
rately for each interesting phase space interval. In practice, the size of 
the studied phase space partitions must be optimized between the statis- 
tical significance of the signal-to-background ratio and the experimental 
resolution. 
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Figure 6.10: (a) The dashed curve denotes the proton beam momentum distribution 
integrated over the whole measurement period. On the horizontal scale, the value of zero 
is set at a nominal beam momentum equal to 2.027 GeV/c. The solid line shows the beam 
momentum distribution after the correction for the mean value which was determined in 
10 second intervals, (b) Missing mass distribution for the pp ppX reaction determined 
by means of the COSY-11 detection system at a beam momentum of 2.027 GeV/c. The 
solid histogram presents the data corrected for effects of the time dependent relative shifts 
between the beam and the target using the method described in section [6.21 The dashed 
histogram shows the result before the correction. 



The accuracy of the missing mass reconstruction depends on two pa- 
rameters: the spread of the beam momentum as well as the precision 
of the measured momenta of the outgoing protons. The latter is pre- 
dominantly due to the geometrical spread of the beam. Since for the 
reconstruction of the momenta we assume that the reaction takes place 
in the middle of the target we cannot correct on an event-by-event ba- 
sis for the momentary spread of the beam. However, we can rectify 
the smearing due to the shifts of the centre of the beam relative to the 
target as well as the average changes of the absolute beam momentum 
during the experiment. Therefore, after the selection of events with two 
registered protons, as a first step of the more refined analysis the data 
were corrected for the mean beam momentum changes (see fig. l6.1Ub ) de- 
termined from the measured Schottky frequency spectra and the known 
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beam optics. In a next step, from the distributions of the elastically scat- 
tered protons, the Schottky frequency spectrum, and the missing mass 
distribution of the pp — ppX reaction, we have estimated that the spread 
of the beam momentum, and the spread of the reaction points in hori- 
zontal and vertical direction amount to cr^pheam) = 0.63 ± 0.03 MeV/c, 
a{x) = 0.22 ± 0.02 cm, and a{y) = 0.38 ± 0.04 cm, respectively. 
Details of this procedure can be found in section 16.21 

Figure 6.11: Missing mass spectrum for 
the pp ppX reaction determined in 
the experiment at a beam momentum 
of 2.0259 GeV/c. The mass resolution 
amounts to 1 MeY /c^{a). The superim- 
posed histograms present the simulation 
for 1.5 • 10* events of the pp — > pprj re- 
action, and 10^" events for the reactions 
pp pp2tt, pp pp3n, and pp ppAn. 
The simulated histograms were fitted to 
the data varying only the magnitude. The 
fit resulted in 24009 ± 210 events for the 
0.54 0.55 0.56 production of the 77 meson. 

missing mass [GeV/c^] 

Further on, a comparison of the experimentally determined momen- 
tum spectra of the elastically scattered protons with the distributions 
simulated with different beam and target conditions allows us to estab- 
lish the position at which the centre of the beam crosses the target with an 
accuracy of 0.25 mm. Accounting for the movement of the beam relative 
to the target we improved the missing mass resolution as demonstrated 
in figure IFTTUb . 

After this betterment, the peak originating from the pp pprj re- 
action became more symmetric and the signal-to-background ratio in- 
creased significantly. The background was simulated taking into account 
pp ppX reactions with X = 27c, Stt, and An. Since we consider here 
only the very edge of the phase space distribution where the protons are 
produced predominantly in the S-wave the shape of the background can 
be reproduced assuming that the homogeneous phase space distribution 
is modified only by the interaction between protons. Indeed, as can be ob- 
served in figure lUTTD the simulation describes the data very well. The cal- 
culated spectrum is hardly distinguishable from the experimental points. 
The position of the peak on the missing mass spectrum and the known 
mass of the rj meson [43] enabled to determine the actual absolute beam 
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momentum to be Pbeam = 2.0259 GeV/c ± 0.0013 GeV/c, which agrees 
within error hmits with the nominal value of = 2.027 GeV/c. 

The real beam momentum corresponds to an excess energy of the pprj 
system equal toQ = 15.5±0.4 MeV. 



6.3.1 Covariance matrix and kinematical fitting 

Beobachtungen, welche sich auf Grossenbesti- 
mmungen aus der Sinnenwelt beziehen, werden 
immer, so sorgfdltig man auch verfahren mag, 
grosseren oder kleineren Fehlern unterworfen 
bleiben [258]. 

Carl Friedrich Gauss 

As already mentioned in the previous section, at the COSY-11 facility 
the identification of the pp — *• pprj reaction is based on the measurement 
of the momentum vectors of the outgoing protons and the utilisation of 
the missing mass technique. Inaccuracy of the momentum determination 
manifests itself in the population of kinematically forbidden regions of the 
phase space, preventing a precise comparison of the theoretically derived 
and experimentally determined differential cross sections. 



Figure 6.12: Dalitz plot distribution of 
the pp — + ppr) reaction simulated at 
Q = 15.5 MeV. The number of entries is 
shown in a logarithmic scale. The solid 
line gives the kinematically allowed area. 
The result was obtained taking into ac- 
count the experimental conditions as de- 
scribed in the text. 
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Figure 16.121 visualizes this effect and clearly demonstrates that the 
data scatter significantly outside the kinematically allowed region (solid 
line), in spite of the fact that the precision of the fractional momentum 
determination in the laboratory system {<7{piab)/piab ~ 7 ■ 10^^) is quite 
high. Therefore, when seeking for small effects like for example the influ- 
ence of the proton-?7 interaction on the population density of the phase 
space, one needs either to fold the theoretical calculations with the exper- 
imental resolution, or to perform a kinematical fitting of the data. Both 
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procedures require the knowledge of the covariance matrix, and thus its 
determination constitutes a necessary step in the differential analysis and 
interpretation of the data. 

In order to derive the covariance matrix we need to recognize and 
quantify all possible sources of errors in the reconstruction of the two 
proton momenta pi and p2. The four dominant effects are: i) finite dis- 
tributions of the beam momentum and of the reaction points, ii) multiple 
scattering in the dipole chamber exit foil, air, and detectors, iii) finite 
resolution of the position determination of the drift chambers, and iv) 
a possible inadequate assignment of hits to the particle tracks in drift 
chambers in the case of very close tracks. Some of these, hke the mul- 
tiple scattering, depend on the outgoing protons' momenta, others, like 
the beam momentum distribution, depend on the specific run conditions 
and therefore must be determined for each run separately. 

In order to estimate the variances and covariances for all possible com- 
binations of the momentum components of two registered protons we 
have generated 1.5 • 10^ pp — * ppt] events and simulated the response of 
the COSY- 11 detection setup taking into account the above listed factors 
and the known resolutions of the detector components. Next, we anal- 
ysed the signals by means of the same reconstruction procedure as used 
in case of the experimental data. Covariances between the i*'* and the 
jth components of the event vector (P = [pix,Piy,Piz,P2x,P2y,P2z]) were 
established as the average of the product of the deviations between the 
reconstruced and generated values. The explicit formula for the sample 
of N reconstructed events reads: 



for the i*'* component of the vector P describing the A;*'* event. 

Because of the inherent symmetries of the covariance matrix {c(nj{i, j) = 
cov{j,i)) and the indistinguishability of the registered protons^ there are 
only 12 independent values which determine the 6x6 error matrix V 
unambiguously. 

Since inaccuracies of the momentum determination depend on the 
particle momentum itself (eg. multiple scattering) and on the relative 
momentum between protons (eg. trajectories reconstruction from sig- 
nals in drift chambers), we have determined the covariance matrix as a 



^ The symmetry of all observables under the exchange of the two protons (pi <-» P2) implies that 
cov(i,j) = cov(i±3,j±3), where the '+' has to be taken for i,j = 1,2,3 and the '-' for i,j = 4,5,6. Thus 
for example cov(2,4) = cov(5,l). 




k=l 



where /^^^^ and f"^ 



denote the generated and reconstructed values 
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function of the absolute momentum of both protons: cov{i,j, \pi\, \p2\)- 
As an example we present the covariance matrix for the mean values of 
\pi\ and IP2I in units of MeV^/c^, as established in the laboratory sys- 
tem with the z-coordinate parallel to the beam axis and y-coordinate 
corresponding to the vertical direction. 



Plx Ply 

5.6 0.0 

- 7.1 



Plz 

-13.7 
0.1 
37.0 



P22 
1.7 



P2y P22 

0.1 -3.0 
-0.2 -0.2 
5.4 



Plx 
Ply 

V = - - 37.0 - - 5.4 pi, (6.3) 

P2x 
P2y 
P2z 

Since the measurements have been performed close to the kinematical 
threshold the ejectile momentum component parallel to the beam is by 
far the largest one and its variance {var{pz) = 37 MeV'^/c^) determines 
in first order the error of the momentum measurement. The second 
largest contribution stems from an anti-correlation between the x— and 
z— momentum components {cov{px,Pz) = — 13.7 MeV"^ /c^), which is 
due to the bending of the proton trajectory - mainly in the horizontal 
direction - inside the COSY-11 dipole magnet (see fig. l6.1|) . There is also 
a significant correlation between the z components of different protons 
which is due to the smearing of the reaction points, namely, if in the 
analysis the assumed reaction point differs from the actual one, a mistake 
made in the reconstruction affects both protons similarly. 
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Figure 6.13: Variances of the protons's momentum components a^{px), a'^{pz), and a'^{py) 
shown as a function of the absolute values of the measured momenta. 



Figure inHSl depicts the variation of var{px), var{py), and var^p^) over 
the momentum plane {\pi\, \p2\)- Taking into account components of 
the covariance matrices V(|pi|, \p2\) and the distribution of the proton 
momenta for the pp pprj reaction at Q = 15.5 MeV results in an average 
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error for the measurement of the proton momentum of about 6 MeV/c. 
This can be also deduced from the distribution of the difference between 
the generated and reconstructed absolute momenta of the protons. The 
corresponding spectrum is plotted as a dashed line in figure 16.141 



Figure 6.14: Spectrum of differences 
between generated and - after simulation 
of the detector response - reconstructed 
absolute momenta of protons, as deter- 
mined before (dashed line) and after the 
kinematical fit (solid line). 
The picture shows results obtained taking 
into account the experimental conditions 
as described in the text. 
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In the experiment we have measured 6 variables and once we assumed 
that the event corresponds to the pp — > ppr] reaction only 5 of them are 
independent. Thus we have varied the values of the event components 
demanding that the missing mass is equal to the mass of the t] meson and 
we have chosen that vector which was the closest to the experimental one 
in the sense of the Mahalanobis distance. The inverse of the covariance 
matrix was used as a metric for the distance calculation. The kinematical 
fit improves the resolution by a factor of about 1.5 as can be concluded 
from comparing the dashed and solid lines in figure 16.141 The finally 
resulting error of the momentum determination amounts to 4 MeV/c. 

6.3.2 Multidimensional acceptance corrections 

Unfehlbar hdtte ich umkommen miissen, wenn ich 
mich nicht mit der Starke meiner Arme an meinem 
eigenen Haarzopf samt meinem Pferd, das ich fest 
zwischen meine Knie schlop, wieder herausgezogen 
hdtte [259]. 

Baron Miinchhausen 

At the excess energy of Q = 15.5 MeV the COSY-11 detection system 
does not cover the full An solid angle in the centre-of-mass system of the 
pp — >• ppr] reaction. Therefore, the detailed study of differential cross sec- 
tions requires corrections for the acceptance. Generally, the acceptance 
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should be expressed as a function of the full set of mutually orthogonal 
variables which describe the studied reaction unambiguously. As intro- 
duced in section 13 .41 to define the relative movement of the particles in 
the reaction plane we have chosen two squares of the invariant masses: 
Spp and Sp^, and to define the orientation of this plane in the center- 
of-mass frame we have taken the three Euler angles: The first two are 
simply the polar 0* and azimuthal 6* angles of the momentum of the r) 
meson and the third angle ip describes the rotation of the reaction plane 
around the axis defined by the momentum vector of the rj meson. 

In the data evaluation we considerably benefit from the basic geomet- 
rical symmetries satisfied by the pp — > ppr] reaction. Due to the axial 
symmetry of the initial channel of the two unpolarized colliding protons 
the event distribution over 0* must be isotropic. Thus, we can safely inte- 
grate over 0*, ignoring that variable in the analysis. Furthermore, taking 
advantage of the symmetry due to the two identical particles in the initial 
channel, without losing the generality, we can express the acceptance as 
a function of Spp,Spr^,\cos{9*)\, and ip- To facilitate the calculations we 
have divided the range of \cos{9*)\ and ip into 10 bins and both Spp and 
Spr) into 40 bins each. In the case of Spp and the choice was made 
such that the width of the interval corresponds to the standard devia- 
tion of the experimental accuracy. For |cos(6'*)| and ip we have taken 
only ten partitions since from the previous experiments we expect only 
a small variation of the cross section over these variables [83,95,96]. In 
this representation, however, the COSY-11 detection system covers only 
50% of the phase space for the pp —>■ ppr] reaction at Q = 15.5 MeV. 
To proceed with the analysis we assumed that the distribution over the 
angle ip is isotropic as it was for example experimentally determined for 
the pp ppuj, pp — i> ppp, or pp ppcp reactions [83,89]. Please note 
that this is the only assumption of the reaction dynamics performed in 
the present evaluation. The validity of this supposition in the case of the 
pp pprj reaction will be discussed later. 

In the calculations we exploit also the symmetry of the cross sections 
under the exchange of the two identical particles in the final state which 
reads: a{sp^ri,ip) = ^"(spjr;, ^ + tt). The resultant acceptance is shown 
in figure 16.151 One recognizes that only a small part (3%) of the phase 
space is not covered by the detection system. In further calculations 
these holes were corrected according to the assumption of a homogeneous 
phase space distribution. Additionally it was checked that the corrections 
under other suppositions e.g. regarding also the proton-proton FSI lead 
to negligible differences. 



112 



Hadronic interaction of 77 and 77' mesons with protons 




<D 10 
O 

c 
CO 

■5.10 

CD 
O 

^ 10 



CD 
O 

^ 10 
Q. 

8 10 

o 

CO 

10 



2000 4000 
phase space bin 



k 



-4 



0. < 



cose* 



<0.1 




10 r 



200 400 

phase space bin 
TOT 




5200 5400 
phase space bin 

Figure 6.15: Acceptance of the COSY-11 detection system for the pp pprj reaction at 
an excess energy of Q = 15.5 MeV presented as a function of Spp,Spr,, and \cos{6*)\. The 
numbers were assigned to the bins in the three dimensional space Spp — Spr, — \cos{0*)\ by 
first incrementing the index of Spp next of Sp^ and in the end that of \cos{9^)\. Partitioning 
of |cos(0*)| into ten bins is easily recognizable. The two lower pictures show the acceptance 
for the first and the last bin of |cos(^*)|. 
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Figure 6.16: (a) Distribution of tlie polar angle of the emission of the rj meson in the 
centre-of-mass system. Experimental data were corrected for the acceptance in the three 
dimensional space (spp, Spn, |cos(6'^)|). Full circles show the result with the assumption that 
the distribution of the tp angle is isotropic, and the open circles are extracted under the 
assumption that ^ is as derived from the data (see text). Both results have been normalized 
to each other in magnitude. 

(b) COSY- 11 detection acceptance as a function of Spp and Sp,,, calculated under the 
assumption that the differential cross sections j^^^^e*) ^^'^ ^ isotropic. 
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Figure 6.17: (a) Dalitz plot distribution simulated for the pp pprj reaction at 
Q — 15.5 MeV. In the calculations the interaction between protons was taken into account, 
(b) Dalitz plot distribution reconstructed from the response of the COSY- 11 detectors simu- 
lated for events from figure (a) taking into account the smearing of beam and target, multiple 
scattering in the materials, and the detectors' resolution. The evaluation included momen- 
tum reconstruction, kinematical fitting, and the acceptance correction exactly in the same 
way as performed for the experimental data. The lines surrounding the Dalitz plots depict 
the kinematical limits. 



Full points in figure 16.16b present the distribution of the polar angle 
of the Tj meson as derived from the data after the acceptance correction. 
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Within the statistical accuracy it is isotropic. Taking into account this 
angular distribution of the cross section we can calculate the acceptance 
as a function of Spp and Sprj only. This is shown in figure lUTT^ . where one 
sees that now the full phase space is covered. This allows us to deter- 
mine the distributions of Spp and Sp^. The correctness of the performed 
procedures for the simulation of the detectors response, the event recon- 
struction programs, the kinematical fitting, and acceptance correction 
can be confirmed by comparing the distribution generated (figure 16.17b ) 
with the ones which underwent the complete analysis chain described in 
this section (figure IHUZb)- 
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Figure 6.18: Distribution of the cross 
section as a function of the angle as 
determined in the first iteration. The su- 
perimposed histogram corresponds to the 
fit of the function ^ = a + b ■ \sin{ip)\. 
The range of the tp angle is shown 
from to TT only, since in the analy- 
sis we take advantage of the symmetry 



Knowing the distribution of the polar angle of the rj meson 6* and 
those for the invariant masses Spp and we can check whether the 
assumption of the isotropy of the cross section distribution versus the 
third Euler's angle ip is corroborated by the data. For that purpose 
we calculated the acceptance as a function of ip and assuming the 
shape of the differential cross sections of and , as determined 

^ dspp dcos(f)*) 

experimentally. The obtained ^ distribution is shown in figure lOSl and 
is not isotropic as assumed at the beginning. A fit of the function of the 
form ^ = ci + b ■ \sin{il)) \ gives the value of b = 0.079 ± 0.014 jjb/sr, 
indeed significantly different from the isotropic solution. This deviation 
cannot be assigned to any unknown behaviour of the background since 
the obtained distribution can be regarded as background free. This is 
because the number of pp —>■ ppr] events was elaborated for each invariant 
mass interval separately. The exemplarily missing mass spectra for the 
first, fourth, seventh, and tenth interval of ip values, corrected for the 
acceptance, are presented in figure 16.191 From this figure one can infer 
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that the shape of the background is well reproduced not only for the 
overall missing mass spectrum as shown previously in figure 16.111 but 
also locally in each region of the phase space. Since the experimental 
data are quite well described by the simulations one can rather exclude 
the possibility of a significant systematic error which could cause the 
observed anisotropy of the differential cross section 
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Figure 6.19: Missing mass distributions for the first, fourtli, seventfi, and tenth bin of 
i/) with the superimposed lines from the simulation of pp — » pprj and the multi-pion back- 
ground pp pp(mn) reactions. Amplitudes of simulated distributions were fitted to the 
experimental points. 



The evaluated distribution is however in disagreement with the work- 
ing assumption that (y{ip) is isotropic. Therefore we performed a full 
acceptance correction procedure from the very beginning assuming that 
the distribution of ^ is as determined from the data. After repeating 
the procedure three times we observed that the input and resultant dis- 
tributions are in good agreement. The result after the third iteration is 
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shown in figure I6.2UI by the full circles. It is only slightly different from 
the one obtained after the first iteration as shown in figure 16.181 
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Figure 6.20: Distribution of the cross section as a function of the angle t/i. Full circles stand 
for the final results of the ^ obtained after three iterations. The superimposed histogram 
(solid line) corresponds to the fit of the function — a + b ■ \sin{ip)\ which resulted 
in a = 0.186 ± 0.004 fib/sr and b = 0.110 ± 0.014 /ife/sr. The dashed line shows the 
entry distribution used for the second series of iterations as described in the text. Open circles 
represent the data from the left upper corner of the Dalitz plot (see for example figure lSTTSb l. 
At that region of the Dalitz plot due to the non-zero four dimensional acceptance over {spp, 
Spri, \cx>s{6^)\, ij)) bins the spectrum (open circles) was corrected without a necessity of any 
assumptions concerning the reaction cross section. 

To raise the confidence of the convergence of the performed iteration 
we accomplished the full procedure once more, but now assuming that 
the distribution of ^ is much more anisotropic than determined from 
the data. As an entry distribution we took the dashed line shown in 
figure I6.2UI Again after two iterations we have got the same result as 
before. To corroborate this observation we have evaluated the distribu- 
tion over the if) angle (see figure inSOI) from the phase space region which 
has no holes in the acceptance expressed as a four- dimensional function 
of the variables Spp, Sp^, |cos(6'*)|, and ip, this is for the values of Spp and 
Spri corresponding to the upper left corner of figure IU?TBb . Again the ob- 
tained distribution presented as open circles in figure inHDl is anisotropic, 
and moreover agrees with the spectrum determined from all events. It is 
important to note that the shape of the Spp, Sp^, and cos(^*) distributions 
keeps unchanged during the whole iteration procedure. Figure lU". 2 II shows 
the distributions of the square of the proton-proton and proton-i] invari- 
ant masses. The spectra after the second and third iterations are shown. 
One recognizes that the form of the spectra remains unaltered. The same 
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conclusion can be drawn for the 
figure 16.16b . 
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distribution as demonstrated in 
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Figure 6.21: Distributions of the invariant masses Spp and Sp^ determined for the two 
different assumptions about the cross section dependence of the t/j angle. 

From that comparison one can conclude that the shapes of the deter- 
mined distributions are - within the statistical accuracy - independent 
of the shape of the ^ cross section and can be treated as derived in a 

completely model independent manner. Similarly, as in the case of the ^ 
distribution, the differential cross sections in all other variables reported 
in this article are not deteriorated by the background. This is because 
the number of pp — > pprj events was determined for each investigated 
phase space interval separately. 
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Figure 6.22: Missing mass distributions determined for the invariant mass bins as depicted 
inside the figures. The spectra were corrected for the acceptance. The histograms show the 
simulations for the multi-pion background pp — > pp{mn) and the pp pprj reactions fitted 
to the data with the amplitudes as the only free parameters. 
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As an example the missing mass spectra for two bins of the proton-?^ 
invariant mass are presented in figure I6.221 As aheady noticed for ^ the 
shape of the background is well reproduced locally in each region of the 
phase space. 

6.3.3 Angular distributions 

Before we declare our consent we must carefully examine 

the shape of the architecture the rhythm of the drums and pipes 

... [260]. Zbigniew Herbert 

In previous subsections we derived the distributions of the cross sec- 
tion in the invariant masses Spp and (figure I6.21|) . and in angle "0 
(figure IU.20j) . For the sake of completeness we will also present the dis- 
tribution in the 005(6**), cos{6*^), and ipN- This will allow us to have 
insight into distributions in all nontrivial variables describing the ppr^ 
system in the case when the polarization of nucleons is ignored. Here, 
as introduced in section 13.41 we consider two different sets of orthogonal 
variables, namely {spp,Spr^,\cos{6^)\,ip) and (spp,Sp^,|cos(6'^)|,?/'Ar). If pos- 
sible the data will be compared to the results of measurements performed 
at the non-magnetic spectrometer COSY-TOF [95]. 

In the previous subsection we have shown that even close-to-threshold 
at Q = 15.5 MeV the production of the pprj system is not fully isotropic. 
In particular, we found that the cosine of angle if) is not uniformly popu- 
lated. The anisotropy of the cross section reflects itself in an anisotropy of 
the orientation of the emission plane, and the latter has a simple physical 
interpretation. 

The determined cross section distribution in function of the polar 
angle 9*^^ of the vector normal to that plane is shown in figure 16.23b . 
The distribution is not isotropic, which is particularly visible for the low 
values of 1005(6*^)1 burdened with small errors. As depicted in figure IX^ 
the I cos{6*^^\ — denotes such configurations of the ejectiles momenta 
in which the emission plane comprises the beam axis. In that case the 
acceptance of the COSY-11 detection system is much larger than for the 
configuration where the emission plane is perpendicular to the beam. 
Due to this reason the error bars in figure 16.23b increase with growth of 
1005(6*^)1. It is worth to stress that the tendency of the pprj system to be 
produced preferentially if the emission plane is perpendicular to the beam 
is in line with the preliminary analysis of the experiment performed by 
the TOF collaboration [261]. Elucidation of that non-trivial behaviour 
can reveal interesting features of the dynamics of the production process. 
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Figure 6.23: (a) Differential cross section as a function of the polar angle of the vector 
normal to the emission plane. 

(b) Differential cross section of the pp — > ppr^ reaction as a function of the rj meson cen- 
tre-of-mass polar angle. Full circles depict experimental results for the pp — > pprj reaction 
measured at Q = 15.5 MeV by the COSY-11 collaboration (this article) and the open circles 
were determined by the TOF collaboration at Q = 15 MeV [95]. The TOF points were 
normalized in amplitude to our result, since for that data the absolute scale is not evaluated. 
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Figure 6.24: Differential cross sec- 
tion in -i/jjv for the pp ppr] reaction 
measured at Q = 15.5 MeV. The vari- 
able i/iiv is defined in section mi 
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The distribution of the polar angle of the rj meson emission in the 
centre-of-mass system is shown in figure 16.23b . Clearly, the COSY-11 
data agree very well with the angular dependence determined by the 
TOF collaboration. As already mentioned, the two identical particles 
in the initial state of the pp ppr] reaction imply that the angular 
distribution of either ejectile must be symmetric around 90 degree in the 
centre-of-mass frame. Here we use that reaction characteristic, yet in 
figure 14.3b (see also reference [18] ) we presented the differential cross 
section of the i] meson centre-of-mass polar angle for the full range of 
cos{9*) and found that this is completely symmetric around cos{9*) = 0. 
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These observations can be regarded as a check of the correctness of the 
acceptance calculation. The cross section distribution in the angle ipN 
- defining the orientation of the ppr] system within the emission plane - 
is shown in figure 16.241 

Up to now we presented invariant mass spectra of the proton-proton 
and proton-r] systems and angular distribution for two sets of non-trivial 
angles which describe the orientation of ejectiles within the emission 
plane and the alignment of the plane itself, namely {cos{9*), ip) and 
{cos{6y), iPn)- However, since one of the important issues discussed in 
this work is the contribution from higher partial waves we evaluated also 
an angular distribution of the relative momentum of two protons seen 
from the proton-proton centre-of-mass subsystem (see figure IHIlSt') • The 
distribution of that angle should deliver information about the partial 
waves of the proton-proton system in the exit channel (see section IX^ . 
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Figure 6.25: (a) 

at Q = 15.5 MeV. 

(b) Distribution of the centre-of-mass polar angle of the relative protons momentum with 
respect to the beam direction determined for the pp —> ppr] reaction at Q = 15.5 MeV. The 
COSY-11 result (closed circles) is compared to the data points determined at Q = 15 MeV 
by the TOF collaboration (open circles) [95]. 

Figure l(i.25b presents the differential cross section as a function of 
angle 9^^ of the relative proton momentum seen from the overall centre- 
of-mass frame, as this is often considered in the theoretical works. This 
figure compares COSY-11 results to the angular distribution extracted 
by the TOF collaboration. Both experiments agree very well within 
the statistical accuracy, and both indicate a slight decrease of the cross 
section with increasing |cos(6'*p)|. 
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6.4 Usage of the spectator model for the study of 
the T] and r]' mesons via the proton-neutron 
interaction 

As regards any subject we propose to investigate, we 
must inquire not what other people have thought, 
or what we ourselves conjecture, but what we can 
clearly and manifestly perceive by intuition or deduce 
with certains. For there is no other way of acquiring 
knowledge [207]. 

Rene Descartes 

In order to measure the pn puMeson reactions by means of the proton 
beam it is necessary to use a nuclear target, since a pure neutron target 
does not exist. Naturally, least complications in the interpretation of the 
experiment will be encountered when using the simplest nuclei. There- 
fore, here deuterons will be considered as a source of neutrons, and for 
the evaluation of the data an impulse approximation will be exploited. 
The main conjecture of this approach is that the bombarding proton in- 
teracts exclusively with one nucleon in the target nucleus and that the 
other nucleons affect the reaction by providing a momentum distribution 
to the struck constituent only (fig. I6.26j) . This assumption is justified if 
the kinetic energy of a projectile is large compared to the binding en- 
ergy of the hit nucleus. In fact, as noticed by Slobodrian [262], also the 
scattering of protons on a hydrogen target, where the protons are bound 
by molecular forces, may serve as an extreme example of the quasi-free 
reaction. In that case, although the hydrogen atoms rotate or vibrate 
in the molecule, their velocities and binding forces are totally negligible 
with respect to the velocity and nuclear forces operating the scattering 
of the relativistic proton [262]. 



d 

Figure 6.26: Spectator model for a particle production reaction via pd pspX. 

The deuteron is also relatively weakly bound with a binding energy of 
^ 2.2 MeV, which is by far smaller - more than two orders of magni- 
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tude in the case of pion and already more than three orders of magnitude 
in the case of meson - than the kinetic energy of the bombarding pro- 
tons needed for the creation of mesons in the proton-neutron interaction. 
However, it has to be considered that even the low binding energy of 
2.2 MeV results in large Fermi momenta of the nucleons which can 
not be neglected. The momentum and kinetic energy distributions of 
the nucleons in the deuteron are shown in figure 16.271 In the considered 
approximation the internucleon force manifests itself only as the Fermi 
motion of the nucleons and hence the struck neutron is treated as a free 
particle in the sense that the matrix element for quasi-free meson produc- 
tion off a bound neutron is identical to that for the free pn — pn Meson 
reaction. The reaction may be symbolically presented as: 

p f ''T' \ — ^ P Meson , ^ 

\P ) Psp 

where Psp denotes the proton from the deuteron regarded as a specta- 
tor which does not interact with the bombarding particle, but rather 
escapes untouched and hits the detectors carrying the Fermi momentum 
possessed at the moment of the collision. 




Figure 6.27: (a) Momentum and (b) kinetic energy distribution of the nucleons in the 
deuteron, generated according to an analytical parametrization of the deuteron wave func- 
tion [263,264] calculated from the Paris potential [136]. 



From the measurement of the momentum vector of the spectator pro- 
ton one can infer the momentum vector of the neutron p„ = — p^p 
and hence calculate the excess energy Q for each event, provided that 
the beam momentum is known. As an example, a distribution of the 
excess energy in the quasi-free pn pni]' reaction is presented in fig- 
ure 16.28b [20]. Due to the large centre-of-mass velocity (/? ~ 0.75) with 
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respect to the colliding nucleons, a few MeV wide spectrum of the neu- 
tron kinetic energy inside a deuteron (fig. 16.27b ) is broadened by more 
than a factor of thirty. 
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Figure 6.28: (a) Distribution of the excess energy Q for the pnq' system originating from 
the reaction pd psppnrj' calculated with a proton beam momentum of 3.350 GeV/c and the 
neutron momentum smeared out according to the Fermi distribution shown in figure IS?77h . 
(b) Spectrum of the off-shell mass of the interacting neutron, as calculated under the as- 
sumption of the impulse approximation. 

Therefore, especially in the case of near-threshold measurements, where 
the cross section grows rapidly with increasing excess energy (see e.g. 
fig. 14. 6|) . the total centre-of-mass energy ^/s has to be determined on an 
event-by-event level. For this purpose, the spectator protons are usually 
registered by means of silicon pad- or /x-strip detectors [20,158,265] which 
allow to determine their kinetic energy (T^p) and polar emission angle 
{$). Thus, it is useful to express the total energy squared s as function 
of these variables: 

s = |Pp + P„p = So - 2 Tsp (ma + Ep) + 2 p^ ^J^% + 2mp^,p cos(^) (6.5) 

with So denoting the squared centre-of-mass energy, assuming a vanishing 
Fermi motion. Measuring both the energy and the emission angle of 
the spectator protons it is possible to study the energy dependence of a 
meson production cross section from data taken at only one fixed beam 
momentum. 

It must be noted, however, that in the framework of the impulse ap- 
proximation, illustrated in figure 16.261 the measured spectator proton is 
a physical particle, yet the reacting neutron is off its mass shell, where 
the explicit expression for its four-momentum vector P„, in the rest frame 
of the deuteron, reads: 




P„ = (mrf - mp - Tsp, -Psp), 



(6.6) 
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with Tsp and p^^ denoting the kinetic energy and the momentum vector 
of a spectator proton, respectively. The mass spectrum of the inter- 
acting neutron (m^ = |P„P) resulting from the distribution of Fermi 
momentum is shown in figure IHUHb- It can be seen that the maximum 
of this spectrum differs only by about 3MeV/c^ from the free neutron 
mass (m„ = 939.57 MeV/c^), however on the average it is off by about 
9MeV/c^. In the framework of the discussed approximation, the struck 
neutron is never on its mass shell and the minimum deviation from the 
real mass occurs for vanishing Fermi-momentum and - as can be inferred 
from equation ()6.6|) - is equal to the binding energy = m^ — m„ — mp. 
Measurements performed at the CELSIUS and TRIUMF accelerators 
for the pp ppi] [121] and pp dn^ [266] reactions, respectively, have 
shown that within the statistical errors there is no difference between the 
total cross section of the free and quasi-free processes. This conjecture 
is confirmed also by the theoretical investigations [267]. 



Figure 6.29: (a) Total cross 
sections for the pp —> pprj reaction as 
a function of excess energy measured 
for the scattering of protons in vac- 
uum (open symbols) [10, 18, 50-53] 
and inside a deuteron (filled sym- 
bols) [121]. 
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In figure 16.291 the production of the 1] meson in free proton-proton 
collisions is compared to the production inside a deuteron and in the 
overlapping regions the data agree within the statistical errors. These 
observations allow to anticipate that indeed the assumption of the iden- 
tity for the transition matrix element for the meson production off free 
and quasi-free nucleons bound in the deuteron is correct, at least on the 
few per cent level. In case of the meson production off the deuteron, one 
can also justify the assumption of the quasi-free scattering with a geomet- 
rical argument, since the average distance between the proton and the 
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neutron is in the order of ^ 3fm. Of course, the other nucleon may scat- 
ter the incoming proton and the outgoing meson. However, these nuclear 
phenomena are rather of minor importance in case of the production on 
the neutron bound in the deuteron, but should be taken into account 
for derivations of total cross sections from experimental data. Indeed, 
calculating an influence on the cross section from a possible rescattering 
of the cu meson in the final state of the pd — > dujpsp reaction, authors 
of reference [269] found that for the spectator-proton momentum below 
100-150 MeV/c the rescattering effects play a minor role. The reduction 
of the beam flux on a neutron, due to the presence of the proton, referred 
to as a shadow effect, decreases for example the total cross section by 
about 4.5 % [270] for the r^-meson production. Similarly, the reduction of 
the total cross section due to the reabsorption of the outgoing 77 meson 
on the spectator proton was found to be only about 3% [270]. The ap- 
praisals were performed according to a formula derived in reference [271] 
which shows that the cross section for the deuteron reduces by a factor 
of: 

R=l- af^^ < r-2 > /4n (6.7) 

compared to the free nucleon cross sections. Here cr^'J^' denotes the r]N 
inelastic cross section and < r~^ > stands for the average of the in- 
verse square nucleon separation in the deuteron taking the nucleon size 
into account [270]. The latter effect for the production of mesons like 
(tt, cu, r)', (f)) is expected to be much smaller, since the s-wave interaction 
of the ?7-meson with nucleons is by far stronger than for any of the men- 
tioned ones. 

6.5 Test measurement of the pn pnr] reaction 

Nothing is so trivial as treating serious subject 
in a trivial manner [272]. 

Erasmus of Rotterdam 

As a general commissioning of the extended COSY-11 facility to investi- 
gate quasi-free pn — > pnX reactions, we have performed a measurement 
of the pn pnr] process at a beam momentum of 2.075 GeV/c [21]. 
The experiment, carried out in June 2002, had been preceded by the in- 
stallation of a spectator [158,273,274] and neutron detectors [275,276], 
and by a series of thorough simulations performed in order to determine 



* The matter radius of the deuteron amounts to 2 fm [268] . 
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the best conditions for measuring quasi-free pn pnt] and pn pnrj' 
reactions [20,28,230]. Figure ESDI presents the COSY-11 detection fa- 
cihty with superimposed tracks of protons and neutron originating from 
the quasi-free pn —* pnX reaction induced by a proton beam [62] im- 
pinging on a deuteron target [240]. The identification of the pn pnrj 
reaction is based on the measurement of the four-momentum vectors of 
the outgoing nucleons and the rj meson is identified via the missing mass 
technique. The slow proton stopped in the first layer of the position sen- 
sitive silicon detector (Sispec) is, in the analysis, considered as a spectator 
non-interacting with the bombarding particle and moving with the Fermi 
momentum as possessed at the moment of the collision. 



EXIT WINDOW 



Figure 6.30: Schematic view of the 
COSY-11 detection setup [61]. Only 
detectors needed for the measurements 
of the reactions pd — > psppn-qij]') and 
pd pprisp are shown. 
Dl, D2 denote the drift chambers; SI, 
S2, S3, S4 and V the scintillation detec- 
tors; N the neutron detector and Si,„on 
and Sispec [158] silicon strip detectors to 
detect elastically scattered and spectator 
protons, respectively. 




As we described in subsection 16.41 from the measurement of the mo- 
mentum vector of the spectator proton Pj,p one can infer the momentum 
vector of the struck neutron p„ = — p^p at the time of the reaction and 
hence calculate the total energy of the colliding nucleons for each event. 
In figures 16.31b and 16.31b the measured and expected distribution of 
the kinetic energy of the spectator proton is presented. Though still a 
very rough energy calibration of the detector units was performed, one 
recognizes a substantial similarity in the shape of both distributions. 
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Figure 6.31: Distributions of the kinetic energy of the spectator protons, 
(a) Experiment. (b) Monte-Carlo simulations taking into account the acceptance of 
the COSY- 11 detection system and an analytical parametrization of the deuteron wave 
function [263] calculated from the PARIS potential [136]. 

Figure 16.321 shows spectra of the excess energy in respect to the pnt] 
system as obtained in the experiment ()6.32b ) and the simulation ()6.32b ) 
for the pn pnr] reaction. The remarkable difference between the distri- 
butions comes from the fact that in reality additionally to the pn pnt] 
reaction also the multi-pion production is registered. The rj and multi- 
pion production cannot be distinguished from each other on the event- 
by-event basis by means of the missing mass technique. 
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Figure 6.32: Distributions of the excess energy Q for the quasi-free pn pnX reaction, 
determined with respect to the pur] threshold, (a) Experiment, (b) Simulation. 

However, we can determine the number of the registered pn pnrj 
reactions from the multi-pion background comparing the missing mass 
distributions for Q values larger and smaller than zero. Knowing that 
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negative values of Q can only be assigned to the multi-pion events we 
can derive the shape of the missing mass distribution corresponding to 
these events. This is shown as the dashed line in figure 16.331 



150- 



c 

o 
o 



100- 



50 - 



pn --> pn X 
above threshold of pn --> pn r| 
below threshold of pn -> pn r\ 




0.1 0.2 0.3 0.4 0.5 

mx [ GeV / c' ] 



75 



I 50 

O 
o 



25 



pn-->pnr| 

SIMULATION 
EXPERIMENT 




I I I I i_!!i I I I I I I I I L 



0.1 0.2 0.3 0.4 0.5 

mx [ GeV/c' ] 



b) 



Figure 6.33: Missing mass spectra as obtained during the June'02 run: 

a) Event distribution for Q < ( dashed line) and for Q > (solid line). 

b) Solid line represents the difference between number of events above and below threshold 
for the pn pnr] reaction, and the dashed line corresponds to the Monte-Carlo simulation. 

A thorough evaluation of the background is in progress, however, 
rough comparison of events for positive and negative Q yields the promis- 
ing results with a clear signal from the pn pnt] reactions, as can be 
deduced by inspection of figures IH75^ and 16.33b . 



6.6 Double quasi-free production 

There is need of a method for investigating 
the truth about things [207]. 

Rene Descartes 

In article [9] we described a method of measuring the close-to-threshold 
meson production in neutron-neutron collisions, where the momenta of 
the colliding neutrons could be determined with the accuracy obtainable 
for the proton-proton reaction. The technique is based on the double 
quasi-free nn — > nnX^ reaction, where deuterons are used as a source of 
neutrons. 

Close-to-threshold meson production in proton-neutron collisions were 
investigated by means of a technique based on a quasi-free scattering of 
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the proton off the neutron bound in the deuteron. Thin windowless 
internal deuterium cluster targets (~ 10^^ atoms/cm^) make a detection 
of an undisturbed spectator proton and a precise determination of the 
reacting neutron momentum - and hence of the excess energy - possible. 
Pioneering experiments of the vr'' meson creation in the proton-neutron 
reaction with the simultaneous tagging of the spectator proton resulted 
in a resolution of the excess energy equal to cr(Q) = 1.8 MeV [158]. 
Similar studies including the production of heavier mesons are carried 
out at COSY-11 and ANKE facilities [28,278], and we intend to continue 
them at the forthcoming facility WASA at COSY. 
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Figure 6.34: Schematic depiction of the double quasi-free nn —> nnX reaction. During the 
collisions of deuterons (left hand side of the figure, with the total momentum (solid arrow) 
resulting from the sum of the beam momentum (dotted arrow) plus the Fermi momentum 
(short arrow)) a double quasi- free neutron- neutron reaction may lead to the creation of 
mesons (small gray circle). The spectator protons (black circles) leave the reaction region 
with their initial momentum plus the Fermi momentum, which they possessed at the moment 
of the reaction. Neutrons are plotted as open circles. Due to the large relative momenta 
between spectators and the outgoing neutrons (~ 1 GeV/c close-to-threshold for the rj meson 
production) a distortion of the nnX system by the accompanied protons can be neglected. 

Experimental investigations of the close-to-threshold production in 
neutron- neutron collisions, however, have not yet been carried out. A 
realisation of such studies - which are characterised by typical cross sec- 
tions of < yub - with high quality neutron beams bombarding a deu- 
terium target is not feasible due to the low neutron beam intensities 
forcing to use liquid or solid deuterium targets which make the precise 
determination of the momentum of the spectator proton impossible. In 
this section a unique possibility of the precise measurement of the close- 
to-threshold meson production in neutron-neutron collisions is pointed 
out. The technique is based on the double quasi- free interaction of neu- 
trons originating from colliding deuterons as depicted in Figure 16.341 
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Utilizing this method, a precision of ~1 MeV can be obtained for deter- 
mining the excess energy, since it depends only on the accuracy of the 
momentum or angle reconstruction for the registered spectator protons. 
At present cooled deuteron beams - available at the facilities CELSIUS 
and COSY - give the possibility of using this method for the studies of 
neutron-neutron scattering. Moreover, the usage of a stored beam cir- 
culating through an internal cluster target permits the study with high 
luminosities ( 10^^cm~^s~^) in spite of very low target densities. These 
investigations can be realized much more effectively once the WASA de- 
tector is installed at COSY. In the double quasi-free interaction, due to 
the small binding energy of the deuteron {Eb = 2.2 MeV), the collid- 
ing neutrons may be approximately treated as free particles in the sense 
that the matrix element for quasi-free meson production from bound 
neutrons is identical to that for the free nn —>■ NNX reaction at the 
same excess energy available in the NNX system. The measurements at 
CELSIUS [121, 197] and TRIUMF [266, 277] have proven that the off- 
shellness of the reacting neutron can be neglected and that the spectator 
proton influences the interaction only in terms of the associated Fermi 
motion [266]. The registration of both spectator protons will allow for 
a precise determination of the excess energy. A possible internal tar- 
get facility based on the COSY-11 setup [61] is presented in figure 
The energy and the emission angle of the "slow" spectator can be mea- 
sured by an appropriately segmented silicon detector, whereas the mo- 
mentum of the "fast" spectator proton can be analysed by the magnetic 
spectrometer. By means of the detection system shown in figure I6.35t 
a resolution of the excess energy of 2 MeV can be achieved for excess 
energies lower than 30 MeV as demonstrated in reference [28]. The dou- 
ble quasi-free nn — *• nnX^ reaction can be identified by the registration 
of both outgoing neutrons. For example, in order to measure the pro- 
duction of the 1] meson, a seven meter distance for the time-of-flight 
measurement would be enough to obtain 8 MeV (FWHM) missing mass 
resolution [28] with a calorimeter segmented by 10 cm * 10 cm and pro- 
viding a 0.5 ns (ex) time resolution, which was obtained in test runs using 
a scintillator/lead sandwich type of detector. The suggested meson pro- 
duction via a double quasi-free neutron-neutron reaction with precisions 
achievable for the proton-proton and proton-neutron reactions, opens 
the possibility of studying for example the charge symmetry breaking 
by comparing cross sections for the pp pprj and nn nnr] reactions, 
similarly to investigations performed via the vr-deuteron reactions [185]. 
The Dalitz-plot analysis of the nn nn Meson would allow for the study 



6. Experiment 



131 



of the neutron-neutron and neutron-Meson [279] scattering lengths, the 
first being still not well estabhshed [64] and the second being unknown. 
In principle when studying the meson production in proton-proton and 

in proton-neutron collisions one has access to all possible isospin combi- 
nations, which can be derived after the correction for the electromagnetic 
interaction. 




Figure 6.35: Schematic view of the extended COSY-11 detection setup [61]. Only detec- 
tors needed for the measurements of the dd — » nnpspPspX reactions are shown. 
D1,D2 denote the drift chambers used for the track reconstruction of the fast spectator pro- 
ton; S1,S3 and V1,V2 are the scintillation detectors used as time-of-flight and veto counters, 
respectively, N1,N2 the neutron detectors, and Si [158] the silicon strip detectors. 

Exceptionally, close-to-threshold meson production via the neutron- 
neutron scattering represents a pure T = 1 isospin channel without ac- 
companying Coulomb interaction and consequently no need for its cor- 
rection. Investigations of neutron-neutron scattering allow also for the 
production of K^K~ pairs in a system with only two charged particles in 
the final state {nn nnK~^K~), simplifying the theoretical calculations 
drastically, which in case of the pp — > ppK^ K~ are not feasible due to 
the difficulty of treating the electromagnetic forces in the system of four 
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charged particles [280]. 

At present the COSY synchrotron can accelerate deuterons up to 
3.5 GeV/c [62] which, utihzing the Fermi momentum, allows for the tt 

and 1] meson production in the nn nnX^ reaction. To investigate the 
neutron-neutron interaction with the production of heavier mesons like 
a;, or 0, a deuteron beam of ~ 7 GeV/c would be required. 
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Using the stochastically cooled proton beam of the cooler synchrotron 
COSY and the COSY-11 facihty we performed measurements of the 
PP ~^ PPV PP ~^ PPV' reaction close to the kinematical threshold. 
We have succeeded in establishing the energy dependence of the total 
cross section for both reactions and the probability density of the phase 
space population for the pp ppr] reaction. At present we are pursu- 
ing the study by elaborating differential cross sections of the pp — > pprj' 
reaction based on the high statistics data from a recently performed ex- 
periment [179]. Measuring four-momentum vectors of protons in the 
entrance and exit channels we determined complete kinematics for each 
registered event and hence derived all in principle possible empirical in- 
formation about the studied processes. Thus, the interpretation of the 
results is limited only by the experimental accuracy and ambiguities of 
the theoretical approaches. 

Inferences of the interaction taking place within the ppr] and pprj' sys- 
tems were based on the comparison between the experimental results 
and the predictions obtained assuming that all kinematically permitted 
momentum combinations of the outgoing particles arc equally proba- 
ble. Indeed, we have discovered statistically incontestable deformations 
of both the excitation function and the phase space abundance. The 
discrepancies can plausibly be assigned to the influence of the hadronic 
interaction occuring among ejectiles since the range of excess energy was 
chosen such that other possible explanations - like variations resulting 
from the dynamics of the primary production or contributions from larger 
than zero angular momenta - can be excluded. 

A unique precision achieved due to the low emittance of the cooled 
beam and the high resolution mass spectrometer enabled us to distiguish 
effects resulting from the ?7-proton interaction from effects caused by 
the almost two orders of magnitude stronger proton-proton force. This 
conclusion was drawn from the phenomenological analysis based on the 
assumption that the primary production can be separated from the flnal 
state interaction and also that the three-body final state can be treated 
as an incoherent system of pairwise interacting objects. We must admit 
that these simplifications enabled us to estimate the effects from the 
proton-77 and proton-77' interaction only qualitatively. Nevertheless, a 
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remarkable difference l:)etween the shape of the excitation functions of 
the pp ppq and pp ppr]' reactions allowed us to conclude that the 
interaction between the r)' meson and the proton is significantly weaker 
than the analogous interaction between the 77 meson and the proton. To 
raise the confidence in the concluded distinction between proton-proton 
and meson-proton interaction we compared the shapes of the excitation 
functions of the reactions pp — > ppr) and pp — > pprj' with the one of 
pp ppn^. We could plausibly assert that the differences in shapes are 
due to the proton-r/, or correspondingly, to the proton-?]' interaction since 
the well known proton- tt" interaction is too weak to manifest itself within 
the current accuracy. Therefore the excitation function for the pp — > ppir^ 
reaction reflects the influence of the interaction between protons only. 
This comparison confirmed - model independently - that the hadronic 
interaction between the f]' meson and the proton is indeed much weaker 
than the 77-proton one. This is the first ever empirical appraisal of this 
hitherto entirely unknown interaction occuring between the rj' meson and 
the proton. 

We minimized a systematic bias of the data performing a multidi- 
mensional acceptance correction and the bin-by-bin subtraction of the 
multi-pion background. In particular, we elaborated invariant mass dis- 
tributions of two-particle subsystems of the ppr] final state independently 
of the reaction model used in the simulations. We must confess, how- 
ever, that with the measurement method used and within the achieved 
statistics it was not possible to determine a background-free occupation 
density over the Dahtz-plot unless a small number of bins was used, which 
in turn caused drastic losses in the experimental accuracy. Therefore as 
a next step of our investigation we proposed a pertinent experiment at a 
WASA@COSY facflity [181]. 

The experimental technique employed and the developed method of 
analysis allowed us to achieve an unique precision for the simultaneous 
determination of the absolute excess energy {a{Q) — 0.4 MeV), missing 
mass {a{m) ~ 0.3 MeV/c^), and momentum of the outgoing particles 
{(t{p) = 4 MeV/c). This unprecedented accuracy enabled us not only 
to find out that the density of the phase space population is strongly 
anisotropic but also to discern distortions originating from the proton- 
proton and proton-r^ interactions. An utterly satisfactory description 
of the deviations observed at low relative proton-proton momenta was 
achieved by taking into account the well known hadronic and electromag- 
netic interactions between the protons. On the contrary, the enhance- 
ment at low relative proton-77 momenta is a few times larger than ex- 
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pected under the assumption that the mutual interaction among particles 
in the ppr; system can be described by the incoherent pairwise interaction 
of proton-proton and proton-r] subsystems. At present many theoretical 
groups are seeking for an explanation of this effect [79,82,122,128,141]. 
Undoubtedly, this and other results presented in this treatise accelerated 
the development of the formalism of the description of the three-body 
system in the complex hadronic potential [123, 128] and led to signifi- 
cant progress in the understanding of the production mechanism on both 
hadronic [81,141,228] and quark-gluon levels [223]. 

Since the interaction between hadrons, their structure, and produc- 
tion dynamics are inseparably connected with each other, in this treatise 
we have discussed not only the hadronic interaction between the rj and 
T)' mesons and the protons but also presented our endeavour to under- 
stand the reaction dynamics and the structure of the studied mesons. 
The observed large difference of the total cross sections between the r] 
and rj' meson production shows that these mesons are created via differ- 
ent mechanisms, since comparable coupling constants are expected for 
both of them, at least in the SU(3) limit. It is well established that 
close-to-threshold the rj meson is produced predominantly through the 
excitation of the 5*11(1535) resonance, and hence the large difference in 
the observed cross sections suggests that the rj' primary production pro- 
cess is nonresonant. Different production mechanisms reflect differences 
in the structure of these mesons, however, at the present stage we cannot 
draw any quantitative conclusion. 

As far as the meson structure is concerned we focused on the r}' me- 
son demonstrating a method which will allow to investigate a gluonic 
component of its wave function. The study of this very interesting as- 
pect - connected to the search of matter built exclusively out of gluons - 
will be based on the comparison of the t]' meson production via the 
pp — > pprj' and pn pnr]' reaction. The measurement of the excitation 
function of the pp pprj' reaction is finished, and in order to accom- 
plish the measurement of the pn — > pnr)' reaction we have extended the 
COSY- 11 detection setup by a neutron and a spectator detectors. The 
neutron detector was designed and built by us in Cracow [275,276], and 
for the registration of the spectator protons we adapted [28, 273, 274] to 
the COSY-11 facihty a sihcon pad detector [158] which was previously 
used by the WASA/PROMICE collaboration at the CELSIUS accelera- 
tor. A successful measurement of the pn — * pnr] reaction has proven that 
it is possible to study quasi-free pn — > pnX reactions using the newly 
accomplished experimental setup. 
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An exploration of the isospin degrees of freedom in the study of the 
creation of mesons led us to the entirely novel idea of using quasi-free 
reactions for the investigations of meson production cross sections as a 
function of the virtuality of the interacting nucleons. We demonstrated 
also a method for measuring the close-to-threshold production of mesons 
in quasi-free neutron-neutron collisions. It is our great hope to conduct 
these investigations in the near future at the facihty WASA@COSY [181]. 

In parallel we began also to explore spin degrees of freedom. In par- 
ticular, we have demonstrated that by measuring spin observables it is 
possible to learn about the very details of the production dynamics of the 
77 meson. Predictions for the analysing power differ so much depending 
on the assumed mechanism that the determination of this quantity alone 
will allow to establish whether vector or pseudoscalar meson exchange 
dominates the production process. 

In this treatise we presented results of investigations aiming to deter- 
mine the interaction of r) and rj' mesons with nucleons, the structure of 
these mesons and the mechanism which governs their creation in the col- 
lision of nucleons. We have not arrived at the very aim yet. But we have 
gained a qualitative understanding of the reaction mechanism, estimated 
the hithertho entirely unknown Ty'-proton interaction, observed a signal 
from the 77-proton interaction, delivered a sample of high quality data on 
both total and differential cross sections, and developed novel methods 
of measurements which should lead to more quantitative statements in 
the near future. 

So even though my speculations pleased me very 
much, I believed that other persons had their own 
speculations which perhaps pleased them even 
more [281]. 

Rene Descartes 

The approbation of the public I consider as the 
greatest reward of my labours; but am determin'd 
to regard its judgment, whatever it be, as my best 
instruction [183]. 

David Hume 
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